rr. 


AD-A042  197 


KCM-WRE/YTO  SEATTLE  WASH  F/G  8/8 

ENVIRONMENTAL  PLANNING  FOR  THE  METROPOLITAN  AREA  CEDAR-GREEN  RI— ETC(U) 
DEC  74  DACW67-73— C-0022 


ADA042197 


+.  * 

Environmental  Management  for  the  Metropolitan  Area 
I Cedar-Green  River  Basins,  Washington 


r 


December  1974 


U.  S.  Army 
Corps  of  Engineers 
Seattle  District 


ENVIRONMENTAL  PLANNING  FOR  THE  METROPOLITAN  AREA 
CEDAR-GREEN  RIVER  BASINS,  WASHINGTON. 

.r  . 

, - — URBAN  DRAINAGE  STUDY  / / 


APPENDIX  B 

URBAN  STORM  DRAINAGE  SIMULATION  MODELS 


Representing 


Members 


Kenneth  Lowthian,  Vice-Chairman 
Robert  Gulino*,  Phil  Buswell 
Nancy  Rising*,  Chris  Smith* 
Arthur  Knutson 

Charles  Ede*,  Warren  Gonnason 
Pat  Nevins 
Mart  Kask 

George  Mi liman*,  Robert  Sloboden 
Robert  Meyer 
Ivan  Day* 

Gustav  Anderson 
James  Smith*,  George  Sherwin 
Jean  DeSpain,  Chairman 
Charles  Gibbs**,  Richard  Page 
Arthur  Dammkoehler 
Shirley  Farley 


City  of  Seattle  / - ■ / 

City  of  Seattle 

City  of  Bellevue,  Cedar  River  basin 
City  of  Kirkland,  Cedar  River  basin 
City  of  Renton,  Green  River  basin 
City  of  Auburn,  Green  River  basin 
Puget  Sound  Governmental  Conference 

Water  District  108  & Cedar-Green  basin  water  districts 
Ronald  Sewer  District  & Cedar  River  basin  sewer  districts 
Lakehaven  Sewer  District  & Green  River  basin  sewer  districts 
Rainier  Vista  Sewer  District  & Green  River  basin  sewer  districts 
Snohomish  County 
King  County 

Municipality  of  Metropolitan  Seattle 
Puget  Sound  Air  Pollution  Control  Agency 
Task  Force  for  Citizen  Participation 


ACCESSION  for 


NTlS 

DDC 

UNANNOUNCED 

JUSTIFICATION 


Ex  officio  members 


Washington  State  Department  of  Ecology 
U.S.  Environmental  Protection  Agency 
U.S.  Army  Corps  of  Engineers 


Robert  Stockman*  Fred  Hahn 
Robert  Burd 

Sydney  Steinborn*  Richard  Sellevold 


former  member 

♦former  member  and  former  chairman 


U.  S.  ARMY  CORPS  OF  ENGINEERS,  SEATTLE  DISTR 


Consulting  Engineers 

KRAMER,  CHIN  & MAYO  — WATER  RESOURCES  ENGINEERS 
YODER,  TROTTER,  ORLOB  & ASSOCIATES 


ACKNOWLEDGMENTS 


The  Corps  of  Engineers  expresses  appreciation  to  the  following  individuals  and 
groups  for  their  assistance  in  conducting  the  Urban  Drainage  Study: 


RIBCO  TASK  FORCE  FOR  CITIZEN  PARTICIPATION 
Shirley  Farley,  Chairwoman 
Brad  Collins,  Vice-Chairman 


Dale  Ashley 
Roy  Avent 
Jay  Becker 
Don  Beuthin 
Larry  Hall 
Jack  Locke 


Edwin  Hendrickson 
Kay  Johnson 
Henry  Keron 
Ted  Mathison 
Michele  Meith 
Bob  Porterfield 
Stew  Sargent 


Marion  Sherman 
Jerry  Tucker 
Jeanette  Veasey 
Forrest  Walls 
Bernice  White 
Ann  Widditsch 
Lon  Woodbury 


RIBCO  URBAN  DRAINAGE  TECHNICAL  REVIEW  SUBCOMMITTEE 


Roy  Avent 
Phil  Buswell 
Harvey  Duff 


MUNICIPALITY  OF  METROPOLITAN 


David  Glaze 
Robert  Gulino 
Dick  Hibbard 
George  Hsieh 


Michele  Meith 
Marvin  Seabrands 
George  Wanamaker 
Jim  Webster 


Nancy  Baggott 
Don  Benson 
Glen  Farris 


Pat  Levine 
Stephanie  Liff 
Theresa  Murphy 


Richard  Page 
Rod  Stroope 
Penny  Wilson 


KRAMER,  CHIN  AND  MAYO 
WATER  RESOURCE  ENGINEERS,  INC. 
YODER,  TROTTER,  ORLOB  AND  ASSOCIATES 


Ark  Chin 
Ron  Mayo 
Larry  Rugaard 
Gerald  Orlob 
Roger  Fry 
Dick  Warren 


President  - KCM 

Executive  Vice-President  - KCM 
General  Manager  - YTO 
President  - WRE 
Project  Manager 
Deputy  Project  Manager 


U.S.  ARMY  CORPS  OF  ENGINEERS 


Colonel  Raymond  Eineigl 
Sydney  Steinborn 
Richard  Sellevold 
Dwain  Hogan 
Walter  Farrar 


District  Engineer 
Chief,  Engineering  Division 
Chief,  Planning  Branch 
Chief,  Urban  Studies  Section 
Study  Manager 


I 


PREFACE 


This  report  is  an  appendix  to  the  Urban  Runoff  and  Basin  Drainage  Study. 
It  contains  program  documentation  for  the  urban  stormwater  computer  sim- 
ulation models  and  can  serve  as  a user's  manual.  The  computer  models 
described  have  been  used  to  simulate  the  quality  and  quantity  of  stormwater 
runoff  from  drainage  sub-basins  within  the  Green  and  Cedar  River  Basins 
(State  of  Washington  Water  Resource  Inventory  Areas  8 & 9)  and  to  develop 
alternatives  for  meeting  existing  and  long  range  drainage  needs. 

The  Urban  Runoff  and  Basin  Drainage  Study  is  part  of  an  environmental 
management  program  for  the  Green  and  Cedar  River  Basins  in  King  and  Snohomish 
Counties,  Washington  and  has  been  conducted  under  the  auspices  of  the  River 
Basin  Coordinating  Committee  (RIBCO).  Four  principal  studies  comprise  the 
RIBCO  Environmental  Management  Program:  Part  I - Water  Resources;  Part  II  - 

Urban  Drainage;  Part  III  - Water  Quality  and  Part  IV  - Solid  Waste. 

The  Urban  Runoff  and  Basin  Drainage  Report  presents  a comprehensive  plan 
for  meeting  the  existing  and  long  range  urban  drainage  needs  within  the  Green 
and  Cedar  River  Basins.  The  study  recommendations  address  drainage  facil- 
ities, capital  cost,  methods  of  financing  and  institutional  arrangements  for 
effective  drainage  management.  The  recommended  plans  are  conceptual  and  are 
intended  for  use  by  local  governments  as  a guide  in  the  future  planning  of 
drainage  systems. 

The  published  report  is  composed  of  the  following  documents: 

Technical  Report 

Appendix  A - Regional  Sub-Basin  Plans 

Volume  1 - Cedar  River  Basin 
Volume  2 - Green  River  Basin 

Appendix  B - Urban  Storm  Drainage  Simulation  Models 

Appendix  C - Storm  Water  Monitoring  Program 

This  report  is  submitted  in  compliance  with  the  terms  of  contract 
DACW67-73-C-0022  between  the  Seattle  District,  U.S.  Army  Corps  of  Engineers 
and  KCM-WRE/YTO. 
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CHAPTER  I 

DESCRIPTION  OF  MODELS 


INTRODUCTION 

HISTORY  AND  PURPOSE 

In  October  of  1971,  under  the  sponsorship  of  the  Environmental 

Protection  Agency,  a consortium  of  private  contractors — Metcalf  and 

> 

Eddy,  Inc.,  the  University  of  Florida  and  Water  Resources  Engineers,  Inc.-- 
presented  a comprehensive  mathematical  model  capable  of  representing  urban 
stormwater  runoff  and  combined  sewer  phenomena.  The  model  was  comprised 
of  a ser  individual  sections  or  blocks  which  could  be  used  either 

separa'  combination.  The  principal  elements  of  the  model  are 

the  r ansport,  storage  and  receiving  water  blocks. 

Subsequent  to  the  initial  development  the  various  blocks  of 
the  model  have  been  extensively  developed  and  modified  in  further 
applications.  In  a major  study  of  stormwater  drainage  for  the  City  of 
San  Francisco,  the  runoff  block  was  updated  to  include  revised  concepts 
of  gutter  flow  and  the  transport  model  was  rewritten  to  develop  a model 
capable  of  handling  looped  networks  typical  of  sewers,  plus  irregular 
piping  systems  and  surcharge.  These  models  have  been  developed  expressly 
to  simulate  urban  stormwater  runoff.  Their  purpose  is  to  provide  a means 
of  assessing  problems  in  an  urban  drainage  system  and  to  provide  a means 
of  examining  potential  solutions. 

The  models  developed  for  the  City  of  San  Francisco  are  those 
which  have  been  updated  and  modified,  particularly  with  respect  to  the 
number  of  quality  constituents,  for  use  in  the  Seattle  Metropolitan  area. 

It  is  a description  of  these  computer  models  which  is  presented  in  tnis 
document. 
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In  the  original  development  work  the  stormwater  models  were 
applied  to  test  watersheds  in  San  Francisco,  Philadelphia,  Cincinnati 
and  Washington,  D.C.  The  San  Francisco  and  Washington,  D.C.  areas  were 
later  more  comprehensively  analyzed  for  the  local  authorities.  Major 
studies  have  also  been  undertaken  in  the  cities  of  Chicago,  Illinois; 
Boston,  Massachusetts;  and  Hamburg  and  Landau  in  West  Germany. 

Individual  blocks  of  the  comprehensive  model  have  proved 
useful  in  more  limited  studies  in  other  areas.  For  example,  the  runoff 
block  has  been  used  to  design  overflow  facilities  in  the  City  of  Corte 
Madera,  California,  and  the  transport  block  in  various  flow  studies  for 
the  City  of  Albany,  the  East  Bay  Municipal  Utility  District  and  Marin 
County  Flood  Control  District,  all  in  California. 


The  models  described  herein  consist  of  four  primary  sections 
or  blocks  which  are  used  to  simulate  the  quality  and  quantity  of  flow 
from  a watershed,  through  a sewer,  and  then  to  display  the  results. 

Specifically,  the  blocks  described  in  this  report  are: 

1.  r.  ' ■'k  - A graphical  output  routine  for 
producing  plots  on  the  line  printer. 

2.  ak  y;’  a '<  - A set  of  computer  routines  which 
simulates  the  quality  and  quantity  of  surface  flow 
from  watersheds. 

3.  . >•; :r>t  ?■  ' -k  - A set  of  computer  routines  which 
dynamically  routes  flow  through  sewer  systems. 

4.  .”?•  7>v.  p‘  duality  b : fk  - A computer  routine  which 
uses  the  flow  simulation  of  the  Transport  Block  to 
route  quality  constituents  in  sewer  systems. 


Please  note  that  the  receiving  water  model  developed  in  the 
above-mentioned  EPA  work  is  not  included  in  this  report. 


GENERAL  DESCRIPTION  OF  MODEL  BLOCKS 


The  Display  Block  is  a set  of  four  subroutines  which  produce 
plots  of  specified  system  variables  on  the  line  printer.  The  items 
which  may  be  displayed  are  the  input  rainfall  hyetographs,  the  output 
flow  hydrographs  from  the  Runoff  and  Transport  Blocks,  and  the  output 
pol 1 utographs  for  quality  constituents. 

The  routines  included  in  the  Display  Block  have  been  programmed 
such  that  they  are  general  in  operation  and  employ  a Calcomp-like  calling 
sequence  in  their  formal  arguments.  Note  that  the  routines  of  the  Display 
Block  produce  continuous  lines  between  data  points,  a feature  not  usually 
included  in  printer  plots.  A more  complete  description  of  the  elements 
in  the  Display  Block  can  be  found  in  Reference  1 and  in  later  sections 
of  this  chapter. 

RUNJFF  BxjOCK 

The  Runoff  Block  consists  of  a set  of  computer  routines  and 
appropriate  data  which  will  simulate  the  rainfall-runoff  characteristics 
of  an  urban  area.  Certain  of  its  internally  programmed  quality  charac- 
teristics have  been  four d from  observations  in  the  Seattle  area  and  thus 
the  model  must  be  considered  somewhat  specific  at  the  present  time.  In 
the  model  flow  is  traced  from  the  onset  of  rainfall  to  the  watershed, 
through  overland  flow,  and  then  to  flow  in  the  gutters  and  minor  sewers. 
Water  quality  mass  emissions  are  generated  as  they  occur  on  the  watershed 
surface  and  in  the  catch  basins  which  comprise  the  total  system. 
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In  the  following  paragraphs  only  general  descriptions  are  given 
to  the  technical  methods  employed  in  the  Runoff  Block.  The  reader  is 
directed  to  Reference  1 and  to  Chapter  V for  a detailed  explanation  of 
the  model's  formulation  and  its  numerical  solution. 

‘u  ‘ Watershed 

It  is  important  to  recognize  that  any  model  is  an  idealization 
of  the  real  system  in  some  sense.  A model  is  intended  to  portray  certain 
important  characteristics  of  the  prototype  system  realistically  and  the 
degree  to  which  it  does  so  is  the  primary  measure  of  its  validity. 
Whenever  possible,  it  seems  desirable  to  represent  the  real  world  by  a 
model  which  conforms  to  our  best  knowledge  of  the  actual  processes 
taking  place.  In  a hydraulics  problem,  for  example,  we  will  worl  with 
the  familiar  entities  of  flow,  inertia,  continuity,  etc.,  albeit  in  ar. 
idealized  fashion.  Sometimes  the  real  system  is  too  complex  for  this 
kind  of  treatment  and  we  resort  to  statistical  measures  or  curve  fitting 
descriptions  of  the  phenomenon. 

The  processes  simulated  by  the  Runoff  Model  are  extremely 
complex  and  interrelated.  The  model  is  based  both  on  the  fundamental 
physics  of  hydromechanics  and  on  the  best  of  carefully  developed  sanitary 
engineering  technology.  Tests  of  the  model  show  that  it  fulfills  the 
criterion  set  forth  above--real i sti c portrayal  of  important  prototype 
beha vior. 


Hydraulic  flow  in  the  Runoff  Model  is  represented  by  the  kinematic 
wave  solution  for  flow  across  a plane  surface.  This  rather  simple  approach 
is  applied  to  a geometric  representation  of  an  urban  area  which  includes 
the  major  hydrologic  sub-units  of  surface  drainage.  In  the  model  each  of 
the  watershed  sub-units  is  treated  as  if  it  were  completely  independent  of 
all  others  in  the  system.  This  idealization  is  basic  to  the  structure  of 
the  Runoff  Model;  when  it  is  found  that  this  is  not  a good  approximation 


corrections  must  be  made  by  changing  watershed  boundaries  or  creating  new 
unit  watersheds.  Linkage  between  major  watersheds  is  usually  accomplished 
in  the  Transport  Model,  not  the  Runoff  Model. 


It  is  convenient  now  to  consider  in  detail  the  methods  employed 
in  modeling  flow  from  one  of  these  independent  watershed  units.  Clearly, 
the  problem  is  extremely  complex.  It  involved  flow  from  building  roofs, 
parking  lots,  lawns,  playfields,  streets,  gutters,  a sewer  network  and 
all  other  features  which  make  up  the  urban  landscape. 

We  shall  begin  our  analysis  by  considering  the  idealized 
watershed  indicated  in  Figure  1-1.  Figure  1-1  shows  the  three  types 
of  drainage  units  found  in  each  watershed.  These  are: 

1.  Pervious  areas  - parks,  lawns,  gardens,  etc. 

2.  Impervious  areas  - roofs,  paved  areas,  etc. 

3.  Minor  sewer  elements  - all  sewers,  gutters,  or 
similar  conveyance  paths  which  collect  the 
runoff  from  the  pervious  and  impervious  areas 

Each  of  tnese  is  treated  as  an  overland  flow  problem.  The  magnitude  of 
flow  is  determined  from  the  kinematic  wave  solution  mentioned  earlier. 

Rain  falls  directly  on  the  pervious  and  impervious  areas.  On  the  former, 
some  is  lost  by  infiltration  to  the  ground  water.  After  deducting  losses 
flow  moves  from  the  watersheds  to  the  minor  sewer  elements  where  the  total 
flow  is  accumulated.  Outflow  from  the  minor  sewer  elements  is  routed  to 
other  minor  sewer  elements  in  the  surface  system,  and  its  magnitude  stored 
by  the  computer  routine  for  later  use  in  the  Transport  Block. 
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PERVIOUS  AREA  ; IMPERVIOUS  AREA 


The  model  assumes  that  the  pervious  and  the  impervious  areas 
are  separate  within  the  watershed,  each  of  which  drains  into  the  adjacent 
street-gutter-sewer  system  which  constitutes  the  minor  sewer  element. 

Tne  minor  sewer  element  collects  the  watershed  flows  and  accounts  for 
the  lag  time  introduced  by  the  distance  from  the  watershed  to  the  major 
sewer  node. 

i .• ' ’ tu  Pvo*i  the  Watershed 

The  quality  of  the  surface  runoff  is  important  in  the  analysis 
of  urban  systems,  and  more  than  20  water  quality  constituents  have  been 
included  in  the  Runoff  Model.  The  contribution  from  each  subarea  for 
each  constituent  is  determined  by  combination  of  land  use  and  equivalent 
street-gutter  length.  It  is  the  responsibility  of  the  user  to  specify 
each  of  these  values  for  all  unit  watersheds  in  the  total  watershed. 

The  constitutive  relationships  for  all  quality  parameters  in 
tne  Runoff  Model  have  been  determined  from  observations  of  actual  sewer 
quality  in  the  Seattle  Metropolitan  area.  For  this  reason  the  quality 
model  must  be  considered  relatively  specific  to  the  Seattle  Metropolitan 
area  for  the  present  time. 

The  Runoff  Quality  Model  considers  two  sources  of  pollution; 
the  watersheds  themselves  and  street  catch  basins.  The  buildup  and  amount 
of  pollutant  in  each  are  made  functions  of  both  the  street  sweeping 
frequency  and  the  number  of  days  without  rain  prior  to  any  given  storm. 

At  the  outset  of  the  simulation  an  estimate  is  made  of  the  total  amount 
of  pollutant  contained  in  and  on  each  of  the  system  elements.  As  rainfall 
begins  pollutants  move  from  the  watersheds  and  catch  basins  to  the  minor 
sewer  areas  of  the  system.  A mass  balance  is  calculated  fu.  each  system 
element  at  each  time  step  and  estimates  are  made  for  quality  concentrations 
expected  over  the  entire  duration  of  a storm. 
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The  rate  of  pollutant  removal  is  proportional  to  the  rate  of 
rainfall  and  the  amount  of  pollutant  remaining  on  the  watershed  or  in 
the  catch  basin.  All  constituents  are  routed  conservatively,  and  output 
is  produced,  printed  or  plotted  at  any  combination  of  preselected  system  nodes. 

The  water  quality  model  is  an  integral  part  of  the  Runoff  Block 
but  it  is  possible  to  run  the  quantity  simulation  without  the  quality 
simul ation. 

Ikvut  D.:ia 


The  geometric  and  hydraulic  properties  of  the  unit  watersheds 
are  required  in  a form  consistent  with  the  idealization  of  the  area  and 
yet  compatible  with  an  ordinary  understanding  of  the  topography  and 
hydrology  of  urban  drainage  areas.  Details  of  the  computer  card  format 
requirements  are  contained  in  Chapter  II.  The  quantities  needed  for 
each  watershed  are: 

1.  Rainfall:  the  rainfall  hyetographs  pertinent 

to  the  watershed. 

2.  Linkage:  a watershed  ’'denti  fication  number  and 

the  sewer  node  to  which  it  drains. 

3.  Geometry:  the  size  and  length  of  each  unit 

watershed. 

4.  Hydraulic:  the  average  ground  slope,  hydraulic 

resistance  factors,  percent  impervious,  surface 
storage  (parking  lot  ponding,  surface  adherence 
to  trees,  etc.),  and  infiltration  coefficients. 

5.  Quality:  the  number  of  dry  days  prior  to  the 

onset  of  rainfall,  the  street  cleaning  frequency, 
and  the  equivalent  gutter  length  and  number  of 
catch  basins  for  each  unit  watershed. 


1-8 


t 


Certain  of  the  required  data  must  be  supplied  by  the  user.  In  other  cases 
the  program  will  select  default  values  if  the  user  elects  not  to  specify 
them.  The  reader  is  referred  to  Chapter  II  and  the  examples  of  Chapters 
III  and  IV  for  more  complete  descriptions  of  data  needs  and  problem  structure. 


The  Transport  Block  is  a set  of  computer  routines  and  data 
which  simulate  the  movement  of  stormwater  through  a sewer  system.  The 
model  uses  as  input  specific  hydrographs  such  as  those  produced  by  the 
Runoff  Block.  The  Transport  Block  can  be  used  to  simulate  quantity  only 
with  the  Transport  Model  or  both  quantity  and  quality  with  serial  operation 
of  the  Transport  and  Transport  Quality  Blocks. 

A fairly  detailed  description  of  the  elements  of  the  Transport 
Clock  is  presented  as  it  is  a completely  different  formulation  of  sewer 
flow  from  that  given  in  Reference  1.  Even  more  detail  on  tne  solution 
technique  can  be  found  in  Chapter  V. 


The  conceptual  representation  of  a storm  sewer  system  is 
illustrated  in  Figure  1-2.  It  employs  the  link-node  concept,  with  the 
basic  elements  listed  in  Table  1-1.  In  the  transport  model  the  property 
of  volume  is  ascribed  to  the  nodes.  The  solution  technique  applies  the 
continuity  equation,  with  storage,  to  each  node  in  the  system  and  computes 
the  hydraulic  head  at  the  manholes  or  junctions.  These  heads,  together 
with  the  friction  and  momentum  forces,  are  used  to  compute  the  flow  through 
the  pipes  or  links.  A summary  of  the  properties  ascribed  to  links  and  nodes 
in  the  model  are  listed  in  Table  1-2. 

By  applying  the  continuity  equation  with  storage  to  the  nodes 
and  then  using  the  resulting  heads  to  compute  the  flow  in  connecting 
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TABLE  1-1 


CLASSES  OF  ELEMENTS  INCLUDED  IN 
THE  TRANSPORT  MODEL 


Element  Class 

Types 

Conduits  or  Links 

Rectangular 

Circular 

Horseshoe 

Baskethandle 

Egqshape 

Trapezoidal 

Junctions  or  'lodes  (Manholes) 

Diversion  Structures 

Orifices 

Transverse  Weirs 
Sideflow  Heirs 

rump  Stations 

On-line  or  off-line  lift  station 

Storaqe  Rasins 

On-line  (enlarned  pipes  or  tunnels) 

Outfall  Structures 

Transverse  weir  with  tide  gate 
Transverse  weir  without  tide  qate 
Sideflow  weir  with  tide  qate 
Sideflow  weir  without  tide  gate 
Free  outfall  with  tide  gate 
Free  outfall  without  tide  qate 

TABLE  1-2 

PROPERTIES  OF 
THE 

THE  LINKS  AND  NODES  IN 
TRANSPORT  MODEL 

Properties  and  Constraints 

NODES 

Constraint 

;Q  * change  in  storage 

Properties 

Volume 
Surface  area 
Head 

L IUKS 

Constraint 

Qin  * Qout 

Propert ies 

Cross  sectional  area 
Length 
Rouahness 
Hydraulic  Radius 
Surface  Width 
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links,  the  backwater  effect  oin  be  propagated  through  the  system. 

Moreover,  the  use  of  this  technique  makes  the  order  of  computing  heads 
at  nodes  and  flows  of  links  unimportant.  All  water  transfers  info  or 
out  of  the  system  occur  at  nodes,  i.e.  weir  or  orifice  diversions  to  dry 
weather  interceptor  lines,  pumping  stations,  and  outfalls,  inflows  from 
surface  runoff,  or  inflows  diverted  from  another  part  of  the  system. 

Definition  of  System  Elements  and  Data  Requirements 

The  various  types  of  sewerage  system  elements  to  be  represented 
in  the  model  were  listed  in  Table  1-1.  In  this  section,  the  aeometric  - 
hydraulic  properties  of  these  elements  as  they  are  used  in  the  model 
wi  1 1 be  defined. 

Conduits  - Conduits  are  the  primary  conveyance  elements  in 
the  system  and  it  is  to  these  elements  that  the  dynamic  non-steady  flow 
equation  is  applied.  Thus,  for  each  pipe  in  the  system  the  following 
aeometric  - hydraulic  parameters  are  required: 

1 . Conduit  length, 

2.  Conduit  rouqhness, 

3.  Cross  sectional  area  as  a function  of  depth  of  flow, 

4.  Hydraulic  radius  as  a function  of  the  depth  of  flow,  and 

5.  Surface  width  as  a function  of  the  depth  of  flow. 

Conduit  length  and  roughness  must  be  supplied  as  input  data  for  the  model. 
However,  the  functional  relationships  between  parameters  3,  4 and  5 and 
the  depth  of  flow  are  retained  internally  in  the  computer  program  for  the 
followina  conduit  types: 

1.  Circular,  4.  Baskethandle, 

2.  Rectangular,  5.  Eggshape,  and 

3.  Horseshoe,  6.  Trapezoidal. 
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For  each  of  the  parameters  indicated,  an  array  of  values  is  stored  for 
each  Dipe  type  which  oives  the  ratio  of  the  value  of  the  parameter  at  a 
given  depth  to  its  maximum  value  for  the  pipe.  While  the  arrays  contain 
values  for  rectangular  conduits,  it  is  more  efficient  to  compute  the 
parameters  for  both  rectangular  and  trapezoidal  sections.  The  maximum 
values  of  these  parameters  are  the  cross-sectional  area  of  the  pipe, 
width  at  the  springline,  and  height  of  the  pipe.  The  idea  of  storing 
arrays  of  fractions  was  originally  conceived  in  the  development  of  tne 
Florida  Transport  Model,  and  it  is  an  extremely  efficient  way  to  maintain 
the  geometric  - hydraulic  characteristics  of  a large  number  of  pipes  of 
different  types  and  sizes. 
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Junati  > ns  - The  junctions,  or  nodes,  in  the  system  represent 
the  on-line,  or  in-pipe,  storage  capacity  of  the  system.  In  addition 
the  hydraulic  heads  at  the  nodes  are  used  to  drive  the  flow  throuah  the 
connecting  nipes  or  links.  Table  1-2  showed  the  properties  associated 
with  a node  to  be:  volume,  head,  and  surface  area.  Referring  again  to 

the  conceptual  representation  of  the  node  illustrated  in  Figure  1-2, 
it  is  seen  that  the  properties  of  the  node  can  be  defined  in  the  following 
manner. 


1.  For  each  conduit  connected  to  the  node  under  consideration 
(call  it  node  j)  the  water  depth  at  the  midpoint  of  the 
conduit  is  computed  as  the  average  of  the  denths  at  the 
ends  of  the  conduit  i.e.  the  average  of  the  depths  at 

the  connecting  nodes. 

2.  The  surface  widths  at  the  midpoint  and  at  the  end  connecting 
to  node  j is  then  found  using  the  depths  of  flow  at  these 
points  and  the  appropriate  lookun  table  as  described  above 
under  Conduits.  The  free  surface  area  of  the  pine  is  then 
computed  as  the  product  of  the  average  of  the  widths  at  the 
midpoint  and  the  end  of  the  conduit  times  the  half-length 

of  the  conduit.  The  surface  area  AS ( j ) of  node  j is  computed 
as  the  sum  of  the  surface  areas  of  the  half-pioes  connected 
to  the  node. 
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3.  Assuminn  that  any  excess  inflow  into  a node  is  distributed 
uniformly  over  the  water  surfaces  of  the  node,  the  change 
in  head  at  node  j is  computed  as  AH  = (IQ). /AS.  so  that 
the  new  head  at  the  node  is  expressed  as  J J 

Hj(t+At)  = HJ(t)  + [lQ(t+At)]j/AS.j(t)  (1-1) 

4.  The  volume  of  water  in  storage  at  the  node  is  never  used 

for  any  computation.  It  is  saved  only  as  a parameter  of 
interest.  At  any  time  the  chanqe  in  volume  at  node  j is 
comDuted  simply  as  AV^  = (£Q)j  so  the  v0^ume  of  water 

contained  at  a node  is  expressed  as 

Vj  (t+At)  = Vj(t)  + [2TQ(t+At)]j  (1-2) 

Of  course  when  all  pioes  are  full  'J . takes  on  the  maximum 
storaqe  value  of  the  node  and  [£Q(t+At)]j  becomes  flood 
water  which  flows  out  of  the  flooded  junction. 

P-rereim  ires  - Diversion  structures  are  used  for  one  of 

two  purposes: 

1.  To  divert  sanitary  sewerage  out  of  the  storm  drainage 
system,  and 

2.  To  relieve  the  stormwater  load  on  sanitary  interceptors. 

These  diversions  always  occur  at  a node  with  the  transfer  beino  made 
between  two  nodes.*  Descriptions  of  the  hydraulic  characteristics  of 
these  diversion  structures  are  given  below. 

The  purpose  of  the  orifice  is  to  divert  sanitary  sewage  out  of 
the  stormwater  system  during  dry  weather  periods  and  to  restrict  the  entry 
of  stormwater  into  the  sanitary  interceptors  durino  periods  of  runoff. 

The  orifice  may  divert  the  flow  to  another  pipe,  a pumping  station,  or  an 
off-line  storage  tank.  Figure  I-3A  shows  a schematic  of  two  types  of 
orifice  diversion  structures.  The  conceptual  representation  as  used 
the  model  is  illustrated  in  Figure  I-3B.  The  orifice  is  completely 
characterized  by  two  parameters: 

♦There  is  no  pipe  connection  between  the  two  nodes. 
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Schematic  of  Two  Types  of  Orifice  Diversions 


PLAN  VIEW 

Conceptual  representation  of  an  Orifice  Diversion 


FIGURE  1-3 

REPRESENTATION  OF  ORIFICE  DIVERSIONS 


1.  Its  cross-sectional  area  A,  and 

2.  A discharge  coefficient,  CQ. 

The  value  of  the  discharge  coefficient,  which  depends  upon  the  type  of 
opening  and  the  length  of  the  orifice  tube,  is  pre-specified  as  input 
data  along  with  the  cross-sectional  area  of  the  orifice.  The  discharge 
through  the  orifice  is  computed  with  the  standard  orifice  equation,  i.e. 

Q0  = CQ  A >2?h  (1-3) 

where  g is  the  oravi tational  constant  and  h is  the  head  on  the  orifice 
i.e.  the  minimum  value  of  (H^-H^)  and  Y^  (See  Figure  1-3).  This  expression 
is  used  even  when  the  orifice  is  not  submerged.  The  resulting  error, 
however,  is  not  considered  significant  since  these  flows  are  small  when 
the  orifice  is  unsubnerged. 

Weir  diversion  in  the  storm  sewer  system  provide  relief  to 
the  sanitary  sewers  durino  periods  of  storm  runoff.  Two  types  of  weir 
diversions  are  treated  in  the  simulation  model: 

1.  Transverse  weirs,  and 

2.  Sideflow  weirs. 

The  only  difference  computationally  is  that  the  weir  coefficients  are 
different  for  the  two  cases. 

A weir  diversion  schematic  and  its  conceptual  representation 

are  illustrated  in  Figure  1-4.  When  the  water  surface  elevation  at  Mode  1 

rises  above  Y , the  diversion  of  flow  to  Mode  2 begins.  The  flow  diversion 
c 

is  governed  by  the  weir  equation 
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where  h is  the  head  on  the  weir,  V the  approach  velocity,  g the  gravi tational 
constant,  L the  weir  length  and  Cw  the  discharge  coefficient.  The  exponent 
a is  set  at  3/2  for  transverse  weirs  and  5/3  for  sideflow  weirs.  In  the 
present  version  of  the  computer  program,  V is  not  used  in  the  computation 
of  the  weir  discharge.  The  value  is  set  to  zero  just  prior  to  the 
computation  of  discharge  by  Equation  1-4.  However,  it  can  be  added  with 
a minimum  amount  of  effort  if  the  need  to  include  it  arises. 

The  discharge  coefficient  Cw  is  prespecified  for  transverse 
weirs  as  input  data  along  with  the  weir  length.  For  sideflow  weirs  this 
input  coefficient  should  be  modified  internally  during  program  execution 
to  reflect  the  effect  of  the  approach  velocity  on  the  rate  of  flow 
diversion.  The  computer  program  contains  a provision  for  this  computation; 
however,  the  present  expression  simply  equates  the  discharge  coefficient 
to  the  input  value.  The  correct  expression  can  be  formulated,  however, 
and  incorporated  into  the  program. 

The  head  over  the  weir,  h is  computed  as  the  minimum  value  of 

(Y.J-Y)  and  jH-j  - H^l . If  is  negative,  a backflow  condition 

exists  and  the  flow  is  from  Node  2 to  Node  1.  In  this  case,  the  transverse 

and  sideflow  weirs  are  both  treated  as  a transverse  weir  with  a = 3/2 

and  C = input  value, 
w 

If  the  weir  is  submerged,  the  Villamonte  weir  correction  is 
anolied  to  the  discharge  Qw  computed  by  Equation  1-4,  i.e. 

0.385 

(1-5) 

If  the  weir  is  surcharged,  it  behaves  like  an  orifice,  and  the  discharge 
is  expressed  as 

Q = C [L  (Y  - Y )]  *4qh  + V2  (1-6) 

vw  w L v max  c/J  3 
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where  (Yrdx  - Y ) is  the  distance  between  the  weir  crest  and  the  top  of 
trie  pipe.  The  only  difference  between  this  expression  and  equation  4 is 
that  the  velocity  head  of  the  approach  flow  has  been  included  for  the  case 
of  the  submerged  weir.  The  discharge  coefficient  for  the  surcharged 
case  is  an  input  value.  However,  if  Equations  1-4  and  1-6  are  equated  at 

Y = Ymax,  and  the  effect  of  the  velocity  heads  is  neglected  (i.e.  set 

O 

V = 0 in  the  two  expressions)  it  can  be  shown  that 


C‘ 

w 


w 

t'Za 


_w 
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The  exDonent  a is  taken  as  3/2  in  this  case.  It  is  recommended  that 


Equation  1-7  be  used  to  define  C in  model  applications. 

w 


Purry  St-.v  ns  - A pump  station  is  conceptually  represented  as 
an  off-line  storaqe  node  (the  wet-well)  from  which  the  contents  are  pumped 
to  another  node  in  the  system  according  to  a programmed  rule  curve. 

Inflows  to  the  storage  node  may  be  diverted  from  the  main  sewer  system 
through  an  orifice,  a weir,  or  a pipe  with  a free  discharge  condition 
at  the  storane  node.  The  Dumping  rule  curve  is  based  on  the  volume  of 
water  in  the  storaqe  node.  It  operates  as  follows: 


1.  Up  to  three  wet-well  volumes  are  prespecified  as  input 

data  for  each  oumn  station:  V-j  < < V^,  where  is  the 

maximum  capacity  of  the  wet  well. 

2.  Three  pumping  rates  as  also  Drespecified  as  input  data  for 
each  station: 

R-j  for  volume  in  wet  well  < V-j 

R„  for  V-j  < volume  in  v/et  well  < V,, 

R3  for  < volume  in  wet  well  < 

3.  A mass  balance  of  outflows  as  specified  by  the  rule  curve 
in  item  2 and  inflows  is  performed  in  the  wet  well  during 
the  model  simulation  period. 
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4.  If  the  wet  well  goes  dry  the  pump  rate  is  reduced  below 
rate  R-j  until  it  just  equals  the  inflow  rate.  When  the 

inflow  rate  again  equals  or  exceeds  R^,  the  pumping  rate 

goes  back  to  operating  on  the  rule  curve. 

5.  If  V ^ is  exceeded  in  the  wet  well,  the  inflow  to  the 

storage  node  is  reduced  until  it  does  not  exceed  the 
maximum  Dumped  flow.  When  the  inflow  falls  below  the 
maximum  pumped  flow,  the  inflow  "gates"  are  opened  again. 


Storage  Basins  - Table  1-1  shows  that  the  transport  model  includes 
on-line  storage.  On-line  storage  is  that  storage  that  is  available  in 
the  system  by  virtue  of  the  fact  that  the  pipe  and  conduits  which  make 
up  this  system  store  water.  On-line  storage  was  described  above  under 

Junctions. 


Outfall  Structures  - The  basic  outfall  structures  that  can  be 
simulated  by  the  transDort  model  are  weir  outfalls  and  free  outfalls, 
both  types  with  or  without  tide  gates.  The  hydraulic  - geometric 
characteristics  of  the  free  outfall  weir  have  been  described  above  under 
Diversion  Structures . It  should  be  noted  that  for  the  free  outfall  weir 
(no  tide  gate ) it  is  tacitly  assumed  that  the  condition  of  backflow  across 
the  weir  will  never  occur'.  If  a backflow  situation  does  occur,  the  flow 
over  the  weir  is  automatical ly  set  to  zero.  This  procedure  in  effect 
prevents  receiving  waters  from  flowing  into  the  sewer  system  over  free 
outfall  weirs. 

The  "free  outfall"  is  a pipe  that  discharges  to  the  receiving 
waters.  It  is  truly  a free  outfall  if  the  receiving  waters  are  low  enough 
that  they  do  not  affect  the  hydraulic  head  at  the  outfall.  For  this  case 
the  water  surface  elevation  at  the  outfall  is  taken  as  the  minimum  value 
of  the  critical  and  normal  depths  at  the  given  discharge.  Discharge 
through  the  outfall  pipe  is  then  computed  in  the  same  manner  as  for  any 
other  link.  When  the  receiving  water  backs  up  into  the  pipe,  the  hydraulic 
head  at  the  outfall  is  taken  as  the  larger  of  the  free  discharge  heads 
and  the  receivina  water  elevation  of  the  outfall. 


When  there  is  a tide  gate  on  an  outfall  weir  or  a free  outfall, 
a check  is  made  to  see  whether  or  not  the  hydraulic  head  inside  the  gate 
exceeds  that  outside  of  the  gate.  If  it  does  not,  the  discharge  through 
the  outfall  is  equated  to  zero.  If  the  driving  head  is  positive,  then 
the  ARMCO  tide  gate  correction  is  applied  and  the  driving  head  is 
reduced  by  an  amount 

Ah  = | V2  exp  (-1.15  V//h  ) (1-8) 

where  h is  the  hydraulic  head  on  the  weir  or  the  free  discharge  head  in 
the  absence  of  the  tide  gate,  V is  the  velocity  over  the  weir  or  in  the 
conduit.  In  the  case  of  the  weir,  the  driving  head  in  Equation  1-4  or  1-6 
is  reduced  by  an  amount  Ah  and  the  discharge  recomputed.  For  the  free 
outfall,  the  hydraulic  head  at  the  outfall  is  raised  an  amount  Ah  thus 
decreasing  the  driving  head  in  the  conduit  by  that  amount. 

The  time  history  of  the  hydraulic  information  is  stored  on 
magnetic  tape  and  interfaced  with  the  Quality  Transport  Block.  It 
would  have  been  logical  to  combine  the  quality  and  transport  models, 
but  it  is  not  possible  because  of  computer  storage  requirements. 

TRANSPORT  QUALITY  BLOCK 


The  Transport  Quality  Block  is  a simple,  advective  flow  model. 

It  uses  the  quality  constituents  generated  by  the  Runoff  Block  and  the 
flows  computed  by  the  Transport  Block  in  combination  to  trace  the  quality 
constituent  through  the  sewer  system.  It  would  have  been  logical  to 
comoine  quality  and  quantity  portions  of  the  Transport  Block  but  this  was 
not  possible  because  of  computer  storage  requirements.  It  may  be  considered 
a post-processor  to  the  transport  model.  It  is  the  shortest  and  simplest 
of  the  model  routines. 

Output  from  both  the  Runoff  Block  and  the  Transport  Block 
serves  as  input  to  the  Quality  Transport  Block.  Data  from  both  models 
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is  transferred  via  magnetic  tape.  There  are  no  subroutines  used  by  this 
program  except  standard  library  routines.  All  operations  are  performed 
in  sequence,  i.e.  there  is  no  major  alternate  path  through  the  program. 

Mathematical  Idealization  and  Solution  Technique 


The  mathematical  basis  of  the  Quality  Transport  model  is 
the  equation  of  mass  conservation, 


where  M = mass  of  constituent  at  a section 

C = concentration  of  constituent  at  a section 
Q = flow  past  the  section 
t = time. 

Converting  mass  into  concentration  yields 

3C_  C TV  ,, 

at  ~ V at  ax 

where  U = velocity  of  flow  and  V = nodal  volume. 


(1-9) 


(1-10) 


The  rinht  hand  side  of  Equation  I - 1 0 represents  the  two  effects  which 
contribute  to  concentration  changes  at  a point.  The  first  term,  j , 
is  the  dilution  taking  place  at  the  section.  The  second  term, 

p 

U , is  the  advective  movement  of  the  concentration  gradient  past  the 
section . 

Figure  1-5  illustrates  the  advective  term  for  a single  conduit, 

n,  during  a small,  finite  time  At.  The  entire  gradient  slides  forward 

due  to  advective  transport  so  that  the  concentration,  C,  which  began  a 

distance  U At  from  node  2,  moves  to  node  2. 
n 
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In  the  model  we  assume  that  only  advective  flow  occurs  in  the 
pipes  and  that  all  dilution  occurs  in  the  nodes.  Thus,  in  model  terminology 
Equation  I - 1 0 becomes: 


AC 

AC 

For 

each 

pipe: 

n 

t 

. .O  • 

n 

C.  AC 

(1-11 

For 

each 

node: 

1 = 

:t 

- n1  AV  + a E ~ 
V-  n . t 

J n=up 

(1-12 

where  the  subscript  "n"  refers  to  the  pipe  and  the  subscript  "j"  refers 
to  the  node.  The  coefficiert  a is  a weighting  factor  based  on  flow, 

0,  which  accounts  for  the  averaging  of  flow  from  various  lines  at  the 
node. 


The  designation  "n  = up"  on  the  summation  in  Equation  1-12  is 
very  important.  It  means  that  only  those  pipes  upstream  from  node  j 
advect  material  to  the  node.  In  the  reverse  flow  condition,  it  refers 
to  pines  flowina  into  the  node.  It  is  assumed  that  node  j is  blind  to 
events  downstream  and  does  not  receive  water  qualitv  inputs  from  downstream. 
This  is  a valid  assumption  for  moderate  to  rapid  conduit  flow,  well 
substantiated  in  the  literature. 

In  the  model.  Equations  1-11  and  1-12  are  solved  alternately 
in  the  conduits  and  at  the  nodes,  so  that  a complete  time  history  of 
concentration  is  obtained  at  each  node. 


COMPUTER  PROGRAM  DESCRI PTI ON- -FORMAT 

Detailed  descriptions  of  each  of  the  computer  programs  which 
comprise  the  four  program  blocks  are  presented  below.  Presented  for  each 
block  is  a short  narrative  description  of  each  of  its  component  subroutines. 
Functional  flowcharts  of  program  logic  plus  a list  of  definition  of 
variables  in  FORTRAN  COMMON  storage  are  included  in  Chapter  V.  These 
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FIGURE  1-6  GENERAL  STRUCTURE  DISPLAY  BLOCK 
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descriptions,  together  with  the  program  listings,  should  enable  anyone 
familiar  with  computer  programming  to  identify  problem  areas  or  places 
for  desired  program  modifications. 

Flowcharts  depicting  the  general  structure  of  the  Display, 
Runoff  and  Transport  Blocks  and  Transport  Quality  Blocks,  are  indicated 
in  Figures  1-6,  1-7,  1-8,  and  1-9,  respectively.  Flowcharts  for  the 
component  routines  of  the  Display  Block  which  are  unchanged  from  the 
original  EPA  development,  can  be  found  in  Reference  1. 


3 LOCK 


The  logic  of  the  Display  Block  is  indicated  in  Figure  1-6. 

Subroutine  CURVE 

The  subroutine  CURVE  performs  the  following  operations: 

1.  Determines  maximum  and  minimum  of  arrays  to  be 
plotted. 

2.  Calculates  the  range  of  values  and  selects  appropriate 
scale  intervals.  In  the  Transport  Block,  a subroutine 
SCALE  is  included  to  perform  this  function. 

3.  Computes  vertical  axis  labels  based  upon  the 
calculated  scales. 

4.  Computes  horizontal  axis  labels  based  upon  the 
calculated  scales. 

5.  Joins  individual  parts  of  the  curve  by  subroutine  PINE. 

6.  Outputs  final  plot. 
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Subroutine  GRAPH 


The  graphing  subroutines  enable  hydrographs  and  pol 1 utographs 
to  be  plotted  on  the  printer  for  selected  locations  on  the  data  file. 

GRAPH  is  the  driving  subroutine,  and  it  calls  CURVE  to  produce  the  actual 
page  of  plotted  output.  In  the  Runoff  Block  subroutine  HCURVE  calls 
CURVE,  and  in  the  Transport  Block  subroutine  OUTPUT  calls  CURVE. 

The  subroutine  GRAPH  (IC)  operates  on  two  modes  which  are 
dependent  upon  the  value  of  IC  in  the  calling  sequence. 

If  IC  = 0 (when  called  by  the  Runoff  Block),  control  information 
is  read  from  cards. 

If  IC  = 1,  both  control  information  and  title  information  are 
read  from  cards. 

Subsequently,  both  options  join  and  the  subroutine  proceeds 
as  one  flow  sequence  as  follows: 

1.  Information  is  read  from  the  data  file  indicating 
the  structure  of  that  file. 

2.  An  array  ITAB  is  set  up  indicating  which  locations 
of  the  data  file  record  are  to  be  plotted. 

3.  All  hydrograph  and  pollutograph  information  is 
read  from  the  data  file. 

4.  For  each  type  of  hydrograph  and  pollutograph, 
individual  curves  are  selected,  transferred  into 
plotting  arrays,  and  outputted  in  a final  plotted 
form  by  subroutine  CURVE. 


FIGURE  1-7  GENERAL  STRUCTURE  RUNOFF  BLOCK 


Subroutine  PINE 

This  subroutine  joins  two  coordinate  locations  with  appropriate 
characters  in  the  output  image  array  A of  PPLOT. 

Subroutine  PPLOT 

This  subroutine  initializes  the  plotting  array,  stores 
individual  locations,  and  outputs  the  final  image  array  A for  the 
printer  plot. 

COMPUTER  PROGRAM  DESCRIPT1  N — RUN  FF  BLOCK 

The  overall  logic  for  the  Runoff  Block  is  shown  in  Figure  1-7. 

.'.‘jtin  Progran  RUNOFF 

This  is  the  main  routine  of  the  Runoff  Block.  This  program 
prints  a heading  and  then  acts  as  the  driver  routine  for  the  block. 

Figure  V-6  is  the  appropriate  flowchart. 

Sub  rout  "ne  GQUAL 

This  subroutine.  Figure  V- 7,  computes  a mass  balance  for  each 
gutter  at  each  time  step.  The  routine  accounts  for  the  mass  rate  of  flow 
into  a gutter  from  upstream  sources  as  well  as  that  from  tributary 
watersheds  and  local  catch  basins.  The  output  from  this  routine  is  an 
estimate  of  the  water  quality  expected  in  the  system  gutters  at  each 
time  step. 
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GUTTER 


The  function  of  subroutine  GUTTER  is  very  similar  to  that  of 
WSHED  and  is  shown  in  Figure  V-8.  It  calculates  a complete  set  of  water 
depth  and  flow  for  gutters  and  pipes. 

The  computation  also  proceeds  one  gutter  at  a time.  For  a 
gutter,  the  inflow  from  tributary  subcatchments  and  gutters  is  first 
computed.  Newton's  iterative  procedure  is  again  used  to  determine  the 
depth  and  outflow  of  gutters  so  that  the  mass  (volume)  of  water  is 
conserved.  The  flow  is  computed  by  Manning's  equation.  The  hydraulic 
radius  of  trapezoidal  gutters  and  circular  pipes  is  calculated  separately 
in  different  paths  of  the  program. 

A pipe  may  surcharge  when  it  is  full  and  the  inflow  is  larger 
than  the  outflow  capacity.  In  this  case,  the  surcharged  amount  will 
be  computed  and  stored  at  the  head  end  of  the  pipe.  A message  will  be 
printed  to  indicate  the  time,  location,  and  total  amount  of  the  surcharge. 
The  oipe  will  remain  full  until  the  stored  water  is  completely  drained. 

Subroutine  HCURVE 

HCURVE,  Figure  V- 9,  performs  some  functions  of  subroutine 
GRAPH,  and  is  designed  to  call  subroutine  CURVE  to  plot  the  hyetographs 
and  the  inlet  hydrograph. 

■1  r <~.'b  HYDRO 

This  subroutine,  shown  in  Figure  V-10,  computes  the  hydrograph 
coordinates  with  the  assistance  of  three  core  subroutines,  i.e.  RHYDRO, 
WSHED,  and  GUTTER,  as  shown  in  Figure  1-7.  It  initializes  all  the 
variables  to  zero  before  calling  RHYDRO  to  read  in  the  rainfall  hyetograph 
and  information  concerning  the  inlet  drainage  basin.  According  to  the 
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upstream  and  downstream  relationship,  the  subroutine  sequences  the 
computational  order  for  gutters/pipes. 

A DO  loop  is  formed  to  compute  the  hydrograph  coordinate  for 
each  incremental  time-step.  In  each  step,  subroutine  WSHED  is  first 
called  to  calculate  the  rate  of  water  flowing  out  of  the  idealized 
subcatchments.  GUTTER  is  then  called  to  route  the  flow,  according  to 
the  input  from  tributary  subcatchments  and  gutters.  Water  flowing  into 
the  inlet  point,  be  it  from  gutters  or  direct  drainage  from  subcatchments, 
is  added  up  for  a hydrograph  coordinate. 

During  the  process  of  computation,  an  accounting  is  made  for 
the  deposition  of  rainfall  water  in  the  form  of  runoff,  detention,  and 
infiltration  loss.  A mass  continuity  can  therefore  be  checked  and 
printed  for  reference. 

Finally,  the  rainfall  hyetograph  and  the  inlet  hydrograph  are 
plotted  as  an  output.  The  control  is  then  returned  to  subroutine  RUNOFF. 

Subroutine  QSHED1 

Subroutine  QSHED1 , shown  in  Figure  V-ll,  is  used  to  estimate 
the  initial  mass  of  pollutant  on  the  unit  watersheds  and  in  the  system 
catch  basins  at  the  beginning  of  a storm.  This  is  done  by  applying 
empirically  determined  pollutant  buildup  factors  for  the  number  of  dry 
days  prior  to  a storm  and  then  reducing  the  total  by  the  amount  taken 
up  by  street  cleaning. 

Subr  ’ utine  QSHED2 

QSHED2 , Figure  V-12,  is  used  to  estimate  the  mass  rate  of 
pollutant  runoff  from  watersheds  and  catch  basins.  The  routine  first 
loops  through  each  unit  watershed  and  makes  an  estimate  of  pollutnat 
runoff  based  on  the  amount  of  pollutant  on  the  watershed  and  the  rate 
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of  flow  from  the  watershed.  A similar  calculation  is  made  through  an 
empirical  relationship  for  each  system  catch  basin.  The  pollutant 
runoff  estimates  are  then  supplied  to  subroutine  GQUAL  for  estimates  of 
gutter  quality.  It  should  be  recognized  that  because  of  the  assumed 
proportionality  between  quality  constituents  and  type  of  land  use  it 
is  necessary  to  route  only  as  many  individual  items  as  there  are  land 
use  categories.  This  combined  routine  is  later  disaggregated  into  the 
numerous  quality  constituents  of  interest. 

S routine  RECAP 

This  subroutine.  Figure  V-13,  reads  the  tapes  created  by  the 
Runoff  Block.  It  then  writes  a summary  report,  either  of  hydraulic 
results  alone  of  hydraulic  and  quality  results,  dependent  on  the 
specification  of  the  run. 

Subroutine  RHYDRO 

This  subroutine  is  called  by  HYDRO  to  read  input  data  related 
to  the  subcatchment  areas  and  to  perform  some  initial  preparatory  work, 
such  as  unit  conversion  and  error  detection.  A normal  execution  of 
RHYDRO  should  provide  all  the  necessary  information  for  the  calculation 
of  a runoff  hydrograph.  Figure  V-14  shows  the  flowchart  for  subroutine 
RHYDRO. 

There  are  four  basic  categories  of  input  data.  The  general 
information  includes  a number  representing  the  subcatchment  area,  period 
of  simulation,  and  a key  indicating  if  the  rainfall  hyetograph  is 
spatially  different  from  that  of  the  previous  basin.  A new  rainfall 
hyetograph  will  be  read  if  it  is  so  indicated.  Otherwise,  that  part 
of  the  read  operation  will  be  skipped  and  the  rainfall  of  the  previous 
inlet  drainage  basin  will  be  used.  The  first  basin  must  have  a rainfall 
input. 
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The  program  proceeds  to  read  subcatchment  data,  e.g.,  the 
size,  width,  ground  slope.  The  gutter  information  is  read  soon  afterward. 

It  must  be  noted  that  the  program  can  detect  only  logical  errors 
such  as  indexing  numbers.  However,  the  input  data  are  tabulated  by  the 
computer  to  allow  a check  against  the  original  for  absolute  correctness. 

Subroutine  WSHED 

This  subroutine  computes  the  depth  and  flow  rate  of  water 
overland.  The  logic  of  subroutine  WSHED  can  be  seen  in  Figure  V - 1 5 . 

As  shown  in  Figure  1-7,  the  subroutine  is  called  by  HYDRO  at  each 
incremental  period  of  integration.  During  that  period,  the  rainfall 
intensity  is  first  interpolated  from  the  designated  rainfall  hyetograph 
for  each  subcatchment.  This  rainfall  intensity  is  assumed  uniform  over 
each  subcatchment. 

A DO  loop  is  set  up  to  treat  the  subcatchments,  one  at  a time. 

For  a subcatchment,  the  amount  of  infiltration  loss  is  calculated  using 
Horton's  equation, 

Infiltration  loss  = fQ  + (f . - fQ)e”xt  (1-13) 


where  f0,  f ^ and  a are  coefficients  and  t is  the  time  from  the  start 
of  rainfall.  The  loss  is  compared  with  the  amount  of  water  existing 
on  the  subcatchment  plus  the  rainfall.  If  the  loss  is  larger,  it  is 
set  equal  to  the  amount  available  and  the  remainder  of  the  computation 
is  skipped. 

The  water  depth  will  thus  increase  without  inducing  an  outflow 
until  it  reaches  the  specified  detention  requirement.  Beyond  that,  the 
outflow  rate  is  calculated  by  Manning's  equation  using  depth  as  the 
hydraulic  radius.  An  iterative  procedure  using  Newton's  technique  is 
established  to  determine  the  water  depth  and  the  outflow  rate  so  that 
the  continuity  of  water  mass  is  satisfied. 
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Upon  completion,  the  subroutine  will  return  with  a set  of 
water  depths  on  each  subcatchment  for  the  next  time-step.  It  also 
produces  the  flow  necessary  for  subsequent  routing  in  the  gutters. 

COMPUTER  PROGRAM  DESCRIPTION— TRANSPORT  BLOCK 

The  overall  logic  of  the  TRANSPORT  BLOCK  is  shown  in  Figure  1-8. 

Main  Program  MAIN 

* This  is  the  driving  or  executive  program  of  the  Transport 
Block.  In  addition,  it  performs  a modified  Euler  solution  of  the 
motion  and  continuity  equations.  The  necessary  data  for  print  requirements 
is  saved  by  this  program.  Figure  V-16  is  a diagram  of  its  logic. 

Subroutine  BOUND 

This  subroutine  computes  the  outflow  "QOU"  for  each  node. 

Both  internal  and  external  transfers  are  computed  as  indicated  in 
Figure  V-17. 

Subroutine  DEPTH 

This  subroutine  finds  the  critical  depth  and  the  normal  depth 
in  a conduit,  corresponding  to  a computed  flow  rate.  The  procedure  is 
diagrammed  in  Figure  V-18. 

Subroutine  HEAD 

This  subroutine  converts  nodal  depths  to  pipe  depths  as 
indicated  in  Figure  V-19.  A storage  surface  area  is  also  assigned 
to  the  proper  nodes. 
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FIGURE  1-8  GENERAL  STRUCTURE  TRANSPORT  BLOCK 
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c' uv routine  HYDRAD 

This  subroutine  shown  in  Figure  V-20  uses  tables  of  normalized 
values  to  compute  the  hydraulic  radius,  the  surface  width,  and  the  "net" 
cross-sectional  area  of  each  conduit. 

Subroutine  IN DATA 

INDATA,  Figure  V-21,  reads  and  prints  all  input  data  except 
for  hydrograph  cards,  which  are  read  by  subroutine  INFLOW  as  required. 

It  also  performs  some  initialization.  All  node-conduit  linkages  are 
set  up  and  converted  to  the  interval  numbering  system.  The  first  two 
hydrographs  from  both  tape  input  and  card  input  are  read  and  converted 
to  internal  numbers. 

SiAbroutine  INFLOW 

This  subroutine  selects  the  input  hydrograph  ordinate  from 
tape  or  card  input.  Figure  V- 22  is  the  appropriate  flowchart. 

Subroutine  OUTPUT 

As  shown  in  Figure  V- 23,  this  subroutine  prints  the  nodal 
information  as  a time  history  of  the  hydraulic  qrade  line.  It  also 
prints  the  time  history  of  conduit  flows  and  velocities.  The  plot 
routines  are  controlled  by  this  subroutine. 

S jjr-ou  * ine  TIDCF 

Subroutine  TIDCF  uses  a least-square  procedure  to  calculate 
the  coefficients  of  the  tidal  function  H(T)  = A1  + A2  SIN(T)  + A3  SIN(2T) 
+ A4  S IN ( 3T ) + A5  CUS(T)  + A6  C0S(2T)  + A7  C0S(3T)  from  input  values  of 
H and  T.  The  logic  is  shown  *n  Figure  V-24. 
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MAIN,  Figure  V-25,  is  the  executive  routine  for  the  Quality 
Transport  Block.  After  initial ization,  the  program  does  a loop  with 
time  as  the  increment.  It  calls  INPUT  and  computer  nodal  volumes  and 
the  velocity  in  each  conduit.  The  pump  quality  changes  are  then  computed, 
and  then  the  advective  transport  through  conduits.  Internal  linkage  mass 
transfers  are  computed,  and  then  all  node  concentrations  are  updated. 
Necessary  values  are  saved  for  subsequent  printing  or  plotting.  After 
the  completion  of  the  requested  time  steps,  subroutine  OUTPUT  is  called 
and  program  is  completed. 

Subroutine  HYDRAD 

This  subroutine  performs  the  same  function  in  this  block  as 
it  does  in  the  Transport  Block.  The  reader  is  referred  to  the  Transport 
Block.  The  logic  is  shown  in  Figure  V-20. 

Subroutine  INDATA 

This  subroutine,  Figure  V- 26,  reads  the  job  description  and 
control  parameters.  It  also  reads  the  node  numbers  for  detailed  printing. 
Initial  information  from  the  input  tape  units  is  read  in.  Initial 
pol 1 utographs  and  hydrographs  are  read  in,  and  checks  are  made  fo 
consistency  in  time  step  requirements. 
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FIGURE  1-9  GENERAL  STRUCTURE  QUALITY  TRANSPORT  BLOCK 
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Subroutine  INPUT 


This  subroutine  does  the  necessary  interpolation  between  data 
points  for  both  quality  data  and  hydraulic  data.  If  new  data  is  required, 
it  is  read  from  the  appropriate  input  tape  unit  prior  to  interpolation. 

The  logic  ii  indicated  in  Figure  V- 27. 

Subroutine  OUTPUT 

This  subroutine  processes  the  stored  computed  nodal  qualities 
and  prints  them  as  shown  in  Figure  V- 28. 
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CHAPTER  II 

MODEL  SPECIFICATIONS 


INTRODUCTION 

This  chapter  outlines  the  specifics  of  the  computer  programs 
as  they  relate  to  data  preparation  and  operation  on  various  computer 
systems.  Estimates  are  given  for  computer  core  requirements  and 
peripheral  equipment,  and  estimates  of  expected  execution  times. 

In  the  preparation  of  these  computer  programs  every  effort  has 
been  made  to  make  them  as  nonspecific  as  possible  in  regard  to  particular 
computer  systems.  However,  if  problems  of  program  installation  develop 
which  cannot  be  resolved  with  the  use  of  this  document,  the  program's 
user  is  urged  to  contact  Water  Resources  Engineers,  Inc.,  Walnut  Creek, 
California,  for  aid  in  resolving  the  difficulty. 


PROGRAM  LANGUAGE 

All  blocks  of  the  model  were  originally  programmed  for  the 
UNIVAC  1108  in  FORTRAN  V.  This  version  of  the  FORTRAN  compiler  is 
essentially  compatible  with  IBM  FORTRAN  IV  level  G and  versions  of  the 
model  have  been  executed  on  IBM  equipment.  This  FORTRAN  is  also 
compatible  with  the  Extended  FORTRAN  compiler  used  on  CDC  6000  series 
equi pment. 
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OPERATING  REQUIREMENTS 


The  model  consists  of  three  discrete  blocks,  each  of  which 
has  separate  operating  requirements.  The  requirements  for  the  Display 
Block,  which  is  included  in  the  other  blocks,  are  included  in  the  figures 
given  below. 

RUNOFF  BLOCK 

Highspeed  core  required  ~ 126,000g  words 

~ 44, 000] g words 

~ T 80 , 000-j g bytes 

4 drum,  disk  or  tape  files:  1 output,  3 scratch 

1 card  reader 

1 line  printer 

TRANSPORT  BLC.  'K 

Highspeed  core  required  ~ 135,000s  words 

~ 47,500] g words 

~ 190,000] g bytes 

2 drum,  disk  or  tape  files:  1 output 

1 input  (generated  by  Runoff  Block) 

1 card  reader 

1 line  printer 

TRANSPORT  QUALITY  BLOCK 

Highspeed  core  required  ~ 126,000g  words 

= 44,000]q  words 

~ 180,000]0  bytes 
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2 drum,  disk  or  tape  files:  1 input  generated  by  Runoff  Block 

1 input  generated  by  Transport  Block 

1 card  reader 
1 line  printer 


TYPICAL  EXECUTION  TIMES 


Typical  program  execution  times  for  the  various  program  blocks, 
including  the  Display  Block,  are  given  as  follows. 

RUNOFF  BLOCK 


Execution  time  is  approximately  proportional  to  the  number  of 
watersheds,  gutters  and  the  number  of  time  steps  desired.  In  an  example 
problem  involving  50  watersheds,  58  gutters,  and  135  time  steps,  76 
seconds  were  required  on  the  UNI  VAC  1108. 


TRANSPORT  BLOCK 


Execution  time  is  approximately  proportional*  to  the  number 
of  conduits  and  the  number  of  time  steps.  As  an  example,  a system 
involving  50  conduits  and  40  time  steps  required  311  seconds  on  UNIVAC  1108. 

TRANSPORT  QUALITY  BLOCK 


Execution  time  is  approximately  proportional*  to  the  number 
of  conduits  and  the  number  of  time  steps.  As  an  example,  a system 
involving  7 conduits  and  450  time  steps  required  30  seconds  on  the 
UNIVAC  1108. 


‘Computation  is  affected  to  some  extent  by  the  volume  of  printing,  hence 
the  relatively  long  execution  time  for  transport  quality. 
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JOB  CONTROL 


The  assignment  of  logical  units  for  system  I/O  requires,  in 
general,  the  provision  for  files  to  be  written  on  specific  physical 
devices.  To  accomplish  this  the  programmer  must  supply  the  necessary 
job  control  cards.  As  a rule,  job  control  is  highly  machine  and 
installation  dependent;  in  fact,  it  often  differs  on  two  identical 
machines  at  different  installations.  Therefore,  the  Runoff  and  Transport 
Models  cannot  include  a job  control  definition  that  is  universally  applicable. 


At  most  installations  the  logical  unit  corresponding  to  the  card 
reader  is  given  the  number  5 and  the  line  printer  is  given  the  number  6. 
The  Runoff  and  Transport  Models  are  programmed  on  the  assumption  that 
units  5 and  6 are  so  used.  Typically,  the  systems  programmers  have 
provided  the  necessary  specification  for  these  units  and  also  for  the 
card  punch. 


Data  file  and  scratch  file  assignments  require  specification 
to  be  supplied  for  each  unit.  The  rules  for  such  assignments  must  be 
ascertained  from  the  systems  programmers  at  the  installation,  since 
there  is  considerable  variation  in  unit  number  availability,  etc.  In 
general  , one  should  only  set  up  the  units  needed  in  a given  run,  since 
there  may  be  a charge  for  file  space  that  is  reserved,  even  if  it  is 
not  used. 


INPUT/OUTPUT  FILES 

The  total  simulation  package  is  designed  to  be  run  as  separate 
blocks,  with  the  Display  Block  used  as  needed  for  plotting.  Interfacing 
is  accomplished  between  blocks  by  magnetic  tape,  drum  or  disk  storage. 
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The  Runoff  Block  can  create  four  tapes,  two  of  which,  JOUT  (1) 
and  JOUT  (2),  are  the  hydraulics  tape  and  quality  tape  respectively  for 
interfacing  with  the  Transport  Block  and  Transport  Quality  Block.  Write 
statements  are  unformatted*  for  these  two  tapes.  The  print  tape  and  plot 
tape  are  JOUT  (3)  and  JOUT  (4),  and  write  statements  are  formatted. 

The  Transport  Block  designates  the  input  hydraulics  tape  from 
the  Runoff  Block  as  N21.  The  output  hydraulics  tape  designation  is  N22. 
The  write  statements  for  N22  are  formatted. 

The  Quality  Transport  Block  designates  the  input  hydraulics 
tape  from  the  Transport  Block  as  N22  and  the  quality  tape  from  the  Runoff 
Block  as  N21 . 


CARD  FORMATS,  INPUT  DATA 

The  following  tables  present  the  specific  card  formats  required 
for  use  of  the  Runoff  Block,  the  Transport  Block  and  the  Transport  Quality 
Block.  Table  1 1 - 1 gives  the  specifications  for  the  Runoff  Block,  Table  1 1 -2 
covers  the  Transport  Block  and  Table  1 1 -3  defines  the  Transport  Quality 
Block . 


The  following  figures  present  deck  setup  schematics  for  the 
models:  Figure  1 1 - 1 is  a schematic  of  a complete  urban  stormwater  model 

deck  setup;  Figure  1 1 -2  is  a schematic  of  a deck  setup  for  the  Runoff 
Block;  Figure  1 1 - 3 is  a deck  setup  schematic  for  the  Transport  Block; 
and  Figure  1 1 -4  is  a deck  setup  schematic  for  the  Transport  Quality  Block. 

All  values  should  be  right  justified  in  the  card  columns  indicated, 
and  the  cards  should  be  arranged  for  input  in  the  same  order  as  described. 
Unless  otherwise  noted  all  cards  defined  for  input  must  appear  in  the  job 
input  stream. 


♦Formating  is  the  procedure  for  placing  information  on  the  tapes. 
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TABU  I l-l 

RUNOFF  BLOCK  CARD  DATA 


Card 

Group 

For-Lit 

Cars 
Co) j-ns 

Description 

Van  able 
Na.  e 

Default 

Value 

1 

6110 

1-10 

Name  of  input  tape  (none  at  present). 

JIN(l) 

None 

11-20 

Name  of  input  tape  (none  at  present). 

J IN  ( 2 ) 

None 

21-33 

Na~e  of  quantity  output  tape. 

J0'JT(1) 

None 

31-43 

f.a-e  of  quality  output  tape. 

J0UT ( 2 ) 

None 

41-50 

Name  of  print  file  tape. 

J0UT( 3) 

None 

51-60 

Name  of  plot  file  tape. 

J0UT ( 4 ) 

None 

2 

20A4 

1-80 

Title  cards:  two  cards  with  heading 

to  be  printed  on  output. 

TITLE 

None 

3 

Control  card:  one  card. 

215 

1-  5 

Number  of  Basin. 

INLET 

None 

6-1 C 

Number  of  time-steps  to  be  calculated. 

NSTEP 

None 

13 

11-13 

Hour  of  start  of  storm  (24-hour  clock). 

NHR 

None 

12 

14-15 

Minutes  of  start  of  stonn. 

NMN 

None 

F5.1 

16-20 

Integration  period  (mi n). 

DELT 

None 

15 

21-25 

Number  of  hyetographs. 

NRGAG 

None 

F5.Q 

26-30 

Percent  of  impervious  area  with  zero 
detention  (immediate  runoff). 

PCTZER 

25.0 

4 

Rainfall  control  card. 

15 

1-  5 

Number  of  data  points  for  each 
hyetograph. 

NHIST0 

None 

F5.0 

6-10 

Time  interval  between  values  (min). 
(Must  be  even  multiple  of  DELT  above) 

TH1ST0 

None 

REPEAT  CARD  GROUP  5 FOR  EACH  HYETOGRAPH. 

5 

Rainfall  hyetograph  cards:  10  intervals 

per  card. 

10F5.0 

1-  5 

Rainfall  intensity,  first  interval 
(in. /hr. ). 

RAIN ( 1 ) 

None 

6-1 0 

Rainfall  intensity,  second  interval 
(in. /hr.). 

RAIN(2) 

None 

11-15 

Rainfall  intensity,  third  interval 
(in. /hr. ). 

RA I N ( 3) 

None 

16-20 

Rainfall  intensity,  fourth  interval 

RAIN ( 4 ) 

None 

(in. /hr. ) 
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TABLE  II-l 
(Continued) 


Card 

Card 

Variable 

Def  Jul  t 

Group 

forr.it 

Colu  ns 

Description 

Name 

Value 

REPEAT  C*  6 F G J T T £ R/  P J P £ . 

6 

Gutter/pipe  cards:  one  card  per  gutter/ 

pipe  (if  rare,  leave  out). 

216 

1-  6 

Gutter  nu  rer. 

NAMEG 

None 

7-12 

Gutter  or  Inlet  number  for  drainage.* 

NGTO 

None 

13 

13-15 

/ = 1 for  gutter 
^ = 2 for  pipe 

NP 

None 

5X 

16-20 

Blank 

7F8.0 

21-28 

Bottom  width  of  gutter  or  pipe 
diameter  (ft). 

GWIDTH=G1 

None 

29-36 

Length  of  gutter  (ft). 

GLEN  =G2 

None 

37-44 

Invert  slope  (ft/ft). 

GSL0PE=G3 

None 

45-52 

Left-hand  slope  (ft/ft),  horizontal/ 

GS1  =G4 

None 

vertical 

53-60 

Right-hand  slope  (ft/ft),  horizontal/ 

GS2  =G5 

None 

verti  cal 

61-68 

Manning's  coefficient. 

GN  =G6 

None 

69-76 

Depth  of  gutter  when  full  (in). 

DFULL  =G7 

10.  in. 

NOTE: 

Blank  card  to  terminate  gutter  cards: 
one  card. 

REPEAT  CARD  7 FOR  EACH  SUBCATCHMENT. 

7 

Subcatch-ent  cards  (216,  13,  10F5.0, 
FI  0.5):  one  card  per  subcatchment. 

216 

1-  6 

Subcatchment  number. 

NAMEW 

None 

7-12 

Gutter  or  inlet  number  for  drainage. 

NGTO 

None 

1-3 

13-15 

Nyetograph  number. 

NHYET= JK 

None 

10F5.0 

16-20 

Width  of  subcatchment  (ft). 

WWI0TH=W1 

None 

21-25 

Area  of  subcatchment  (acres). 

WAREA  =W2 

None 

26-30 

Percent  imperviousness  of  subcatchment. 

PCIMP  =W3 

None 

31-35 

Ground  slope  (ft/ft). 

WSL0PE=W4 

0.030 

36-40 

Impervious  area  > 

W5  =W5 

0.013 

\ Resistance  Factor 

41-45 

Pervious  area  J 

W6  =W6 

0.250 

46-50 

Impervious  area  h Retention  storage 

WST0RE=W7 

0.062 

51-55 

„ • /(in.) 

Pervious  area  ■> 

W$T0RE=W8 

0.184 

‘Gutter  or  inlet  to  which  gutter  drains. 
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TAB  It  11-1 
(Continued) 


Card 

Card 

Variable 

Defa  jl  t 

Group 

Format 

Columns 

Description 

Name 

Value 

7 

10F5.0 

(cont'd) 

56-60 

Maxir^m  infiltration  rate  {in. /hr.). 

ULMAX=„9 

3. DO 

61-65 

Minimum  infiltration  rate  (in. /hr.). 

ULMlN=W10 

0.52 

F10. 5 

66-75 

Decay  rate  of  infiltration  (1/sec). 

DECAY=W1 1 

0.00115 

F5.0 

75-80 

Maximum  alio. .'able  infiltration  (inches). 

DEPIN-U12 

12. inches 

NOTE: 

Blank  card  to  terminate  subcatchment  cards 
one  card. 

8 

110 

1-10 

Nurber  of  quality  constituents. 

NQS 

None 

2E10.0 

11-20 

Number  of  dry  days  since  last 
precipi  tation 

DRY DAY 

None 

21-30 

Street  cleaning  frequency  (days). 

CLFREQ 

None 

no 

31  -40 

Number  of  passes  per  cleaning. 

NPASS 

None 

E1J.0 

41-50 

Catchbasin  volume  (cubic  feet). 

CBVOL 

None 

OPTIONAL  CARDS,  SKIP  IF  NQ's  ON  CARD  GROUP 

8 ABOVE  IS 

ZERO 

9 

REPEAT  CARD  GROUP  9 FOR  EACH  SUBCATCHMENT. 

2110 

1-10 

Subcatchment  number. 

N 

None 

11-20 

Land  use  type. 

KLAND=KL 

None 

=1  for  single  family  residential 
=2  for  multi-family  res:  lent i a 1 
= 3 for  commercial 
=4  for  industrial 
=5  for  undeveloped  or  park  lands. 

2E10.0 

21-30 

Number  of  feet  of  gutter  in 

GQLEN= 

None 

subcatchments. 

GQ(100  FT) 

31-40 

Equivalent  number  of  standard 
catchbasins. 

BASINS=BA 

None 

10 

lnlet/gutter  print  control : one  card. 

215 

1-  5 

Number  of  inlets/gutters  for  which 
flows  are  to  be  printed. 

NPRNT 

None 

6-10 

Number  of  time-steps  between  printings. 

INTERV 

None 

IF  NPRNT'O,  SKIP  CARDS  11. 

Inlet/Gutter  print  cards:  16  values  per  card. 


1615 

1-5 

Inlet/Gutter  numbers  for  which 
flows  are  to  be  printed. 

IPRNT(l) 

None 

6-10 

(Maximum  of  200) 

I PRNT ( 2 ) 

None 

11-15 

I PRNT  ( 3 ) 

None 

IPP.NT(NPRNT) 

None 
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TABLE  II-l 
(Continued) 


Card 

Group 

For  at 

Card 

Col-  ns 

Cescri ption 

Variable 

Name 

Default 

Value 

12 

13 

1-5 

‘.-■Ler  of  gutter  and  inlets  to  be 
plotted. 

NPLOT 

None 

13 

If  r.PLOT’O,  Sr.IP  CAPOS  13  and  14 

Gdtter  a".;  inlet  plot  cards: 

16  values  per  carp. 

1615 

1-5 

Inlets  and  gjtters  wmch  are  to  be 
plotted. 

I PLOT ( 1 ) 

None 

6-10 

I PLOT  (2)1 

None 

11-15 

I PLOT ( 3 ) 

None 

IPLOT(NPLOT) 

None 

14 

Constituents  to  be  plotted,  0 means 
no  plot,  1 means  plot. 

2511 

I 

F low 

None 

None 

2 

Settleable  Solids 

None 

None 

3 

Suspended  Solids 

None 

None 

4 

Total  Dissolved  Solids 

None 

None 

5 

BOD 

None 

None 

6 

COD 

None 

None 

7 

Chlorides 

None 

None 

8 

S04 

None 

None 

9 

Grease 

None 

None 

10 

Total  Col i form 

None 

None 

11 

Fecal  Col i form 

None 

None 

12 

Ammonia 

None 

None 

13 

Organic  N 

None 

None 

14 

Nitrate  plus  nitrite 

None 

None 

15 

Total  hydrolized  phosphorous 

None 

None 

16 

Orthophosphate 

None 

None 

17 

Mercury 

None 

None 

18 

Copper 

None 

None 

19 

Zinc 

None 

None 

20 

Lead 

None 

None 

21 

Chromi  urn 

None 

None 

22 

Cadmi  urn 

None 

None 

23 

Arseni  c 

None 

None 

NOTE:  The  arrangement  of  constituents  to  be 
plotted  will  apply  to  all  gutters  and 
inlets.  No  changes  can  be  made  between 
plots. 
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TABLE  11-2 

TRANSPORT  BLOCK  PROGRAM  DATA 


m 

Card 

Variable 

Defaul  t 

Format 

Columns 

Description 

Name 

Value 

1 

1 5A4 

1-60 

Description  of  computer  run,  will  be 
printed  on  output.  (2  cards). 

ALPHA 

None 

2 

15,  2F5.0, 
8I5.4X.A6, 
15 

1-5 

Number  of  integration  steps  or  time 
cycles  desired.* 

NTCYC 

None 

6-10 

Length  of  integration  step,  seconds. 

DELT 

0.0 

11-13 

Start  time,  hours. 

NHR 

0 

14-15 

Start  time,  minutes 

NMN 

0 

16-20 

Number  of  nodes  for  detailed  printed 
information  for  each  cycle  (20  maximum) 

NHPRT 

None 

21-25 

Number  of  conduits  for  detailed  printed 
information  for  each  cycle  (20  maximum) 

NQPRT 

None 

26-30 

Number  of  nodes  to  be  plotted 
(18  maximum) 

NPLT 

None 

31-35 

Number  of  conduit  flows  to  be  plotted 
(13  maximum). 

LPLT 

None 

36-40 

Start  print  cycle 

NS  TART 

1 

41-45 

Interval  between  print  cycles  (maximum 
number  of  cycles  printed  is  100, 

NTCYC  - N START  ,ln) 

INTER 

1 

ITnTE  IUU) 

46-50 

Number  of  input  nodes,  if  card  input 
hydrographs  are  used. 

NJSW 

None 

51-55 

Hydrograph  input  tape 

N21 

0 

60-65 

Name  of  saved  hydrograph  input  file  (N21) 

NAME 

None 

66-70 

Quality  input  tape 

N22 

0 

3 

3110 

1-10 

First  node  number  for  detailed  printing 

JPRT  (1) 

None 

11-20 

Second  node  number,  up  to  number  of 
nodes  defined  by  NHPRT. 

JPRT  (2) 

None 

♦Suggested  for  initial  simulation  run. 


NTCYC 


(NSTEP  x PELT)  of  runoff  block  x 60  sec,  per  min. 
DELt  of  transport  block 
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TABLE  1 1-2 
( Continued ) 


Card 

Card 

Variable 

Default 

Group 

Format 

Columns 

Description 

Name 

Value 

4 

8110 

1-10 

First  conduit  number  for  detailed  printing. 

CPRT  (1) 

None 

11-20 

Second  conduit  number  up  to  number  of 
nodes  defined  by  NQPRT , 

CPRT  (2) 

None 

(FOR  ELEVATIONS  AND/OR  FLOWS) 

OPTIONAL 

5 

8110 

1-10 

First  node  number  for  plotting 

JPLT  (1) 

None 

11-20 

Second  node  number,  up  to  number  of  nodes 
defined  by  NPLT.  (This  option  is  for 
junction  on  water  surface  elevations) 

JPLT  (2) 

None 

6 

8110 

1-10 

First  conduit  number  for  plotting 

KPLT  (1) 

None 

11-20 

Second  conduit  number  for  plotting,  up 
to  number  of  nodes  defined  by  LPLT 
(This  option  is  for  conduit  flow  rate). 

KPLT  (2) 

None 

1 CARD/CONDUIT 

7 

415, 

1-5 

Conduit  number 

NC0ND  (N) 

None 

7F5.0 

6-10 

Node  number  at  one  end  of  conduit  upstream 

NJUNC (N ,1 ) 

None 

11-15 

Node  number  at  other  end  of  conduit 
downstream 

NJUNC ( N ,2 ) 

None 

16-20 

Type  of  conduit  shaDe 

NKLASS(N) 

None 

1 = circular 

2 = rectangular 

3 = horseshoe 

4 = egg 

5 = baskethandle 

6 = trapezoid  ( Open  Channel) 

21-25 

Cross  sectional  area  of  conduit,  Sq.  Ft. 
(necessary  only  for  types  3,4,  and  5). 

AFULL  (N) 

None 

26-30 

Vertical  depth  of  conduit,  feet 

DEEP  (N) 

None 

31-35 

Maximum  width  of  conduit,  feet 

WIDE  (N) 

0.0 

(for  Trapezoid,  Bottom  Width; 
for  Pipe,  Dia.) 

36-40 

Length  of  conduit,  feet 

LEN  (N) 

None 

41-45 

Distance  of  conduit  invert  above  node 

ZP  ( N, 1 ) 

0.0 

invert  at  NJUNC  (N,l)  (upstream) 
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TABLE  1 1-2 
(Continued) 


Card 

Card 

Variable 

Defaul  t 

Group 

Format 

Columns 

Description 

Name 

Value 

7 

415, 

46-50 

Distance  of  conduit  invert  above  node 

ZP  (N.2) 

None 

(Con't) 

7F5.C 

invert  at  NJUNC  (N,2)  (downstream) 

51-55 

Mannings  coefficient 

ROUGH  (N) 

0.014 

66-70 

Slope  on  one  trapezoid  side 
(horizontal/vertical ) 

STHETA  (N) 

0.0 

71-75 

Slope  on  other  trapezoid  side 
(horizontal/vertical ) 

SPHI  (N) 

0.0 

(Last  card  must  have  99999  in  columns  1 to  5 

1 CARD/ NODE 

8 

15, 

1-5 

Node  number  (none  greater  than  9999). 

JUN  (J) 

None 

3F5.0 

6-10 

Ground  elevation,  feet. 

GRELEV  (J) 

0.0 

11-15 

Invert  elevation,  feet. 

Z (J) 

0.0 

16-20 

Net  constant  flow  into  nodes, 
cfs  (optional) 

QINST  (J) 

0.0 

(Last  card  must  have  a 99999  in  columns  1 to  5) 

1 CARD/ORIFICE 

9 

215, 

2F5.0 

1-5 

Junction  containing  orifice. 

NJUNC  (N,1) 
NJUNC  ( N , 2 ) 

None 

6-10 

Junction  to  which  orifice  discharges 

None 

11-15 

Orifice  area  in  sq.  ft. 

A0RIF  (I) 

0.0 

16-20 

Orifice  discharge  coefficient. 

C0RIF  (I) 

0.0 

(Last  card  must  have  99999  in  columns  1 to  5) 

1 CARD/WEIR 

NJUNC  (N ,1 ) 

10 

315, 

1-5 

Node  at  which  weir  is  located. 

None 

5F5.0 

6-10 

Conduit  to  which  weir  discharges 

NJUNC  (N ,2) 

None 

(NOTE:  If  free  outfall  weir,  set 

NJUNC  (N,2)  equal  to  zero). 
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TABLE  II -2 
(Continued) 


A 


Variable 

Default 

Group 

Format 

Columns 

Description 

Name 

Value 

10 

315 

11-15 

Type  of  weir 

KUEIR  (I) 

None 

(Con't) 

5F5.0 

1.  transverse 

2.  transverse  with  tide  gates 

3.  side  flow 

4.  side  flow  with  tide  gates 

16-20 

Height  of  weir  above  invert,  feet. 

YCREST  (I) 

0.0 

21-25 

Height  to  top  of  weir  opening,  above 
invert,  feet. 

YT0P  (I) 

0.0 

26-30 

Weir  length,  feet. 

WLEN  (I) 

0.0 

31-35 

Coefficient  of  discharge  for  weir. 

C0EF  (I) 

0.0 

36-40 

Coefficient  of  discharge  for  surcharged 
condition. 

C0EFS  (I) 

0.0 

(Last  card  must  have  a 99999  in  columns  1 to  5) 

1 CARD/ PUMP 

11 

215, 

1-5 

Node  being  pumped 

NJUNC  (H,l ) 

None 

7F5.0 

6-10 

Pump  flow  goes  to  this  node. 

NJtJNC  (N,2) 

None 

11-15 

Initial  wet  well  volume,  cu.  ft. 

VWELL  (I) 

0.0 

16-20 

Lower  pumping  rate,  cfs. 

PRATE  (1,1) 

0.0 

21-25 

Mid  pumping  rate,  cfs. 

PRATE  (1,2) 

0.0 

26-30 

High  pumping  rate,  cfs. 

PRATE  (1,3) 

0.0 

31-35 

Wet  well  volume  for  mid  rate  pumps  to 
start,  cu.  ft. 

VRATE  (1,1) 

0.0 

36-40 

Wet  well  volume  for  high  rate  pumps  to 
start  pumping,  cu.  ft. 

VRATE  (1,2) 

0.0 

41-45 

Total  well  volume,  cu.ft. 

VRATE  (1,3) 

0.0 

(Last  card  must  have  99999  in  columns  1 to  5) 

1 CARD/FREE  OUTFALL* 

12 

15 

1-5 

Node  for  free  outfall. 

JFREE  (I) 

None 

*0nly  one  connecting  conduit  (Link) 
is  permitted  to  a free  outfall  node 

(Last  card  must  have  99999  in  columns  1 - 5) 
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TABLE  11-2 
(Continued) 


Card 

Group 

Format 

Card 

Columns 

Description 

Variable 

Name 

Default 

Value 

1 CARD/NON-WEIR  TIDE  GATE 

13 

15 

1-5 

Node  at  which  gate  is  located 
(Last  card  must  have  a 99999  in  column  1-5) 

JGATE  (I) 

None 

14 

15,  7F10.0 
F5.0 

. 1-5 

1 no  water  surface  at  outfalls  NTIDE 

2 outfall  control  water  surface  at 
constant  elevation,  A 1 

3 tide  coefficients  provided 

4 program  will  compute  tide  coefficients 

None 

6-10 

First  tide  coefficient 

A 1 

0.0 

11-15 

Second  tide  coefficients 

A 2 

0.0 

16-20 

Third  tide  coefficient 

A 3 

0.0 

21-25 

Fourth  tide  coeffficient 

A 4 

O 

o 

26-30 

Fifth  tide  coefficient 

A 5 

0.0 

31-35 

Sixth  tide  coefficient 

A 6 

0.0 

36-40 

Seventh  tide  coefficient 

A 7 

0.0 

41-45 

Tidal  period  in  hours 

W 

0.0 

REQUIRED  IF  NTIDE  = 4 

15 

315 

1-5 

If  one,  there  are  four  information  points 
program  will  develop  tide  coefficients. 

, KO 

None 

6-10 

Number  of  information  points  (4  if  KO 
above  equals  1). 

NI 

None 

11-15 

If  one,  will  print  information  on  tide 
coefficient  development. 

NCHTID 

0 

REQUIRED  IF  NTIDE  = 4 

16 

8F10.0 

1-10 

Time,  first  information  points. 

TT  (1) 

0.0 

11-20 

Tidal  stage,  at  time  above. 

YY  (1) 

0.0 

21-30 

Time,  second  information  points. 

TT  (2) 

0.0 

Tidal  stage,  at  time  above,  up  to  YY  (2) 

number  of  points  as  defined  by  NI . 


31-40 


0.0 


TABLE  1 1 -2 
(Continued) 


Card 

Card 

Variable 

Oefaul  t 

Group 

Format 

Columns 

Description 

Name 

Value 

Required  if  NJSWil  (in  Card  Group  2) 

17 

1615 

1-5 

First  Input  node  for  card  hydrograph 

JSW  (1) 

None 

6-10 

Second  input  node  for  card  hydrograph, 
up  to  HJSW  nodes. 

JSW  (2) 

None 

18 

8F10.0 

1-10 

Initial  dry  weather  flows  (cfs) 

Q(l) 

None 

11-20 

Initial  dry  weather  velocities  (fps) 

V(l) 

None 

21-30 

Q(2) 

None 

31-40 

V ( 2 ) 

None 

41-  . 

• 

(4  conduits  per  card,  up  to  NTC  conduits. 

Q(NTC) 

None 

Includes  internal  links.) 

V(NTC) 

Hone 

19 

8F10.0 

1-10 

Initial  dry  weather  heads  (ft) 

Y(l) 

None 

11-20 

Y(2) 

None 

Y(NJ) 

None 

Required  if  NJSkCl  (in  Card  Group  2) 

20 

8F10.0 

1-10 

Time,  referenced  to  start  time,  hours 
(decimal ) . 

TEO 

0.0 

11-20 

Flow  rate,  cfs.,  first  input  node, 
JSW  (1). 

QCARD  (L,l ) 

0.0 

21-30 

Flow  rate,  cfs.,  second  input  node, 
JSW  (2),  up  tc  N JSW  nodes. 

QCARD  (L,2) 

0.0 

(Repeat  group  20  cards,  final  time  on 
last  group  20  card  mjst  be  greater 
than  end  time  of  run). 

TABLE  1 1-3 

TRANSPORT  QUALITY  BLOCK  CARD  DATA 


Card 

Group 

Format 

Card 
Col umns 

Description 

Variable 

Name 

De  faul  t 
Val  ue 

1 

1 5A4 

1-60 

(Two  cards) 
Title  cards 

BETA 

None 

2* 

15 

1-5 

Number  of  integration  cycles 

NQCYC 

None 

2F5.0 

6-10 

Length  of  integration  step 
(seconds) 

DELT 

None 

11-13 

Start  time  (hours) 

NHR 

None 

14-15 

Start  time  (minutes) 

MIN 

None 

415 

16-20 

Number  of  nodes  for  detailed 
printout  (20  maximum) 

NPRNT 

None 

21-25 

Cycle  for  start  of  printing 

NSTART 

None 

26-30 

Interval  between  printing  cycles 

INTER 

None 

31-35 

Quality  input  tape 

N21 

None 

A6 

40-45 

Not  used  at  this  time 

NAMER 

15 

46-50 

Hydraulic  input  tape 

N22 

None 

A6 

55-60 

Not  used  at  this  time 

NAMES 

15 

61-65 

Number  of  nodes  for  printer  plot 

NPL0T 

None 

If  NPRNT 

= 0,  skip  cards  3 and  4 

3 

8110 

1-10 

Node  numbers  for  detailed 
printout  (maximum  of  20) 

JQPRT(l) 

None 

11-20 

JQPRT(2) 

None 

21-30 

JQPRT ( 3 ) 

None 

•' 

JQPRT(NPRNT) 

None 

4 

Constituents  to  be  printed, 

0 means  no  plot,  1 means  plot. 

2511 

1 

Ignored 

IC0DL(1) 

None 

2 

Settleable  Solids 

IC0DL(2) 

None 

3 

Suspended  Solids 

IC0DL(3) 

None 

4 

Total  Dissolved  Solids 

IC0DL(4) 

None 

*To  be 

same  as 

Transport  Block 

11-16 
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TABLE  1 1-3 
(Continued) 


Card 

Group  Format 

Card 

Columns 

Description 

Variable 

Name 

Defaul t 
Value 

4 (Continued) 

5 

BOD 

I CODL ( 5 ) 

None 

6 

COD 

I COOL ( 6 ) 

None 

7 

Chlorides 

IC0DL(7) 

None 

8 

S04 

IC0DL(8) 

None 

9 

Grease 

I CODL ( 9 ) 

None 

10 

Total  Coliform 

I CODL ( 10) 

None 

11 

Fecal  Coliform 

ICODL(ll) 

None 

12 

Ammonia 

I CODL (12) 

None 

13 

Organic  N 

ICODL ( 1 3) 

None 

14 

Nitrate  plus  nitrite 

I CODL (14) 

None 

15 

Total  hydrolized  phosphrous 

I CODL (15) 

None 

16 

Orthophosphate 

IC0DL( 16) 

None 

17 

Mercury 

I CODL (17) 

None 

18 

Copper 

I CODL  (18) 

None 

19 

Zinc 

I COOL (19) 

None 

20 

Lead 

IC0DL(20) 

None 

21 

Chromium 

I CODL (21 ) 

None 

22 

Cadmium 

I CODL { 22 ) 

None 

23 

Arsenic 

I COOL ( 23 ) 

None 

NOTE: 

The  arrangement  of  constituents  to  be 
printed  will  apply  to  all  junctions. 
No  changes  can  be  made  between  plots. 

5 

If  NPL0T=0,  skip  cards 

5 and  6 

Junction  cards:  8 values  per  card 

1615 

1-10 

Junctions  which  are  to 
plotted 

be 

JQPRT 
(NPRNT+1 ) 

None 

11-17 
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TABLE  1 1-3 
(Continued) 

Card 

Card 

Variable 

Defaul  t 

Group 

Format  Columns 

Description 

Name 

•Value 

5 

(Continued) 

11-20 

JQPRT(NPRNT+2) 

None 

21-30 

JQPRT(NPRNT+3) 

None 

JQPRNT ( NPRNT+NPLOT ) None 

6 

Constituents  to  be  plotted,  0 

means 

no  plot,  1 means  plot. 

2511  1 

FI  ow 

I CODE  ( 1 ) 

None 

2 

Settleable  Solids 

IC0DE(2) 

None 

3 

Suspended  Solids 

IC0DE(3) 

None 

4 

Total  Dissolved  Solids 

IC0DE(4) 

None 

5 

BOO 

IC0DE(5) 

None 

6 

COD 

I CODE ( 6 ) 

None 

7 

Chlorides 

I CODE ( 7 ) 

None 

8 

S04 

IC0DE(3) 

None 

9 

Grease 

IC0DE(9) 

None 

10 

Total  Col  i form 

I CODE (10) 

None 

11 

Fecal  Col i form 

I CODE ( 1 1 ) 

None 

12 

Ammonia 

I CODE (12) 

None 

13 

Organic  N 

ICODE ( 1 3) 

None 

14 

Nitrate  plus  nitrite 

I CODE ( 14) 

None 

15 

Total  hydrolized  phosphorous 

I CODE ( 1 5) 

None 

16 

Orthophosphate 

1 CODE (16) 

None 

17 

Mercury 

I CODE (17) 

None 

18 

Copper 

I CODE ( 18) 

None 

19 

Zinc 

I CODE (19) 

None 

11-18 
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TABLE  1 1-3 
(Continued) 
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Card 

Card 

Variable 

Defaul  t 

Group  Format 

Col umns 

Description 

Name 

Val  ue 

6 (Continued) 

20 

Lead 

I CODE (20) 

None 

21 

Chromium 

I CODE (21 ) 

None 

22 

Cadmi urn 

I CODE ( 22 ) 

None 

23 

Arsenic 

I CODE ( 23) 

None 

NOTE: 

The  arrangement  of  constituents  to  be 
plotted  will  apply  to  all  junctions. 
No  changes  can  be  made  between  plots. 
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FIGURE  II-!.  SCHEMATIC  OF  'IRFAN  STORMWATER  SETUP  DECK 


FIGURE  1 1-2.  SCHEMATIC  OF  RUNOFF  HOC'  SE^0  DEC 


FIGURE  1 1-3.  SCHEMATIC  OF  TRANSPORT  BLOCK  SETUP  DECK 


FIGURE  1 1-4.  SCHEMATIC  OF  TRANSPORT  QUALITY  BLOCK  SETUP  DEO 


CHAPTER  III 

EXAMPLES  OF  PROGRAM  APPLICATION 


INTRODUCTION 

Two  example  problems  are  presented  in  this  chapter  to  provide 
the  reader  with  a reference  to  an  actual  program  application.  The  first 
example,  Example  A,  shows  the  complete  simulation  of  a typical  area  in 
the  Seattle  vicinity,  the  South  Seattle  Industrial  Park.  This  problem 
demonstrates  the  majority  of  problems  which  one  encounters  in  the 
application  of  the  models  and  is  explained  in  its  entirety  with  full 
reproduc'on  of  the  data  used  for  input  and  the  computer  produced  results. 

The  second  problem,  Example  B,  is  also  fully  presented  and  is 
an  extension  of  the  first  problem  in  that  it  covers  the  situation 
encountered  during  surface  flooding.  This  important  problem  has  been  of 
considerable  concern  in  the  development  of  these  models,  and  it  is  intended 
that  this  example  clear  up  any  remaining  questions  regarding  the  manner 
in  which  this  problem  should  be  handled. 

Any  potential  user  of  the  models  described  in  this  report  is 
urged  to  study  these  example  problems  fully.  By  following  carefully  the 
procedure  set  forth  in  the  examples  many  questions  will  be  answered  and 
the  problem  of  ‘How  do  I get  started?"  should  be  resolved.  Obviously 
question1  will  arise  which  are  not  included  in  these  examples,  but  a good 
understanding  hero  will  go  a long  way  toward  easing  the  startup  pains 
often  encountered  in  the  use  of  computer  programs  as  intricate  and 
comprehensive  as  those  provided  here. 
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As  a further  aid  to  understanding  both  the  models  in  general 
and  the  example  problems  in  particular,  the  reader  is  referred  to 
Figure  1 1 1 - 1 which  is  an  overview  of  a complete  urban  drainage  system. 
While  this  representation  contains  more  elements  than  the  models  can 
accommodate,  one  should  be  able  to  locate  the  areas  to  which  the  various 
model  blocks  are  applicable.  The  Runoff  Block  would  be  used  in  the 
catchment  and  overland  flow  areas  of  the  drainage  system  and  would  be 
used  to  transform  specified  rainfall  patterns  into  expected  hydrographs 
at  locations  on  the  surface  system.  Normally  no  pressure  or  confined 
flow  would  be  included  in  this  description,  and  quality  would  be 
determined  from  the  characteri sties  of  the  watershed. 

The  Transport  and  Transport  Quality  Blocks  would  be  used  to 
extend  the  routing  of  the  surface  hydrographs  into  the  surface  sewer 
system  and  ultimately  to  the  treatment  or  overflow  facilities.  The 
sewer  system  could  be  expected  to  flow  under  both  free  surface  and 
pressure  conditions  and  would  have  a number  of  control  structures  and 
overflow  features  not  found  in  overland  flow.  Momentum  and  backwater 
effects  are  assumed  to  have  more  significance  in  the  routing  of  the 
Transport  Block,  and  this  model  should  be  used  wherever  these  effects 
are  expected  to  dominate  as  compared  to  the  kinematic  wave  effects 
formulated  into  the  Runoff  Block. 


APPLICATION  TO  A TYPICAL  URBAN  WATERSHED 

The  example  problem  is  taken  from  a small  urban  watershed  known 
as  the  South  Seattle  Industrial  Park.  This  area,  which  is  a little  over 
27  acres  in  size,  is  representative  of  a number  of  watersheds  in  the 
Seattle  area  and  is  also  one  of  the  model  calibration  watersheds. 

Figure  1 1 1-2  shows  the  basic  topography  of  the  area,  and  the  necessary 
layout  for  both  the  Runoff  and  Transport  Blocks  have  been  included  on 
the  figure.  The  modeling  elements  associated  with  the  Runoff  Model 
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RAINFALL 


FIGURE  1 1 1-1  URBAN  DRAINAGE  SYSTEM  REPRESENTATION 


are  the  subcatchments  and  gutters,  while  the  Transport  system  features 
are  the  conduits  and  junctions. 

The  example  problem  is  presented  by  first  giving  a general 
narrative  description  of  the  data  and  parameters  which  must  be  supplied 
to  both  the  Runoff  and  Transport  Slocks  and  then  by  indicating  the 
specific  values  chosen  for  the  example  problem.  To  the  extent  possible 
the  initial  narrative  description  has  been  made  general  and  most  comments 
are  applicable  to  problems  beyond  those  of  this  specific  example. 

SETUP  OF  RUUOFF  EL  , 

Use  of  the  surface  flows  portion  of  the  Runoff  Block  requires 
three  basic  steps: 

Step  1 - Geometric  representation  of  the  drainage  basin 

Step  2 - Estimate  of  coefficients 

Step  3 - Preparation  of  data  cards  for  the  computer  program 

.-  -f  Fiszvetizaf- ok 

Discretization  is  a procedure  for  the  mathematical  abstraction 
of  the  physical  drainage  system.  For  the  computation  of  hydrographs  the 
drainage  basin  may  be  conceptually  represented  by  a network  of  hydraulic 
elements,  i.e.,  subcatchments,  gutters,  and  pipes.  Hydraulic  properties 
ft  each  element  are  then  characterized  by  various  parameters,  such  as 
size,  slope,  and  roughness  coefficient. 

Di screti zation  begins  with  the  identification  of  the  drainage 
basin  boundary,  the  location  of  major  inlets,  and  the  selection  of  those 
gutters/pipes  to  be  included  in  the  runoff  model  system.  This  is  best 
shown  by  an  example.  Figure  111-2  shows  the  discretization  of  the  South 
Seattle  Industrial  Park  calibration  drainage  area  which  is  tributary  to 
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a sampling  manhole,  junction  1.  The  Transport  Block  is  used  for  the 
conduit  system  in  this  sample  problem  with  hydrographs  and  pol 1 utographs 
at  inlets  to  the  conduit  system  being  generated  by  the  Runoff  Block. 

It  is  obvious  that  judgment  must  be  exercised  in  setting  up 
the  problems  as  far  as  determining  that  portion  of  the  drainage  system 
that  can  best  be  represented  by  either  the  Runoff  or  the  Transport  Block. 

The  desired  form  of  the  results  usually  indicates  which  model  should  be 
used  for  which  portions  of  a drainage  system.  If  the  objective  is  that 
of  a reconnaissance  level  analysis  of  a drainage  system  and/or  the  impact 
on  the  system  of  alternative  land  use  patterns,  then  the  Runoff  Block 
can  be  used  exclusively  to  represent  the  entire  drainage  system.  For  a 
detailed  analysis  of  the  conduit  drainage  system,  a greater  portion  of 
the  system  should  be  represented  by  the  Transport  Block.  If  backwater 
effects  are  significant  or  if  hydraulic  elements  other  than  conduits  and 
gutters,  such  as  weirs,  are  used,  then  the  Transport  Block  must  be  used. 

Subcatchments  or  unit  watersheds  are  idealized  rectangular  areas 
with  uniform  slope  and  groundcover,  i.e.  asphalt,  concrete,  or  turf.  Each 
subcatchment  has  unique  properties  in  terms  of  slope  and  groundcover.  Thus, 
tiie  roof  of  a house  may  be  represented  by  two  subcatchments  because  the 
water  drains  in  two  different  directions,  even  though  both  units  have  tne 
same  groundcover  and  absolute  ground  slope.  Likewise,  bare  dirt  and 
pavement  can  be  treated  separately  because  of  the  difference  in  groundcover. 

While  the  subdivision  of  any  area  can  be  taken  to  infinitesimal 
detail  in  theory,  computation  time,  data,  and  manpower  requirements  impose 
real  limitations  in  practice.  No  ready  rule  for  the  subdivision  can  be 
offered,  but  if  detailed  analysis  of  the  results  is  desired,  then  correspondi ng 
detail  of  representation  of  the  drainage  basin  should  be  undertaken.  The 
runoff  from  a basin  responds  more  quickly  with  respect  to  changes  in  rainfall 
intensities  as  the  representation  of  the  basin  is  made  more  detailed. 
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Coefficients  and  parameters  necessary  to  characterize  the 
hydraulic  properties  of  a subcatchment  include  surface  area,  width, 
percent  imperviousness,  ground  slope,  roughness  coefficient,  detention 
depth,  infiltration  rate,  and  soil  infiltration  capacity.  Since  real 
subcatchments  are  not  rectangular  areas  experiencing  uniform  overland 
flow,  average  values  must  be  selected  for  computation  purposes. 

For  the  roughness  coefficient,  one  can  use  the  values  given  in 
Table  1 1 1 - 1 , as  suggested  by  Crawford  and  Linsley.  Detention  depths  are 
taken  by  the  program  as  l/16th-inch  for  impervious  areas  and  1/4-inch 
for  pervious  areas,  unless  specified  at  other  values  by  the  user.  The 
infiltration  rate  can  be  specified  by  the  user  from  field  test  data  or 
estimated  from  "standard  infiltration  capacity  curves"  shown  in  Figure  1 1 1 - 3 , 
which  was  produced  by  the  American  Society  of  Civil  Engineers  (ASCE).  A 
limit  on  the  amount  of  rainfall  that  can  infiltrate  into  the  soil  can  be 
placed  on  the  pervious  areas  by  specifying  the  maximum  infiltration  amount 
in  inches.  Twelve  inches  is  used  unless  specified  otherwise  by  the  user. 
Resistance  factors  for  the  pervious  and  impervious  parts  of  a subcatchment 
are  specified  separately  with  default  values  of  .250  and  .013  (Manning's 
n for  overland  flow)  being  taken  in  the  absence  of  other  information. 


The  data  cards  should  be  prepared  according  to  the  input  data 
card  formats  presented  in  Chapter  II  of  this  report.  The  first  step  in 
the  data  preparation  is  the  determination  of  the  number  of  time-steps  to 
be  used  and  the  length  of  each  time-step.  The  time-step  length  is  usually 
2 minutes  but  may  range  from  1 to  30  minutes,  depending  on  the  length  and 
intensity  of  storm  and  the  degree  of  accuracy  required.  The  time  of 
simulation  should  extend  past  the  storm  termination  sufficiently  to  account 
for  the  routing  of  storm  runoff.  Along  with  the  input  of  time-steps,  the 
number  of  unique  hyetographs  for  the  drainage  basin  is  needed. 


1 1 1-7 


r 


TABLE  III-l  ESTIMATE  OF  MANNING'S  ROUGHNESS  COEFFICIENTS 


Ground  Cover 

Manning's  n for 
Overland  Flow 

Smooth  asphalt 

0.012 

Asphalt  or  concrete  paving 

0.014 

Packed  clay 

0.03 

Light  turf 

0.20 

Dense  turf 

0.35 

Dense  shrubbery 
and  forest  litter 

0.4 

Source:  N.  H.  Crawford  and  R.  K.  Linsley,  "Digital 
Simulation  in  Hydrology,  Stanford  Watershed  Model  IV”. 


Source:  American  Society  of  Civil  Engineers, 

Manual  of  Engineering  Practice  No.  37, 
1960. 

FIGURE  III -3  STANDARD  INFILTRATION-CAPACITY 
CURVES  FOR  PERVIOUS  SURFACE 
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The  rainfall  data  cards  are  then  prepared  for  each  hyetograph 
from  rainfall  records  or  are  assumed  if  a hypothetical  test  case  is  neing 
run.  The  time  interval  for  rainfall  input  need  not  be  the  same  as  in  the 
flow  and  quality  portion  of  the  model,  although  the  computation  time  step 
and  the  rainfall  interval  should  be  even  multiples  of  each  other.  The 
major  effort  in  data  preparation  is  forming  a tree-shaped  structure 
drainage  system  and  dividing  the  drainage  basin  into  subcatchments.  The 
drainage  system  data  is  obtained  from  available  maps.  Pipes  smaller  than 
2-3  feet  with  no  backwater  effects,  flow  dividers,  or  lift  stations  are 
usually  designated  as  gutter/pipes  for  computation  by  the  Runoff  Block. 
These  pipes  are  not  connected  to  one  another  by  manholes  but  join 
directly  and  lead  to  an  inlet  for  further  routing  by  the  Transport 
model.  Once  the  drainage  system  is  labeled  with  numbers  (see  Figure 
1 1 1 -2 ) , the  subcatchment  areas  are  formed  reflecting  the  existing 
drainage  system,  ground  cover,  and  land  slope.  Data  cards  are  then 
made  up  for  each  numbered  subcatchment,  defined  by  its  width,  area, 
slope,  percent  imperviousness,  etc.,  along  with  the  gutter/pipe  or 
inlet  into  which  the  flows  are  routed.  Next,  the  gutter/pipe  cards 
are  prepared  giving  the  required  information. 

. A l ‘ * l *€  q A X * J 


Data  to  simulate  surface  quality  are  prepared  at  the  same  time 
as  the  rest  of  the  Runoff  model.  Thus  when  a drainage  basin  is  selected, 
it  should  be  subdivided  into  subcatchments  containing  a single  type  of 
land  use.  Five  land  uses  which  may  modeled  are:  single  family, 

residential,  multi-family  residential,  commercial,  industrial,  and 
undeveloped  or  parklands. 

The  program  may  be  used  with  runoff  from  a design  storm  or  an 
actual  storm.  If  an  actual  storm  is  being  modeled,  the  number  of  dry  days 
prior  to  that  storm  is  determined  from  rainfall  records.  Otherwise,  the 
number  of  dry  days  is  part  of  the  information  associated  with  a design 
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storm.  In  determining  dry  days  from  actual  storms  the  real  number  of 

continuous  antecedent  days  without  rainfall  should  be  increased  to  allow 

for  residual  surface  solids  from  the  earlier  storms.  A suggested  starting 
estimate  for  dry  days  is  the  total  consecutive  antecedent  days  until  the 
sum  of  daily  rainfalls  equals  or  exceeds  1.0  inch.  If  a sizable  storm 
(rainfall  greater  than  0.3  inch)  occurs  within  the  four  days  prior  to  the 
test  storm  the  earlier  storm  should  also  be  modeled.  The  equivalent  dry 
days  should  then  be  calculated  using  the  actual  surface  residual  plus  the 
between-storm  accumulation. 

The  data  needed  on  the  frequency  of  street  cleaning,  which  is 
used  to  reduce  the  basic  buildup  of  pollutants  in  the  watershed  plus  the 
number  of  passes  made  by  the  street  sweeper  can  be  found  from  a public  works 

department.  The  volume  of  liquid  remaining  in  the  catch  basins  may  be  found 

by  analysis  of  the  construction  drawings.  The  quality  of  the  remaining 
liquid  can  be  estimated  or  measured.  Presently  a lack  of  data  has  precluded 
reasonable  estimates  of  the  quality  of  catch  basin  liquids  and  no  parameters 
for  catch  basin  quality  are  included  in  the  program. 

The  data  cards  defining  an  individual  subcatchment  provide  the 
model  with  the  subcatchment,  number,  type  of  land  use,  the  length  of 
gutters  in  the  subcatchment,  and  the  number  of  catch  basins  in  the 
subcatchment.  Land  use  information  may  be  obtained  from  a governmental 
planning  department,  direct  observation,  or  by  other  means.  The  length 
of  gutters  within  each  subarea  may  be  obtained  by  scaling  them  from  a 
street  map.  The  number  of  catch  basins  (gutter  inlets)  per  acre  may  be 
estimated  from  visual  observation  or  obtained  from  a public  works  department. 

The  final  data  cards  for  the  Runoff  Block  are  for  the  printing 
out  of  gutter  and  inlet  hydrographs  and  pollutographs  for  the  user.  On 
the  first  card,  the  total  number  of  gutters  and  inlets  for  which  the 
results  will  be  printed,  and  the  time-step  interval  between  the  printing 
of  the  results  are  specified.  The  gutters  and  inlets  for  which  the  results 
are  to  be  printed  are  then  listed  with  a blank  card  terminating  the  model  run. 
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TRANSPORT  BLOCK 

The  Transport  model  is  usually  set  up  after  the  set  up  of  the 
Runoff  model . Therefore,  decisions  on  the  extent  of  the  drainage  system 
to  be  simulated  by  the  Transport  Block  have  previously  been  made.  The 
data  preparation  is  relatively  straightforward  with  required  information 
coming  from  field  data  and  available  drainage  facilities  maps  and  plans. 

The  South  Seattle  Industrial  Park  calibration  area  example  problem  map. 
Figure  1 1 1 -2,  shows  a good  example  of  a drainage  system  map  with  the 
necessary  data  on  it  for  a simple  pipe  system. 

In  setting  up  the  Transport  model  for  a simulation  run  the  data 
cards  should  be  prepared  according  to  the  input  data  card  formats  presented 
in  Chapter  II  of  this  report.  The  system  control  parameters  are  the  first 
items  to  be  considered  and  these  include: 


1 . 

Simulation  start  time 

2. 

Simulation  end  time 

3. 

Time  interval  (or  step) 

to 

be  used  for  the 

numerical  integration 

4. 

Input  hydrographs  location 

information 

5. 

Control  information  for 

the 

printing  and 

plotting  of  resul ts 

These  parameters  control  the  solution  and  provide  the  information  required 
for  temporally  correlating  the  runoff  hydrographs  computed  in  the  watershed 
model  for  input  to  the  Transport  model.  The  user  specifies  the  junctions 
and  conduits  for  which  he  desires  detailed  results  to  be  printed  and  plotted. 

The  spatial  control  parameters  are  contained  within  the  physical 
system  data  which  are  delineated  next.  The  conduits  and  junctions  are 
represented  by  specifying  the  following  data: 
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1.  Conduit  number 

2.  Junctions  which  the  conduit  connects 

3.  Conduit  type  (circular,  rectangular,  etc.) 

4.  Cross-sectional  area 

5.  Conduit  height 

6.  Maximum  width  of  conduit 

7.  Conduit  length 

8.  Manning's  "n"  value 

9.  Distance  of  conduit  invert  above  nodal 
elevation  if  different  from  zero 

10.  Head  loss  coefficients  at  junctions 

Junctions 

1.  Junction  number 

2.  Ground  elevation  at  the  node 

3.  Invert  elevation  at  the  node 

4.  Net  value  of  constant  inflow  to  the  node,  if  any 

The  net  inflow  value  is  used  for  specifying  the  net  rate  of  infiltration- 
exfiltration  associated  with  the  nodes  as  well  as  any  other  constant  inflow 
or  outflow  that  might  be  associated  with  the  node. 

These  conduit  and  junction  data  specify  the  basic  configuration 
of  the  system.  Other  system  elements  serve  as  transfer  or  diversion 
mechanisms  for  internodal  transfers  within  the  basic  system  defined  by 
the  conduit  and  node  data.  Data  pertaining  to  these  transfer  elements 
are  prepared  next  with  the  following  information  specified  for  each 
transfer  element. 


1 1 1-12 


! 


Orifices 


1.  Junction  containing  the  orifice 

2.  Junction  to  which  the  orifice  discharges 

3.  Cross-sectional  area  of  the  orifice 

4.  Discharge  coefficient,  C0 

W&7,PS 

1.  Junction  containing  the  weir 

2.  Conduit  into  which  the  weir  discharges* 

3.  Weir  type  (transverse,  side  flow,  tide 
gate,  no  tide  gate) 

4.  Vertical  distance  from  node  invert  to  weir  crest 

5.  Vertical  distance  from  the  weir  invert  to  the 
top  of  the  pipe  in  which  the  weir  is  located 

6.  Weir  length 

7.  Discharge  coefficient,  C£ 

8.  Discharge  coefficient  for  surcharged  condition,  C^ 
(recommended  value  of  C^  is  Cw/8) 

Pumps 

1.  Node  in  which  pumping  station  is  located 

2.  Initial  volume  in  the  wet  well 

3.  The  three  pumping  rates  that  specify  the  rule 
curve  for  the  pump  station  (see  Pumps  in  the 
preceding  section) 

4.  The  three  volumes  associated  with  the  three 
pumping  rates 

' » ’ i-t er  of  0 is  specified,  this  signals  that  the  weir  is 

' * f nf  | f , 
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1.  Nodes  which  are  free  outfalls 


1.  Free  outfall  nodes  containing  tide  gates 

. ut  : -is 


1.  Water  surface  elevation/tidal  conditions  at  outfall 

. ■ ’ ' ‘ 'z  ;* ’n  •;*  vs  ar.i  Heads 

1.  Initial  values  of  flow  and  velocity  for  each  conduit 
actually  specified  and  internally  generated 

2.  Initial  value  of  head  for  each  junction  actually 
specified  and  interally  generated 

IHAilSPCRT  QUALITY  BLOCK 


The  Transport  Quality  model  is  used  only  if  the  Runoff  model 
simulates  runoff  quality.  The  data  required  for  the  Transport  Quality 
model  simulation  are  system  control  parameters  and  the  data  cards  should 
be  prepared  according  to  the  input  data  card  formats  presented  in 
Chapter  II  of  this  report.  The  system  control  parameters  include: 

1.  Simulation  start  time 

2.  Simulation  end  time 

3.  Time  interval  (or  step)  to  be  used  for 
the  numerical  integration 


i 
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APPLICATION  TO  SOUTH  SEATTLE  INDUSTRIAL  PARK 


NORMAL  RjNOFF  CONDITIONS 
Runoff  PC' ok 


An  example  problem  is  given  for  the  Runoff  model.  The 
discretization  of  the  South  Seattle  Industrial  Park  calibration  area  is 
used  as  an  example;  the  area  map  is  shown  in  Figure  III-2.  The  listing 
of  the  data  cards  for  the  example  is  given  in  Table  1 1 1 -2 . The  procedure 
for  determining  the  data  for  the  model  representation  of  the  drainage 
basin  is  as  follows: 

1.  Define  the  drainage  basin  on  a contour  map  of 
sufficient  scale  to  readily  show  the  drainage 
system  for  the  example  problem;  in  this  case 

it  is  the  area  tributary  to  the  sampling  manhole 
which  is  designated  as  junction  1. 

2.  Determine  the  rainfall  hyetograph  for  the  simulation 
run  for  the  basin--for  the  example  problem,  actual 
rainfall  measurements  from  a nearby  gage  were  used. 

The  length  of  time  of  the  hyetograph  plus  sixty 
minutes  added  at  the  end  of  the  rainfall  to  allow 
for  routing  of  the  storm  runoff  was  used  for  the 
length  of  the  simulation. 

3.  Identify  the  drainage  system  to  be  used  in  the 
Runoff  model --for  the  example  problem,  surface 
gutters  flowing  to  drainage  inlets  were  used. 

4.  Identify  the  inlets  for  the  Runoff  model  drainage 
system--for  the  example  problem,  the  drainage 
inlets  were  designated  as  inlets  for  the  model  to 
store  hydrographs.  These  are  numbered  1 through  7. 
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TABLE  1 1 1-2 

RUNOFF  BLOCK  EXAMPLE  PROBLEM  INPUT  DATA  CARDS 


A 


5.  Identify  the  tributary  area  and  drainage  system  to 

each  inlet--for  the  example  problem,  eleven  subcatchments 
and  eleven  gutters  were  used. 

6.  Determine  values  of  parameters  to  represent  drainage 
system  gutters/pipes--for  example  problem,  lengths, 
street  section,  slopes  were  measured  and  calculated. 

An  "n"  value  of  0.015  was  assumed  for  the  gutters. 

7.  Determine  values  for  parameters  and  coefficients  to 
represent  the  subcatchments--for  the  example  problem, 
area,  width  of  subcatchment  perpendicular  to  overland 
flow,  percent  imperviousness,  average  ground  slope 
were  determined  for  each  subcatchment.  Coefficients 
for  resistance  factors,  retention  storage,  and 
infiltration  were  assumed  to  be  those  default  values 
of  the  program. 

8.  Determine  if  quality  constituents  are  to  be  modeled. 

If  so,  find  the  number  of  dry  days  prior  to  precipitation 
to  be  modeled,  the  street  sweeping  frequency,  the  number 
of  passes  per  cleaning  and  the  average  catch  basin 
volume.  For  the  example  problem  actual  rainfall  records 
indicated  that  there  were  six  dry  days  prior  to  the 
storm;  from  the  public  works  department  it  was  learned 
that  the  streets  were  swept  once  a week  with  one  pass 
per  cleaning;  the  average  volume  of  a catch  basin  was 
assumed  to  be  25  cubic  feet. 

9.  Within  each  subcatchment  determine  the  length  of 
gutters  and  the  number  of  catch  basins--for  the 
example  problem,  the  length  of  gutter  was  measured 
and  there  is  one  catch  basin  per  subcatchment. 

10.  Identify  those  gutters  for  which  printout  of  hydrographs 
and  pol 1 utographs  are  desired  and  the  time  interval  for 
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the  results  to  be  printed--for  the  example  problem,  two 
gutters,  101  and  104,  were  selected  for  printout  and  an 
interval  of  10  minutes  (5  time  steps)  was  selected. 

11.  Punch  the  cards  and  execute  the  program. 

A reproduction  of  all  computer  output  from  the  Runoff  Block  for  the 
example  problem  is  presented  in  Tables  1 1 1-3  through  1 1 1 -9  and  Figures 
1 1 1-4  and  1 1 1-5. 


Tables  1 1 1 - 3 through  1 1 1 -6  show  an  echo  of  the  input  data  and 
Tables  1 1 1 - 7 , 1 1 1 -8  and  III -9  indicate  the  total  mass  balance  and  the 
time  history  of  quality  and  quantity  at  gutters  101  and  104,  respectively. 
Figure  1 1 1 -4  indicates  the  shape  of  the  input  hyetograph  and  Figure  1 1 1 - 5 
plots  the  shape  of  the  surface  flow  hydrograph  at  junction  1.  Figures  III-6 
to  1 1 1 -9  are  examples  of  the  hydrograph  and  pollutograph  plots  produced 
by  the  Runoff  Block.  Figure  III -6  is  a plot  of  the  time  history  of  flow 
at  gutter  101  while  Figures  1 1 1 - 7 to  III -9  are  the  time  histories  of 
three  quality  parameters  at  the  same  location. 
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NRE  URBAN  RUNOFF  SIMUALT ION  MODEL 
SURFACE  drainage  QUANTITY  AND  QUALITY 
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PRINTED  TAPE  AND  RAINFALL  INPUT  DATA 
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TABLE  1 1 1 -4 

PRINTED  CONDUIT  AND  WATERSHED  INPUT  DATA 
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TABLE  1 1 1-6 

QUALITY  SIMULATION  PARAMETERS 
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TABLE  1 1 1-8 

LOCATION  101  FLOW  AND  QUALITY  SUMMARY 
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TABLE  1 1 1-9 
(Continued) 


"AINPALL  hyetoorapm 


FIGURE  1 1 1-4  EXAMPLE  PROBLEM,  HYETOGRAPH  PLOT 


ON  131NI  MO  xJuno  "013 


FIGURE  1 1 1-6  TIME  HISTORY  OF  GUTTER  FLOW 


SUSPENDED  SOLIDS 


FIGURE  1 1 1-7  TIME  HISTORY  OF  GUTTER  QUALITY— SUSPENDED  SOLIDS 
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The  next  example  is  for  the  Transport  model.  The  South  Seattle 
Industrial  Park  calibration  area's  storm  drain  conduit  system  is  used  as 
an  example  and  is  shown  in  the  area  map  in  Figure  1 1 1 -2 . The  hydrographs 
and  pol 1 utographs  generated  by  the  Runoff  Block  simulation  described  in 
the  preceding  example  problem  are  used  as  input  data  to  the  Transport 
model  example.  The  procedure  for  determining  the  data  for  the  model 
representation  of  drainage  system  is  as  follows: 

1.  Determine  the  elements,  conduits  and  junctions  of 
the  drainage  system — for  the  example  problem,  there 
are  six  conduits  and  seven  junctions.  The  conduits 
should  be  of  a relatively  consistent  length  to  help 
the  stability  of  the  solution  technique  and  to  keep 
the  integration  period  as  long  as  possible.  The 
length  of  the  shortest  conduit  determines  the  maximum 
integration  period.  The  time  step  should  be  short 
enough  so  that  the  water  flow  in  the  conduit  cannot 
pass  entirely  through  the  length  of  the  conduit. 

An  integration  period  of  10  seconds  is  used  in 
the  example. 

2.  Specify  the  control  parameters  for  the  run--for  the 
example  problem,  the  simulation  start,  length  and 
tape  inputs  must  correspond  to  those  of  the  Runoff 
model.  Since  the  example  problem  is  small  in  size, 
the  results  for  all  the  conduits  and  junctions  are 
printed  out. 

3.  Specify  the  parameters  for  the  representation  of  the 
conduits  and  junctions--for  the  example  problem, 
conduit  data  and  junction  data  were  read  or  measured 
from  the  area  map.  An  "n"  value  of  0.015  was  assumed 
for  the  conduits. 


1 1 1-34 


I 


4.  Specify  the  parameters  for  the  representation  of 
orifices,  weirs,  pumps,  and  outfalls.  For  the 
example  problem,  only  one  free  outfall  was 
specified. 

5.  Specify  the  initial  dry  weather  flows  and  heads  for 
the  actual  and  internally  generated  conduits  and 
junctions.  For  the  example  problem  three  blank 
cards  for  the  initial  heads  are  included  in  the 

data  deck  (i.e.,  initialize  flows  and  heads  to  zero). 

A listing  of  the  actual  input  data  cards  is  given  in  Table  III-10.  The 
results  of  transport  model  example  problem  are  presented  in  Tables 
III-ll  t 11-15  and  Figures  1 1 1 - 1 0 and  1 1 1 - 1 1 . 

.es  III-ll  and  1 1 1 - 1 2 are  an  echo  check  of  the  card  input 
data,  while  1 1 1 - 1 3 shows  the  input  hydrographs  produced  as  output  from 
the  Runoff  Block.  Table  1 1 1 - 1 4 gives  the  position  of  the  hydraulic  grade 
line  at  selected  junctions  in  the  system  and  Table  1 1 1 - 1 5 provides  a 
time  history  of  flow  and  velocity  in  various  system  conduits.  Typical 
printer  plots  of  relevant  system  information  are  shown  in  Figures  1 1 1 - 1 0 
and  III-ll.  Figure  1 1 1 - 1 0 indicates  the  water  surface  of  junction  2, 
while  the  time  history  of  flow  in  conduit  201  is  displayed  in  Figure  III-ll. 
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TABLE  III-10 

TRANSPORT  BLOCK  EXAMPLE  PROBLEM  DATA  CARDS 
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TIME  HISTORY  OF  HYDRAULIC  GRADE  LINE 


30UT"  SEATTLE  INDUSTRIAL  HARK  - TRANSHmT  MMEl  - RUN  ST.  1 »A  T£R  RESOURCES  ENGT’.r.ERS,  INC 

E»A«PLE  PROMLE**  USING  STORM  Of  MARCH  10,1973  AT  DIAGONAL  aAlNUT  C'<EEK,  CALIH  <NIA 

STORM  DRaINACF  MCOtl. 


TABLE  1 1 1 - 1 4 
(Continued) 


SOUTH  SEATTLE  INDUSTRIAL  PARK  - TRANSPORT  HODEL  » RUN  NO,  1 • WATER  RESOURCES  ENGINEERS,  INC 

EXAMPLE  PROBLEM  USING  STORM  OE  march  tO> 1973  AT  DIAGONAL  wALNUT  CREEK,  CALIFORNIA 

STORM  DRAINAGE  MODEL 
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TABLE  III-15 

TIME  HISTORY  OF  FLOWS  AND  VELOCITIES 


SOUTH  SEATTLE  INDUSTRIAL  PAR*  - TRANSPORT  MODEL  - HUN  NO.  1 NATEH  RESOURCES  ENGlNFrnS,  INC 

EXAMPLE  PMODLEM  USING  STORM  OF  MARCH  16.  1973  AT  DIAGONAL  nalNUI  CHIC*.  CALITOmMA 

S1URM  UHAINAGE  MODEL 


FIGURE  III-11  TIME  HISTORY  OF  CONDUIT  FLOW 
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An  example  problem  is  now  given  for  the  Transport  Quality 
Block.  The  South  Seattle  Industrial  Park  calibration  area  drainage  system 
is  used  as  an  example.  The  pollutographs  generated  from  the  Runoff  Block 
are  transferred  through  the  Transport  Block  and,  along  with  the  hydraulic 
results  from  the  Transport  model,  are  used  as  input  data  to  the  Transport 
Quality  model.  The  Transport  Quality  model  does  not  require  additional 
specific  data  except  for  run  control  parameters.  The  data  control 
parameters  must  correlate  to  those  used  in  the  Transport  Block. 

The  data  card  listing  is  presented  in  Table  1 1 1 - 1 6 and  examples 
of  program  output  are  presented  in  Tables  1 1 1 - 1 7 through  II I -22 . 

Tables  1 1 1 -1 6 and  1 1 1 - 1 7 recap  the  input  data  and  Tables  1 1 1 - 1 8 through 
III -22  indicate  typical  time  histories  of  quality. 

If  the  user  desires  a complete  history  of  quality  behavior  the 
constituents  are  plotted  in  the  following  order: 


1 . 

Settleable  Sol  ids 

12. 

Organic  Nitrogen 

2. 

Suspended  Solids 

13. 

Nitrate  + Nitrite 

3. 

Total  Dissolved  Solids 

14. 

Total  Phosphate 

4. 

5-day  BOD 

15. 

Ortho  Phosphate 

5. 

Chemical  Oxygen  Demand 

16. 

Mercury 

6. 

Chlorines 

17. 

Copper 

7. 

Sul  fates 

18. 

Zinc 

8. 

Grease 

19. 

Lead 

9. 

Total  Coliform 

20. 

Chromi urn 

10. 

Fecal  Coliform 

21 . 

Cadmi urn 

11  . 

Ammonia 

22. 

Arsenic 

Typical  printer  plots  for  the  mass  emission  rates  of  settleable 
solids  and  BOD  are  indicated  for  junction  1 in  Figures  1 1 1 - 1 2 and  1 1 1 - 1 3 , 
respecti vely. 
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TABLE  1 1 1 -1 6 
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TABLE  1 11-18 

COMPUTED  TIME  HISTORY,  SETTLEABLE  SOLIDS 
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TABLE  III-19 

COMPUTED  TIME  HISTORY,  SUSPENDED  SOLIDS 
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TABLE  111-22 

COMPUTED  TIME  HISTORY,  COD 


FIGURE  1 1 1-12  COMPUTED  TIME  HISTORY,  SETTLEABLE  SOLIDS 


JUNCTION  no 
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SURFACE  FLOODING  CONDITIONS 
Fun  1 - Induced  Flooding 

An  example  problem  of  the  Transport  model's  capability  of 
surface  flooding  is  given.  The  South  Seattle  Industrial  Park  storm  drain 
pipe  system,  as  shown  in  Figure  III-2,  is  used  as  the  example  problem. 

Two  runs  are  made,  one  to  show  flooding  and  a second  to  show  a solution 
to  the  flooding  problem.  For  the  flooding  example.  Run  1,  the  capacity 
of  the  existing  drainage  system  is  overwhelmed  with  an  input  hydrograph 
placed  at  junction  7 which  creates  flooding  at  two  junctions,  6 and  7. 

The  data  listing  and  results  are  shown  in  Tables  1 1 1 -23  through  III -28 
and  Figures  1 1 1 - 1 4 through  1 1 1 - 1 7 . 

Table  1 1 1 -2  3 shows  input  card  listings,  while  Tables  1 1 1 - 24 
and  1 1 1 - 25  display  the  program's  echo  check  of  the  input  data.  Table  1 11-26 
indicates  the  output  produced  in  a flooding  situation  and  shows  the  components 
of  flow  entering  each  flooded  junction  plus  the  total  flooded  volume  in  cubic 
feet.  Table  1 1 1 - 27  shows  the  history  of  the  hydraulic  grade  line  and  Table 
III -28  indicates  the  time  history  of  flows  in  the  system  conduits. 

In  Figures  1 1 1 - 1 4 and  1 1 1 - 1 5 we  see  plots  of  water  surface 
elevations  at  junctions  2 and  7,  and  Figures  1 1 1 - 1 6 and  1 1 1 -1 7 provide 
graphical  displays  of  the  flow  in  conduits  201  and  204.  These  results 
show  that  the  maximum  flooding  volume  at  junction  7 is  1673  cubic  feet 
and  at  junction  6 the  maximum  is  53  cubic  feet.  Using  the  maximum  amount 
of  flooding  and  a detailed  contour  map  or  volume-elevation  table  at  a 
node,  the  user  can  determine  the  extent  of  flooding.  With  the  volume 
of  flooding  known,  the  user  can  formulate  plans  to  eliminate  flooding 
or  contain  the  flooding  in  a prescribed  holding  basin. 
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TABLE  1 11-27  TIME  HISTORY  OF  HYDRAULIC  GRADE  LINE 
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TABLE  1 1 1-28  TIME  HISTORY  OF  CONDUIT  FLOWS  AND  VELOCITIES 
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FIGURE  1 1 1 - 1 4 TIME  HISTOR'  OF  WATER  SURFACE  ELEVATION,  FLOODING  RUN  I 
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FIGURE  1 1 1 - 1 5 TIME  HISTORY  OF  WATER  SURFACE  ELEVATION,  FLOODING  RUN  1 


FIGURE  1 1 1 -16  TIME  HISTORY  OF  CONDUIT  FLOW , FLOODING  RUN  1 


Fun  2 - Possible  Solution  to  Flooding  Problem 

Run  2 of  the  flooding  example  illustrates  a manner  in  which 
the  Transport  model  can  simulate  a solution  to  the  flooding  problem 
that  occurs  in  Run  1.  Flooding  occurs  at  junction  6 because  pipe  205 
has  a smaller  capacity  than  the  upstream  pipe  206.  For  Run  2,  pipe  206 
is  increased  in  size  from  1.5  feet  to  1.75  feet  in  diameter.  This 
increase  will  eliminate  flooding  at  junction  6 for  Run  2. 

It  is  assumed  for  the  example  problem  that  any  flooding  at 
junction  7 will  be  allowed  as  long  as  it  is  maintained  between  the  curbs 
of  the  street.  The  ground  elevation  of  7 feet  at  junction  7 in  Run  1 is 
assumed  to  be  the  gutter  elevation  and  therefore  the  top  of  curb  elevation 
is  assumed  to  be  7.50  feet.  For  Run  2 the  ground  elevation  at  junction  7 
is  raised  to  7.50  feet.  Also,  connected  to  junction  7 are  three  conduits 
to  represent  storage  available  in  the  street  section.  Conduit  401  represents 
the  street  section  as  50  feet  wide,  0.5  feet  deep  and  200  feet  long.  Conduit 
301  is  a short,  fictitious  conduit  to  connect  conduit  401  and  junction  101 
to  junction  7.  Conduit  402  is  another  fictitious  conduit  that  connects 
conduit  401  to  a system  outfall.  For  the  example  problem,  junction  1 is 
designated  as  the  downstream  end  of  conduit  402  for  ease  of  data  preparation. 
See  Figure  1 1 1 - 1 8 for  a graphic  display  of  the  street  section  representation 
and  Table  1 1 1 - 29  for  the  data  listing  of  Run  2. 

The  results  of  Run  2 are  shown  in  Tables  1 1 1 - 30  through  1 1 1 - 33 
and  Figures  1 1 1 - 1 9 through  1 1 1 - 22 . As  before,  Tables  1 1 1 - 30  and  1 1 1 - 31 
give  an  echo  check  of  the  input  data,  while  Table  1 1 1 - 32  shows  the  position 
of  the  hydraulic  grade  line  and  1 1 1 - 33  the  flow  in  various  system  conduits. 
Figures  1 1 1 - 1 9 and  1 1 1 - 20  plot  the  water  surface  elevation  at  junctions  2 
and  7 while  Figures  1 1 1 - 2 1 and  II 1-22  display  time  histories  of  conduit  flow. 

There  is  no  flooding  in  Run  2 for  the  revised  system  using  the 
same  input  hydrograph  as  in  Run  1.  The  water  surface  elevation  at  junction  7 
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rises  to  a peak  of  7.20  feet  and  recedes  with  the  reduction  of  flow  from 
the  input  hydrograph.  The  printout  of  the  flow  in  conduit  301  shows  water 
filling  the  street  section  initially  and  then  later  draining  from  the 
street  section.  A negative  flow  in  a conduit  means  that  the  flow  is  from 
the  junction  with  a lower  invert  elevation  to  a junction  with  a higher 
invert  elevation;  in  the  example  this  is  from  junction  7 to  junction  101. 
Zero  flow  is  shown  for  a conduit  401  and  402  because  junction  102  is  not 
wetted  during  the  example  problem. 

Before  leaving  the  subject  of  surface  flooding,  one  additional 
item  should  be  mentioned:  stationary  surface  flooding.  If  a gutter  or 

pipe  becomes  surcharged  in  a low  area  water  will  leave  the  system  and 
it  is  often  of  interest  to  know  the  depth  and  extent  of  the  expected 
flooding.  To  handle  this  problem  the  user  need  only  specify  a "dead 
end"  storage  element  at  the  node  in  question  and  the  model  will  simulate 
the  inundation  and  draining  of  the  surface  area  in  accordance  with  the 
capacity  and  discharge  of  the  basic  sewer  system.  One  can  specify  the 
stage-volume  relationship  of  the  flooded  area  by  the  appropriate  selection 
of  one  or  more  equivalent  trapezoidal  sections  and  the  model  will  produce 
the  time  history  of  water  surface  (surcharge)  at  the  given  nodes.  With 
a little  experience  the  user  will  find  this  method  of  estimating  surface 
flooding  both  easy  and  efficient  as  no  new  types  of  data  or  changes  in 
the  basic  computer  codes  is  required. 
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TABLE  1 1 1-29 

DEMONSTRATION  FLOODING  PROBLEM  INPUT  DATA,  RUN  2 
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PRINTED  CONTROL  DATA 
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FIGURE  1 11-22  TIME  HISTORY  OF  CONDUIT  FLOW,  FLOODING  RUN  2 


CHAPTER  IV 

MODEL  CALIBRATION  PROCEDURE 


INTRODUCTION 

In  order  to  fill  the  need  for  prototype  data  representative 
of  the  Seattle  Metropolitan  area,  a program  of  data  collection  and  model 
calibration  was  undertaken  in  connection  with  this  project.  The  details 
and  results  of  the  data  collection  program,  conducted  by  Metro,  can  be 
found  in  Metro's  "Water  Quality  and  Quantity  Monitoring  Program  for  the 
RIBCO  Urban  Runoff  and  Basin  Drainage  Study."  It  is  the  intent  of  this 
adapter  to  outline  the  procedure  used  in  the  model  calibration  process 
and  to  indicate  the  results  of  our  calibration  effort.  The  complete  set 
of  cal ibration  computations  are  both  tedious  and  voluminous  and  thus 
only  examples  of  the  calibration  calculations  are  included  in  this  chapter 
to  give  the  reader  an  understanding  of  the  procedures  employed  to  develop 
the  model  coefficients  for  the  Seattle  Metropolitan  area. 


QUANTITY  CALIBRATION 

Model  calibration  is  defined  to  be  the  process  of  determining 
local  and  appropriate  values  for  the  empirical  coefficients  which  appear 
in  the  various  models.  It  should  be  recognized  that  the  calibration 
process  itself  is  a somewhat  empirical  process,  and  that  no  quantifiable 
measure  is  made  of  how  well  the  models  have  been  calibrated.  In  the 
paragraphs  which  follow  the  adequacy  of  the  calibration  is  judged  by 
visual  comparisons  of  model /prototype  behavior  and  the  point  of  "best" 
fit  located  by  judgmental  determination. 


IV-I 


The  calibration  process  was  primarily  concerned  with  the 
simulation  of  surface  runoff.  The  calibration  was  performed  using  data 
on  rainfall,  runoff  and  runoff  quality  which  was  collected  for  seven 
calibration  areas.  The  seven  areas  varied  in  size  from  24  acres  to  over 
600  acres.  The  areas  represented  four  different  land  uses:  1)  single 

family  residential,  2)  multi-family  residential,  3)  industrial,  and 
4)  commercial.  Computer  simulation  runs,  attempting  to  duplicate  the 
gaged  hydrographs  and  pol lutographs,  were  performed  on  the  seven  areas 
using  a number  of  gaged  rainfall  events.  Table  IV-1  summarizes  the 
areas  used  for  calibration  in  this  project. 

Rainfall  gages  were  installed  in  or  near  each  of  the  seven 
calibration  areas.  The  rainfall  hyetographs  from  these  gages  were  used 
directly  in  making  calibration  runs.  When  a rain  gage  was  not  located 
within  a calibration  area  the  hyetograph  from  the  nearest  rain  gage 
representative  of  the  gaged  runoff  hydrograph  was  used  for  the  calibration 
runs.  The  rainfall  hyetograph  is  the  most  important  factor  in  determining 
the  shape  and  magnitude  of  a simulated  hydrograph. 

The  parameters  selected  for  the  representation  of  the  subcatchment 
were  carefully  chosen  to  depict  the  runoff  from  the  calibration  area.  The 
percent  imperviousness  of  the  subcatchment  is  a critical  parameter  for 
calibrating  the  quantity  of  runoff  from  a subcatchment.  The  other  parameters 
of  the  subcatchment  affect  the  runoff  but  none  as  importantly  as  the  percent 
imperviousness. 

Once  the  rainfall  hyetographs  and  the  subcatchment  parameters 
were  determined,  simulation  runs  on  the  calibration  areas  were  performed 
to  best  duplicate  the  gaged  runoff  hydrograph.  Success  in  duplicating 
the  gaged  hydrographs  with  respect  to  total  volume  of  runoff  and  the  shape 
of  the  hydrograph  was  obtained  for  most  areas  but  not  for  every  rainfall 
event.  The  lack  of  accurate  rainfall  and  runoff  data  for  some  areas  made 
the  quantity  calibration  of  those  areas  more  difficult  to  perform. 


Calibration  Area 


Land  Classification 


Area  in  Acres 


1. 

View  Ridge  1 

Single  Family 

Resi dential 

630 

2. 

View  Ridge  2 

Multi-Family  1 

Residential 

105 

3. 

South  Seattle 
Ind.  Park 

Industrial 

27.5 

4. 

Southcenter 

Commercial 

24 

5. 

Lake  Hills 

Single  Family 

Residential 

139.9 

6. 

Hignland 

Single  Family 

Residential 

107.4 

7. 

Central 

Business  District 

Commercial 

27.8 

TABLE  IV-1 
CALIBRATION  AREAS 
SEATTLE  METROPOLITAN  AREA 


The  results  of  the  runoff  model  calibration  process  indicate 
good  representation  of  the  various  types  of  drainage  basins  associated 
with  the  seven  calibration  areas.  Excellent  results  of  simulating  gaged 
hydrographs  were  obtained  for  both  the  largest  calibration  areas.  View 
Ridge  1 and  View  Ridge  2 (735  acres  total)  and  the  smaller  areas,  Central 
Business  District  and  South  Seattle  Industrial  Park  (25  acres  each). 

Typical  of  the  calibration  results  are  the  data  shown  in  Figure  I V- 1 . 

This  figure  shows  the  observed  and  modeled  outflow  hydrographs  for  the 
View  Ridge  1 and  Central  Business  District  calibration  areas  for  the 
storm  of  March  10,  1973.  The  agreement  between  model  and  prototype  is 
judged  satisfactory  for  these  areas,  with  total  runoff  and  general 
hydrography  shape  in  reasonably  close  agreement  for  both  areas.  Of 
particular  importance  in  obtaining  a good  simulation  is  an  accurate 
rainfall  hyetograph  which  represents  the  actual  rainfall  event  that 
occurred  over  the  calibration  area.  This  was  particularly  evident  for 
the  25-acre  areas.  A rain  gage  located  even  half  a mile  to  a mile  from 
an  area  does  not  necessarily  represent  that  area's  rainfall,  especially 
the  rainfall  intensity.  The  calibration  areas  of  Southcenter,  Lake  Hills, 
and  Highlands  had  a limited  amount  of  data  that  was  suitable  for  the 
calibration  process  due  to  1)  the  rain  gage  was  not  representative  for 
the  area,  2)  the  runoff  hydrograph  and  the  rainfall  records  did  not 
correspond  to  each  other,  3)  the  data  was  not  available  as  a result  of 
recorder  malfunction.  In  every  calibration  area  at  least  one  good 
comparison  between  observed  and  simulated  hydrographs  was  obtained  with 
the  exception  of  the  Highlands  area.  The  result  of  the  calibration 
exercise  is  that  the  user  can  be  confident  that  the  Runoff  model  reasonably 
represents  the  runoff  process  for  watersheds  of  varying  sizes  and  land  uses 
in  the  Seattle  Metropolitan  Area,  assuming  good  rainfall  data  can  be 
suppl ied  as  an  input. 
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RAINFALL  HYETOORAPH  AT  U OF  W RAINFALL  H YE  TOG  R a PH  AT  DENNY  WAY 


FIGURE  IV-1  COMPARISON  OF  GAGED  AND  SIMULATED  HYDROGRAPHS 
FOR  THE  STORM  OF  MARCH  10,  1973 


QUALITY  CALIBRATION 


:y.yj tituent  relationships 


Subsequent  to  obtaining  satisfactory  quantity  simulations  as 
determined  by  a comparison  of  simulated  and  observed  hydrographs,  the 
quality  calibration  was  begun.  Quality  data  was  obtained  for  a few 
select  storms  at  the  gaging  manholes  of  the  calibration  areas.  The 
quality  data  consisted  of  the  analysis  of  the  manhole  grab  samples  for 
a total  of  23  water  quality  parameters.  The  grab  samples  were  taken 
at  a minimum  interval  of  15  minutes  for  the  duration  of  the  rainfall 
event.  An  example  of  the  water  quality  data  is  shown  in  Table  IV-2. 

Using  the  quality  data  and  the  gaged  hydrographs,  quality  constituent 
relationships  for  each  of  the  seven  calibration  areas  were  developed  by 
determining  the  number  of  pounds  in  trie  runoff  of  each  constituent  Tor 
several  rainfall  events,  and  then  determining  the  average  relative  amounts 
with  respect  to  total  solids  (settleable  solids  + suspended  solids  + total 
dissolved  solids).  Mass  emissions  rates  in  pounds  per  second  for  all 
constituents  were  determined  by  multiplying  the  constituent  concentration 
by  the  gaged  flow  at  the  time  of  the  grab  sampling.  The  mass  emission 
rates  were  then  integrated  over  the  duration  of  the  storm  to  arrive  at 
the  total  poundage  of  each  constituent  assumed  to  have  washed  off  of  the 
calibration  area  during  the  rainfall  event.  Table  IV-3  is  an  example 
summary  of  the  aforementioned  calculations  showing  the  amount  of  each 
constituent  and  their  relationships  to  the  amount  of  total  solids  for 
three  rainfall  events  at  the  South  Seattle  Industrial  Park  calibration 
area.  The  percentages  of  the  constituents  based  on  the  total  solids 
were  then  averaged  to  develop  the  23  constituent  relationships  based  on 
land  use.  Since  there  were  only  a few  storms  that  were  accurately 
monitored  for  complete  quantity  and  quality  data,  judgment  was  used  to 
eliminate  obvious  extreme  values  which  would  erroneously  affect  the 
average.  These  average  values  were  then  incorporated  into  the  model 
as  the  variable  QFACT  and  are  the  values  that  make  the  model  specific 
to  the  Seattle  area.  The  values  of  QFACT  are  shown  in  Table  IV-4. 


I V- 6 


i 


OUALIT*  DATA  F3i«  SOUTH  SEATTLE  INDUSTRIAL  PARK  ■»  STORM  OF  MARCH  16,1973 
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TABLE  IV-2  EXAMPLE  OF  MONITORED  WATER  QUALITY  DATA 


Storm  of  March  10  Storm  of  March  16  Storm  of  June 
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TABLE  I V-3 

STORMWATER  QUALITY  CONSTITUENT  AMOUNTS  FOR  THE 
SOUTH  SEATTLE  INDUSTRIAL  PARK  CALIBRATION  AREA 


Land  Use  Classification 


Consti tuents 

Single  Family 

Multi-Family 

Commercial 

Industrial 

Undeveloped 
and  Parks 

Set  Solids 

180 

280 

240 

176 

153 

Sus  Solids 

340 

250 

445 

249 

323 

TDS 

480 

470 

315 

575 

525 

BOO 

35.7 

53.1 

35.9 

14.0 

11.0 

COD 

338 

427 

300 

240 

222 

Chlorides 

26.3 

29.8 

22.2 

37.9 

54.7 

S04 

31.7 

48.4 

49.8 

55.7 

92.2 

Grease 

45.1 

78.6 

30.2 

35.4 

28.0 

Tot  Col i form 

7.1  x 105 

2.9  x 106 

2.9  x 104 

6.0  x 104 

6.2  x 104 

Fecal  Col i form 

CD 

ro 

X 

o 

1.9  x 105 

8.2  x 102 

7.7  x 103 

2.0  x 104 

NH3 

0.66 

2.89 

0.93 

0.61 

0.37 

Org.  Nit. 

5.00 

0.61 

4.05 

2.82 

4.13 

N03  + N02 

2.51 

4.51 

2.27 

1.70 

7.29 

Total  Hyd  P 

1.06 

1.08 

0.65 

0.69 

1.19 

Ortho  P04 

0.25 

0.47 

0.19 

0.18 

0.36 

Hg 

0.0053 

0.0030 

0.0012 

0.0009 

0.0019 

Cu 

0.29 

0.17 

0.21 

0.19 

0.78 

Zn 

0.40 

0.33 

0.78 

0.98 

0.26 

Pb 

1.21 

1.57 

1.07 

0.51 

0.68 

Cr 

0.056 

0.052 

0.069 

0.032 

0.061 

Cd 

0.025 

0.022 

0.015 

0.015 

0.030 

As 

0.178 

0.139 

0.117 

0.171 

— 

’Loadings  in  mg/gram  of  Total 

Solids  or  organisms/gram  for 

Col i form. 

TABLE  I V- 4 


STORMWATER  QUALITY  CONSTITUENTS  RELATIONSHIPS1 


POLLUTANT  BUILDUP 

The  rate  of  pollutant  buildup  on  the  calibration  areas  was  not 
measured  in  the  calibration  process  and  must  be  considered  an  unknown  at 
present.  The  model  uses  a straight  line  accumulation  of  dust  and  dirt 
(assumed  equivalent  to  total  solids  in  runoff)  based  on  a variable  amount 
of  pounds  in  dust  and  dirt  per  day  per  100  feet  of  curb,  depending  on 
the  type  of  land  use.  The  original  EPA  Storm  Water  Model  used  pollutant 
loading  rates  determined  in  a separate  study  in  the  Chicago  area.  These 
rates  did  not  seem  applicable  for  the  Seattle  area  in  light  of  data 
obtained  from  additional  studies  recently  performed  in  the  local  area. 
Based  on  the  EPA  report  titled,  "Water  Pollution  Aspects  of  Street 
Surf  :e  Contaminants,"  dated  November  1972,  an  average  loading  factor 
of  0.72  pounds  per  day  per  100  feet  of  curb  was  adopted.  This  rate 
of  0.72  is  less  than  one  half  of  the  average  loading  used  in  the 
original  model.  The  loading  rates  for  the  five  land  use  classifications 
were  then  determined  based  on  the  original  model  rates  and  the  average 
Seattle  loading  rate  as  shown  in  Table  IV-5. 


Land  Use  Classification 

EPA  Stormwater  Model 
Loading  Rate 

RIBC0  Loading 
Rate 

Overal  1 

1.5 

0.72 

Single  Family  Residential 

0.7 

0.4 

Multi-Family  Residential 

2.3 

1 .1 

Commercial 

3.3 

1.6 

Industrial 

4.6 

2.2 

Undeveloped  or  park  land 

1.5 

0.7 

1.  Loading  rates  in  units  of  pounds  per  day  per  100  feet  of  curb 


TABLE  IV-5 

POLLUTANT  LOADING  RATES 
SEATTLE  METROPOLITAN  AREA 


IV-10 


The  loading  rates  were  then  incorporated  into  the  Runoff  Block 
as  DDFACT.  Quality  simulation  runs  were  performed  for  the  calibration 
areas.  Lack  of  accurate  and  consistent  data  for  rainfall,  runoff,  and 
quality  for  the  few  select  storms  made  the  quality  calibration  extremely 
difficult  for  most  areas. 

The  amount  of  buildup  of  dust  and  dirt  is  dependent  on  the  number 
of  dry  days  preceding  the  rainfall  event  to  be  analyzed  and  the  frequency 
of  street  sweeping  activities.  For  the  quality  calibration  process, 
street  sweeping  schedule  information  for  three  of  the  calibration  areas 
was  obtained  from  City  of  Seattle,  Department  of  Engineering.  The 
frequency  of  street  sweeping  varied  from  areas  being  swept  once  a day 
to  areas  not  being  swept  on  a regular  basis.  The  number  of  dry  days 
preceding  a rainfall  event  was  determined  by  inspection  of  rainfall  data. 
The  guideline  presented  in  the  original  Stormwater  Management  Model  work 
which  stated  that  nearly  all  the  pollutant  is  washed  off  the  drainage 
basin  with  a rainfall  intensity  of  0.5  inches  per  hour  for  a duration 
of  one  hour  was  loosely  used  in  the  calibration  process  to  estimate  the 
number  of  dry  days  preceding  monitored  rainfall  events.  The  monitoring 
program  spanned  one  of  the  driest  springs  and  summers  on  record  and 
therefore  considerable  judgment  was  relied  upon  to  arrive  at  the  number 
of  dry  days  preceding  a storm.  A value  of  6 dry  days  was  used  for  the 
March  16,  1973,  storm. 

The  hydrographs  and  pollutographs  generated  from  the  Runoff 
model  are  used  as  input  to  the  Transport  model  and  Transport  Quality 
model.  The  values  for  the  runoff  quantity  and  quality  are  determined 
by  the  Runoff  model  simulation  and  the  other  two  models  have  only  a 
slight  effect  on  altering  either  except  for  the  attenuation  of  flow 
due  to  backwater  effects  which  can  be  modeled  by  the  Transport  model. 

At  the  present  time  it  is  not  possible  to  make  a quantitative 
statement  concerning  the  absolute  accuracy  of  the  quality  simulation 
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and  its  various  component  parameters.  In  many  cases  the  quality 
simulation  seems  reasonably  good  but  in  other  cases  the  results  may 
deviate  by  factors  of  two  or  more  from  observed  values.  It  is  probably 
safe  to  assume  that  the  models  are  reasonably  adequate  for  assessing 
the  quality  differences  between  various  drainage  alternati ves,  but  that 
more  data  needs  to  be  collected  on  the  rate  of  pollutant  buildup  on  and 
washoff  from  the  individual  watersheds  before  the  accuracy  of  surface 
quality  model  will  rival  that  of  surface  quantity  model. 
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CHAPTER  V 

ADDITIONAL  MODEL  DETAILS 


MATHEMATICAL  RELATIONSHIPS 

The  following  sections  give  the  mathematical  basis  of  the  three 
major  sections  of  the  urban  stormwater  model,  the  Runoff  Block,  the 
Transport  Block,  and  the  Transport  Quality  Block.  Details  of  the 
numerical  discretization,  underlying  assumptions  and  input  requirements 
are  also  described. 

RUNOFF  BLOCK 

A typical  subwatershed  as  represented  in  the  Runoff  Block  is 
shown  in  Figure  V-l . It  is  assumed  that  conditions  of  gradually  varied, 
uniform  flow  obtain  on  the  watershed,  and  that  flow  computations  can  be 
made  based  on  kinematic  wave  theory.  Three  flow  depths  are  shown  on 
the  figure: 

dQ  = depth  at  time  t, 

d^  = depth  at  time  t+At,  and 

d$  = maximum  depth  of  detention  storage. 

The  objective  of  the  calculations  which  pertain  to  this  element  is  to 
find  the  new  depth,  d^ , determining,  in  the  process,  the  outflow,  Q0, 
and  maintaining  mass  continuity  at  all  times.  To  accomplish  this,  two 
equations  must  be  solved  simultaneously.  First,  the  continuity  or 
storage  equation 

Ad  _ dl  - do  _ n t x Qi  " ^0 

At  At  K As 


V-l 


INFILTRATION: 


I =k,  + (k2-  kjJe'V 


FLOW:  , l/_  /j  . j .5, 


o0=-^s2w (4^!  _dJ 


STORAGE 


«■[»- 


Qj  -Q0 
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FIGURE  V-l  ELEMENTARY  FLOW  CALCULATIONS  TYPICAL  ELEMENT 
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where 


r 


R = rainfal 1 during  At, 

I = infiltration  to  groundwater  during  At, 

Q-j  = inflow  from  upstream  catchments, 

Q0  = outflow  from  catchment, 

As  = surface  area  of  element, 

and  second,  the  flow  equation,  which  is  a form  of  Manning's  equation 


where 


1.49  1/2 

n 


w 


5/3 


n = Manning  coefficient 
w = width  of  the  plane 


Here  we  have  two  equations  in  two  unknowns,  Q0  and  d-j  . Note  that  the 
flow  computation  is  based  on  the  average  depth  during  At,  and  that 
surface  detention  is  not  included  in  the  effective  depth  of  flow. 


(V-2) 


Rainfall  intensity  is  an  input  quantity,  variable  in  time,  but 
considered  constant  during  each  time  interval  At. 


Infiltration  is  computed  by  a modified  Horton  formula 

I = K1  + (K2  - K-| ) e_l<3t  ( V-3) 

where 

«i , Ko  = minimum  and  maximum  infiltration  rates, 
respectively,  and 

= exponential  rate  of  loss  in  infiltration  capacity. 

During  periods  of  light  or  zero  rainfall,  the  net  precipitation  value,  R-I, 
could  become  negative.  Traps  in  the  program  prevent  this  occurrence,  thus 
modifying  the  Horton  approach  somewhat. 
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The  simultaneous  solution  of  Equations  V-l  and  V-2  is  performed 
by  a Newton's  method,  which  easily  handles  the  nonlinearity  involved. 
First  the  equations  are  combined  and  rearranged  in  the  form 


F = o = Ad  - At  (kd5/3  + R ) 


where 


F = Newton's  function 


1.49  1/2 


w) 


/As, 


do  + dl 


d = - d$  - d0  - ds  + and 

Rnet  " (R  ' » * 

Then,  differentiating  yields 


dF  , 5 . n 2/3 

0Ta3T  " 1 - At  F k d 


(V-4) 


(V-5) 


The  numerical  method  is  a recursive  process  for  finding  the  value  of  Ad, 


(Ad) 


n + 1 


where  the  subscripts  refer  to  the  ntn  and  (n+l)^  iterations.  Repeated 
application  of  this  expression  converges  upon  F = 0,  which  is  the  desired 
results . 


( V-6) 


Catch  basins  collect  and  store  pollutants  during  dry  weather, 
flushing  rather  rapidly  in  the  early  hours  of  a storm.  Analysis  of 
their  behavior  is  mostly  empirical,  and  field  data  should  be  obtained 
whenever  possible.  The  program  computes  the  washout  of  material  from 
catch  basins  by  the  expression 

Mr 


Qc  t 


M_ 


1.6  V, 


( V-7 ) 


V-4 
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where 


M£  = mass  emission  rate  from  catch  basins  (mass/unit  time), 

Mc  = total  mass  in  catch  basins  at  beginning  of  time  interval, 

Qc  = flow  through  catch  basins  (runoff  from  pervious  and 
impervious  areas), 

Vc  = total  catch  basin  volume. 

The  total  available  mass,  Mc,  includes  that  remaining  from  the  original 
contents  at  the  start  of  the  storm  event,  plus  material  washed  into  the 
catch  basins  during  the  storm  by  che  surface  runoff. 

Surface  runoff  of  polluta.  fs  is  treated  much  the  same  as 
catch  basin  washout.  The  buildup  of  material  on  the  surface  is  based 
upon  data  contained  in  the  study  by  the  American  Public  Works  Association, 
January  1969,  entitled,  "Water  Pollution  Aspects  of  Urban  Runoff."  In 
this  study 

PQ  = f (K,  D,  G,  N)  (V-8) 

where 

Pg  = total  mass  on  watershed  at  start  of  storm, 

K = mass  of  constituent  per  unit  mass  of  dust  and 
dirt  (a  function  of  land  use) 

0 = rate  of  dust  and  dirt  accumulation  (also  a 

function  of  land  use) 

G = total  length  of  gutters  in  the  area, 

N = number  of  days  since  last  storm  or  street  cleaning. 

Auxiliary  computations  account  for  the  efficiency  of  the  street  cleaning 
process  and  other  minor  variations. 

Once  the  mass  of  the  pollutant,  P,  on  the  watershed  has  been 
determined  the  rate  of  washoff  to  the  local  drainage  course  (gutter)  can 


where 

k = runoff  coefficient 

P = the  mass  of  pollutant  on  the  subwatershed 

Now,  if  we  assume  that  the  runoff  coefficient  varies  directly  with  the 
rate  of  surface  runoff  we  can  write 

K = br  (V-10) 

where 

b = proportionality  coefficient 
r = rate  of  surface  runoff 


To  evaluate  b we  shall  assume  that  a uniform  runoff  of  0.5 
in/hour  would  wash  away  90  percent  of  the  pollutant  in  one  hour.  This 
yields  a value  of  b = 4.6  and  Equation  V-9  takes  the  form 

^ = 4.6  rt  (V-ll) 

Further,  on  the  watershed  itself  the  time  history  of  the  pollutant  is 
given  by  the  simple  exponential  decay  as 

P * Pq  e'4,6  rt  ( V - 1 2 ) 


To  route  a typical  pollutant-  through  the  surface  drainage  system 
we  can  write  a mass  balance  for  each  gutter  as 


dM 

3t 


i = l 


si 


(V-13) 


where 

M = the  total  mass  in  the  gutter  volume 

s^  = a flux  of  mass  in  or  out  of  the  control  volume 
n = the  number  of  mass  fluxes  associated  with 
the  gutter.  This  will  usually  be  upstream 
gutter  inflow  and  tributary  watershed  washoff 


* 
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If  we  further  define 


■I 


where 

V-  = the 

1 1 

volume  of 

¥ v 

the 

gutter 

C = the 

concentrati on 

of  pollutant 

and 

dM 

d(VC)  _ 

r 

dV  + dC 

ar  " 

dt 

U 

at  + V Ht 

or 

y.dC  = 

* at 

E si  - C 

dV- 

ar 

For  simplicity  let  dc/dt  = C and  d¥/dt  = V-  and  Equation  V-16  becomes 

n 

VC  = E s-  C -V 
i-1  1 

To  achieve  an  integration  of  equation  V-17  with  respect  to 
time  we  shall  make  the  following  assumption  concerning  the  behavior  of 
concentration  in  time. 

Ct+At  = Ct  + ~J  ^t  + ^t+At  ^ 

where 

Cf,  Ct+.t  = pollutant  concentration  at  a time,  t, 
and  at  a later  time  t+At 

C^.,  = the  rate  of  change  of  concentration  at 

times  t and  t+At,  respectively 

Solving  Equation  V- 13  for  Ct  we  obtain 


( V- 14) 


(V-15) 


(V-16) 


(V-17) 


(V-18) 


^t+At  = At  Ct+At  ' (At  Ct  + 

In  the  general  case  Ct  and  Ct  are  known  from  a previous  solution  or 
initial  conditions  and  can  be  treated  as  constants  for  any  time  step. 
Accordingly  we  shall  define 


(V-19) 


V-  7 


= — r + c 

At  Lt 


(V- 20) 


a = 


At 


( v-2i : 


which  substituted  in  Equation  V-19  gives 


Ct+At  a Ct+At  ' B 


(V- 22) 


As  Equations  V-17  and  V- 22  are  general  in  nature  and  refer  to 
conditions  at  a particular  instant  in  time,  they  can  be  combined  as  follows 


CV-  + ¥ (aC  - 0)  = E s. 

i = l 1 

or 

. n 

C¥  + Ca¥  - £V-  = E s. 

i = l 1 
n 

E S 

Usually  the  term  ^ i contains  several  mass  inflows,  and 
one  mass  outflow,  QC-,  from  the  gutter.  This  being  the  case  Equation  V- 23 
i s final ly  wri tten  as 

n 

CV-  + CaV-  - 3^-  = E s.  - QC 
i = l 1 

which  can  be  solved  for  the  concentration  at  the  end  of  the  time  step 


(V- 23) 


At  in  the  form 


C = 


n 

Z s.  + B¥ 
i = l 

V-  + aV-  + Q 


(V-24) 


(V-25) 


Equation  V-25  is  the  form  used  in  the  quality  routing  of  the 
Runoff  Block,  and  is  solved  sequentially  for  each  gutter  moving  down  the 
tree-shaped  structure  required  by  the  model.  The  outflow  from  each  gutter 
becomes  the  inflow  to  the  next,  and  total  mass  flux  is  summed  using  gutter, 
watershed  and  catch  basin  contributions. 
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TRANSPORT  SLOCK 

The  Transport  Block  uses  a link-node  description  of  a sewer 
system  which  lends  itself  very  nicely  to  the  discrete  representation  of 
the  physical  system  and  to  its  mathematical  solution.  The  conduit  system 
is  idealized  as  a series  of  links  (pipes)  which  are  connected  to  each 
other  at  nodes.  Links  and  nodes  each  have  very  precisely  defined  qualities 
which,  taken  together,  represent  the  entire  pipe  network.  Moreover,  the 
link-node  approach  is  very  useful  in  representing  the  various  flow 
control  devices. 

Links  transmit  flow  from  point  to  point.  Properties  associated 
with  the  links  are  roughness,  length,  cross-sectional  area,  hydraulic 
radius  and  surface  width.  The  last  three  properties  are  functions  of  the 
instantaneous  depth  of  flow.  Any  shape  cross-section  can  be  incorporated 
easily  into  the  computer  program;  at  present  it  will  accept  pipes  whicn 
are  rectangular,  circular,  horseshoe,  baskethandle  or  eggshape.  A 
trapezoidal  open  channel  is  also  permitted. 

The  primary  dependent  variable  in  the  links  is  the  discharge,  Q. 

It  is  assumed  that  Q is  constant  (in  space,  not  time)  in  the  link.  This 
is  a common  assumption  in  the  link-node  method.  It  should  be  noted  that 
the  velocity  and  the  cross-sectiona1  area  of  flow,  or  depth,  are  not 
constant  in  the  link. 

Modes  are  the  storage  elements  of  the  system.  The  variables 
associated  with  a node  are  volume,  head,  surface  area  and  inflow/outflow. 

The  primary  dependent  variable  is  the  head,  H,  which  is  assumed  to  be 

nstant  (in  space,  not  time),  except  for  some  special  cases  as  noted  below. 

A sketch  of  the  link-node  idealization  is  shown  in  Figure  V-2. 

A node  resembles  a manhole,  i.e.,  it  is  located  at  a point  and  has  a head 
associated  with  it.  Inflow  and  outflow  occur  at  a node.  However,  the 
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volume  of  the  node  is  roughly  equivalent  to  the  water  volume  in  the  half 
r ipe  lengths  which  enter  the  node.  The  link,  as  shown  in  the  figure,  is 
the  flow  path  transferring  fluid  from  one  node  to  the  next. 

This  conceptual  representation  of  a sewer  system,  which  is 
physically  as  well  as  mathematically  descriptive,  proves  to  be  very  useful. 
Dynamic  flow  in  sewers  is  an  extremely  complicated  process,  with  many 
subtleties  surrounding  the  actual  numerical  solution.  The  development  of 
the  model  showed  that  a thorough  understanding  of  both  tne  physical  and 
mathematical  aspects  is  essential  to  produce  a working  model. 

The  basic  differential  equations  for  the  sewer  flow  problem  are 
the  equations  of  motion  and  continuity.  The  general  motion  equation  for 
discharge  in  a conduit  not  flowing  full  is 


where 


+ 2u 


2 9A  « 

“ 57  -9-' 


5H 

3x 


Q = discharge, 
u = velocity, 

A = cross-sectional  area  of  the  flow, 
H = hydraulic  head, 

S^=  friction  slope, 
x = distance  along  conduit,  and 
t = time. 


The  friction  slope  can  be  expressed  as 


S 


f 


_k_ 

gAR/j 


QM 


( V-26) 
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where 


n = Manning's 
R = hydraulic 


The  absolute  value  sign 
the  flow. 


roughness  coefficient,  and 
radi us . 

ensures  that  the  frictional 


force  always  opposes 


There  are  a number  of  ways  to  express  the  motion  equation  in 
finite  form.  This  program  uses  the  expression 


u 

to  evaluate  the  derivative  of  Q at  any  time  t.  The  term  C ^ 1S  the 
entrance  and  exit  losses  in  the  conduit.  The  coefficients  C(  and  C2  are 
input  values.  The  subscripts  i,  j,  refer  to  the  nodal  values  at  the  ends 
of  the  link.  The  barred  symbols  for  R,  u and  A indicate  a weighted  average 
along  the  1 ink,  e.g. , 

*rr('i  ( 

where  the  subscript  m refers  to  the  midpoint  of  the  link. 


The  continuity  equation  for  a node  can  be  written  directly  in 


finite  form.  It  is,  for  time  t, 

[ 3H\  = LQt 

\9t't  Ast 


( V-30) 
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where 


£Qt  = algebraic  sum  of  all  conduit  flows  and 
I/O  flows  to  the  node, 

A ^ = water  surface  area  associated  with  the  node. 


When  all  of  the  pipes  entering  a node  are  flowing  full,  the 
node  is  said  to  be  surcharged.  The  usual  form  of  the  continuity  equation 
is  not  strictly  applicable  in  this  case,  because  no  additional  storage 
is  available  at  the  node.  Continuity  then  becomes: 


EQt  = 0 ( V- 31 ) 

and  the  earlier  form  is  indeterminate.  The  exact  change  in  hydraulic  head 
cannot  be  determined  directly;  the  situation  is  now  a dynamic  version  of  the 
traditional  pipe  network  problem.  However,  an  approximate,  first  order 
correction  based  on  the  Hardy  Cross  method  can  be  computed.  The  Hardy  Cross 
head  correction  for  a node  is 


-n£Q 

Ah  = £(a}  ( V- 32 ) 

where 

Ah  = change  in  nodal  head  for  one  iteration, 

Q = flow  in  a line  entering  the  node, 

H = pressure  drop  along  that  line,  and 
£ = summation  over  all  lines  entering  node. 

Successive  application  of  Equation  V- 32  to  the  nodes  of  the  pipe  network 
leads  to  the  satisfaction  of  Equation  V-31 . A derivation  similar  to  that 
used  to  obtain  Equation  V-32,  but  using  the  dynamic  equation  of  motion 
V- 28  yields 


g(At)2  £* 


( V-33) 
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where  K = Constant,  usually  taken  as  0.25  to  introduce  some  underrelaxation 
in  the  system.  Note  that  A is  the  cross-sectional  area  of  the  lines;  the 

term  K is  a measure  of  the  capacity  of  the  system  to  absorb  excess 

g(At)2zA 

flow.  It  is  a kind  of  equivalent  surface  area.  For  a surcharged  node, 
the  calculations  are  carried  out  using  Equation  V-33,  which  is  identical 
in  form  to  Equation  V-30.  The  process  could  be  an  iterative  one,  but 
for  small  time  steps,  where  AH  remains  small,  the  first  order  approximation 
has  been  found  to  work  sufficiently  well. 

Flooding  is  a special  case  of  surcharge,  when  the  hydraulic 
grade  line  breaks  the  ground  surface.  Traps  in  the  program  limit  the 
hydraulic  grade  line  to  a value  no  greater  than  the  ground  elevation. 

The  numerical  integration  of  equations  can  be  accomplished  by 
a modified  Euler  method.  The  results  are  accurate  and,  when  certain 
constraints  are  recognized,  stable.  Figure  V-3  shows  how  the  process 
would  work  if  only  one  equation  were  involved.  The  first  three  operations 
determine  the  slope  at  the  "half-step"  time.  Tnis_  is  used  in  operation  (4) 
to  project  the  "full -step".  In  other  words,  it  is  assumed  that  the  slope 
at  time  t + ^ is  the  mean  slope  during  the  interval.  The  method  is  extended 

easily  to  more  than  one  equation,  although  graphic  representation  is  then 
very  difficult. 

The  modified  Euler  method  yields  a completely  explicit  approach. 

The  motion  equation  is  applied  to  each  link  and  the  continuity  equation  to 
each  node,  entirely  without  implicit  coupling.  It  is  well  known  that 
explicit  methods  involve  fairly  simple  arithmetic  and  req  e little  storage 
space,  compared  to  implicit  methods.  However,  they  are  generally  less 
stable,  and  often  require  very  short  time  steps. 
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First  examination  of  the  results  seemed  to  validate  the  usual 
requirements  of  explicit  methods.  A closer  look  revealed  some  rather 
distinctive  features  which  are  worth  noting. 

The  solution  does  become  unstable  for  time  steps  which  exceed 
certain  criteria,  but  not  for  reasons  commonly  encountered.  The  typical 
stability  relationship  for  solutions  of  the  motion  equation  is  of  the 
form 

L 

S /gy  (V-34) 

that  is,  the  time  step  must  be  less  than  the  time  required  to  propagate 
a wave  through  the  reach.  Experience  showed  that  the  restraint  operating 
in  this  model  was  more  severe.  A more  appropriate  stability  relationship  is 


A+  C'As  Ymax 

at  - — m — 


tV-35) 


where 

C'  = constant  determined  by  experience,  usually  about  0.10. 

Physically,  the  cause  of  instability  is  too  rapid  a rise  in  the  water 
surface.  Equation  V-35,  with  C'  = 0.10,  says  that  the  rise  in  water 
surface  during  At  should  not  exceed  0.10  of  the  depth  available  in  the 
pipe.  This  criterion  clearly  is  approximate  and  must  be  approached  with 
caution. 


Application  of  the  finite  difference  form  of  the  motion  equation. 
Equation  V-28,  is  not  always  strai ghtforward  because  of  some  special 
conditions  which  may  exist  all  or  part  of  the  time  in  a sewer.  First, 
the  invert  elevations  of  pipes  which  join  at  a node  may  be  different; 
sewers  are  built  frequently  with  invert  discontinuities.  Second,  critical 
depth  may  occur  in  the  conduit.  Third,  normal  depth  may  control.  Finally, 
the  pipe  may  be  dry.  In  all  of  these  cases,  or  combinations  thereof,  the 
flow  must  be  computed  by  special  techniques.  Figure  V-4  shows  each  of  the 
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FIGURE  V-4  HYDRAULIC  CASES  HANDLED  BY  TRANSPORT  BLOCK 


possibilities  and  describes  the  way  in  which  surface  area  is  assigned 
to  the  nodes.  The  options  are: 


1.  Uormal  case.  Flow  computed  from  motion  equation; 
half  of  surface  area  assigned  to  each  node. 

2.  Critical  depth  downstream.  Use  lesser  of  critical 
or  normal  depth  downstream;  assign  all  surface  area 
to  upstream  node. 

3.  Critical  depth  upstream.  Use  critical  depth;  assign 
all  surface  area  to  downstream  node. 

4.  Flow  computed  exceeds  normal  flow.  Set  flow  to 
normal  value;  assign  surface  area  in  usual  manner. 

6.  Dry  pipe.  Set  flow  to  zero;  if  any  surface  area 
exists,  assign  to  downstream  node. 

Once  these  corrections  are  applied,  the  computations  can  proceed  in  the 
normal  manner.  Note  that  any  of  these  special  situations  may  begin  and 
end  at  various  times  and  places  during  simulation.  The  program  detects 
these  automatically. 

The  link-node  computations  can  be  extended  to  include  devices 
which  divert  sanitary  sewage  out  of  the  storm  drainage  system  or  relieve 
the  storm  load  on  sanitary  interceptors.  As  treated  here,  these 
diversions  are  assumed  to  take  place  at  a node,  and  are  handled  as  a 
special,  internodal  transfer. 

Weir  diversions  provide  relief  to  the  sanitary  system  during 
periods  of  storm  runoff.  Flow  over  a weir  is  computed  by  the  equation 


«„  ■ W<h  * k]i  ■ %>a] 


f 


where 


C = discharge  coefficient, 
w 

L = weir  length, 
w 

h = driving  head  on  the  weir, 
v = approach  velocity,  and 
a = weir  exponent;  3/2  for  transverse  weirs, 
s/3  for  si  deflow  weirs 


Both  C and  L are  input  values  for  transverse  weirs.  For  si  deflow  weirs, 
w w 

C should  be  a function  of  the  approach  velocity,  but  the  present  version 
of  the  program  does  not  include  this. 


Normally,  the  driving  head  on  the  weir  is  computed  to  be: 

h = y,  - yc 


( V- 37 ) 


where 

y = water  depth  on  upstream  side  of  weir, 
yc  = height  of  weir  crest  above  invert. 


However,  several  other  conditions  can  exist  which  modify  the  situation. 

If  the  downstream  water  depth,  y , also  exceeds  the  weir  crest  height,  the 
weir  is  submerged,  and  the  flow  computed  by  Equation  V-36  is  corrected 
according  to  Villemonte,  i.e.. 


-,0.385 


Qw  = Qw 

r i v*  - ycY  1 

L I*.  - v J 

( V-38) 


If  the  upstream  pipe  is  surcharged,  the  weir  will  behave  like  an  orifice, 
and  the  flow  is  computed  by  the  equation 


where 


Qw  = 


C'L  , 
w w 


(y 

VJmax 


• *c} 


+ V' 


Cv'  = discharge  coefficient  for  submerged  case,  and 

y = maximum  flow  depth. 

•'max  K 
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This  is  almost  equivalent  to  having  the  weir  behave  as  an  orifice.  The 

coefficient  C'  is  an  input  variable;  as  rule  of  thumb,  it  can  be  taken 
w 1 

approximately  equal  to  g Cw. 

For  all  of  the  above  cases,  the  direction  of  flow  may  be  reversed 
under  certain  conditions  which  cause  the  "downstream"  water  surface  to 
rise  above  the  "upstream".  This  does  not  affect  the  equations,  except  that 
yt  and  y2  are  switched.  However,  the  configuration  of  a sideflow  is  such 
that  reversed  flow  is  more  nearly  like  a transverse  weir;  the  weir  coefficient, 
a,  is  adjusted  to  f2  for  this  situation. 

Frequently,  weirs  are  installed  together  with  a tide  gate  at 
points  of  overflow  into  the  receiving  waters.  Flow  across  the  weir  is 
restricted  by  the  tide  gate,  which  may  be  partially  closed  at  times.  This 
is  accounted  for  by  reducing  the  effective  driving  head  across  the  weir 
according  to  an  empirical  factor  published  by  Armco  in  their  design  manual 
"Armco  Water  Control  Gates" 

-1 ,15v 

h'  = h - | v2e  ^ ( V-40) 

An  orifice  is  usually  installed  as  a dropout  to  divert  sanitary 
sewage  to  an  interceptor,  out  of  the  storm  system,  and  to  limit  the  flow  of 
stormwater  into  the  sanitary  system.  Flow  through  an  orifice  is  governed  by 

Qo  = C0A^gh  ( V-41 ) 

where 

Co  = discharge  coefficient, 

A = cross-section  area,  and 
h = driving  head  on  the  orifice. 

The  driving  head  is  usually  the  depth  of  flow  in  the  upstream  pipe,  yt. 

For  a few  instal lations , when  the  orifice  is  side  flowing  instead  of  dropout, 
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the  head  difference  Ht  - H2  may  control.  The  variables  y and  H are  defined 
identically  for  orifices  as  for  weirs,  and  backflow  can  occur. 

A pump  is  represented  conceptually  as  an  off-line  storage  node 
(the  wet  well)  from  which  the  contents  are  pumped  to  another  node  according 
to  a programmed  rule  curve.  One,  two  or  three  stage  pumping  is  permitted, 
and  the  turning  on  and  off  of  the  pumps  is  handled  automatically. 

The  program  simulates  weir  outfalls  and  free  outfalls,  both 
types  with  or  without  tide  gates.  The  characteristics  of  the  weir  outfall 
were  described  above.  The  "free"  outfall  may  be  truly  free,  if  the  elevation 
of  the  receiving  waters  is  low  enough,  or  it  may  consist  of  a backwater 
condition.  In  the  former  case,  the  water  surface  at  the  free  outfall  is 
taken  as  critical  or  normal  depth,  whichever  is  less.  If  backwater  exists, 
the  receiving  water  elevation  is  taken  as  the  value. 

If  a tide  gate  is  present,  the  Armco  correction  Equation  V-40, 
is  applied  to  either  the  weir  or  free  outfall  case.  The  outflow  is  zero, 
of  course,  when  the  tidal  elevation  exceeds  the  water  level  in  the  pipe, 
because  the  gate  is  closed. 

TRANSPORT  QUALITY  BLOCK 

Water  quality  constituents  are  routed  through  the  principal 
sewers  by  the  Transport  Quality  Model.  It  is  the  simplest  of  the  three 
models,  and  would  be  combined  with  the  Transport  Block  except  that 
computer  storage  requirements  will  not  permit  it. 

Most  of  the  movement  of  water  quality  constituents  in  the 
sewer  system  is  by  advective  transport.  The  differential  equation 
which  describes  this  process  is  seen  usually  in  one  of  two  forms: 
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In  terms  of  mass  rate  of  change. 


where 


3M 

ST  = 


QC 


M = mass,  and 
C = concentration. 


In  terms  of  concentration  rate  of  change, 

3C  _ C 3v  . „ 3C 
3t  " ' V aT  u a3T 

where 

V = volume. 


( V-42 ) 


(V-43) 


These  equations  are  completely  interchangeable  for  our 
purposes;  some  aspects  of  both  are  used  in  the  numerical  solution. 

The  conceptual  representation  of  advective  transport  is  shown 
in  Figure  V- 5.  Conduit  "n"  with  velocity  un,  connects  nodes  1 and  2, 
which  have  pollutant  concentrations  of  and  C^,  respectively,  at  time  t. 

During  time  interval  At,  the  initial  concentration  gradient  will  move  at 
velocity  un,  so  that  it  occupies  the  position  shown  by  the  dashed  line 
at  the  time  t+At.  The  value  C,  which  was  the  distance  upAt  away  from 
node  2,  becomes  the  new  value  at  node  2, 

This  part  of  the  advective  transport  process,  which  does  not 
account  for  volume  changes,  can  be  expressed  by  writing  only  the  last 
part  of  Equation  V-43  in  finite  form;  applicable  to  a conduit 

ACJ  _ ,.ACn  ( V-44) 

At  " ur 
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where 


L = length  of  conduit 

ACn  = change  in  concentration  along  the  conduit 
ACt  = change  in  nodal  concentration,  due  only  to 

J 

advection  in  a single  conduit. 


This  equation  is  applied  to  each  conduit  entering  a node,  provided 
that  the  flou  in  that  conduit  is  into  the  node.  This  implies  that 
the  concentration  at  the  node  is  blind  to  events  which  occur 
downstream.  It  is  an  important  concept  not  only  from  the  theoretical 
viewpoint  but  also  from  the  numerical  one,  for  it  greatly  reduces  the 
numerical  propagation  of  errors. 


Combining  the  conduit  computations  at  a node,  and  adding  the 
effects  of  volume  change,  source-sink  contributions,  pumps,  weirs  and 
orifices  produces  the  total  change  in  concentration  at  the  node 


where 


AC. 

_J 

At 


C.  AV. 
J _J 
VT  At 
J 


+ — E Qn  —A 
an  u n 


At 


Vj  = nodal  volume, 

= summation  over  all  upstream  conduits, 
an  = u^n 

Qn  = flow  in  conduit 

l = summation  of  all  source-sinks,  pumps,  etc. 


The  latter  term  in  Equation  V-45  must  be  computed  via  a mass  balance 
technique. 


(V 


The  time  history  of  concentration  at  all  nodes  can  be  obtained 
by  successively  integrating  Equation  V-45  forward  through  time.  Although 
many  integration  methods  could  be  used,  including  the  modified  Euler  method, 
the  equation  is  quite  well  behaved  and  it  is  sufficient  to  perform  the 
single  step  operation 


At 
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ACj  ( t+At)  = ACj(t)  + 


The  integration  step,  At,  must  be  restricted  so  that 


for  each  conduit,  or  an  unstable  solution  may  occur.  In  practice,  this 
restraint  can  be  violated  occasionally  without  causing  the  solution  to 
become  unstable,  although  some  error  must  be  introduced  thereby.  The 
program  contains  a trap  to  guard  against  the  occurrence. 


The  Transport  Block  routes  all  of  the  constituents  produced 
by  the  Runoff  Block.  It  performs  only  a mass  balance,  i.e.  products 
output  by  the  Runoff  Block  are  transferred  without  addition  or  deletion. 
Note  especially  that  BOD  does  not  decay  during  routing.  This  is  a good 
assumption,  because  the  total  transport  time  can  be  expected  to  be  quite 


A more  serious  assumption,  which  results  in  some  difficulty  in 
simulation,  is  that  sewer  systems  normally  contain  large  quantities  of 
stored  material  in  the  form  of  sludge  deposits,  grits,  etc.;  these  are 
ignored  by  this  program.  Consequently,  the  "plug  flow"  phenomenon  often 
observed  is  not  reproduced. 


PROGRAM  FLOWCHARTS 

The  following  figures  are  the  program  flowcharts  for  the 
subroutines  included  in  the  Runoff  Block,  the  Transport  Block,  and 
the  Transport  Quality  Block. 
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RUNOFF  BLOCK 


The  Runoff  Block  contains  the  following  figures 


Figure  V-6 

RUNOFF 

Figure  V-7 

GQUAL 

Figure  V-8 

GUTTER 

Figure  V-9 

HCURVE 

Figure  V-10 

HYDRO 

Figure  V-ll 

QSHED1 

Figure  V- 1 2 

QSHED2 

Figure  V-13 

RECAP 

Figure  V- 1 4 

RHYDRO 

Figure  V- 1 5 

WSHED 
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FIGURE  V- 


6 Subroutine  RUNOFF 
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FIGURE  V-7  (cont'd)  Subroutine  GQUAL 
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FIGURE  V- 9 SUBROUTINE  HCURVE 


V-31 


FIGURE  V-12  Subroutine  QSHED2 
V- 35 
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FIGURE  V-1 3 Subroutine  RECAP 


FIGURE  V- 1 5 Subroutine  WSHED 
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TRANSPORT  BLOCK 


The  Transport  Block  contains  the  following  figures 


Figure 

V- 1 6 

MAIN 

Figure 

V-l  7 

BOUND 

Figure 

V- 1 8 

DEPTH 

Figure 

V-l  9 

HEAD 

Fi gure 

V-20 

HYDRAD 

Fi gure 

V -21 

IN DATA 

Figure 

V -22 

INFLOW 

Figure 

V- 23 

OUTPUT 

Figure 

V- 24 

TIDCF 
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FIGURE  V - 1 6 (cont'd)  Subroutine  MAIN 

v-/n 


L 
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DOES  SYSTEM 
CONTAIN  PUMPS 


ARE  THERE 
FREE  OUTFALL 
NOOES 


ARE  THERE 
TIDAL  NOOES 


FIGURE  V-l 7 Subroutine  BOUND 
V-44 
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INTERPOLATE  AND  FIND 
CRITICAL  OEPTH 


FIGURE  V-18  Subroutine  DEPTH 
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FIGURE  V 


18  (cont'd)  Subroutine  DEPTH 
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CALL  HYRAD 


COMPUTE  VELOCITY, SURFACE 
AREA,  AND 
HYDRAULIC  RADIUS 


FIGURE  V-19  Subroutine  HEAD 
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FIGURE  V-20  Subroutine  HYDRAD 
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V-49 


ARE  PLOTS  REQUESTED 


FIGURE  V- 23  Subroutine  OUTPUT 
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TRANSPORT  QUALITY  BLOCK 

The  Transport  Quality  Block  contains  the  following  figures: 

MAIN 
IN DATA 
INPUT 
OUTPUT 


Figure  V- 25 
Figure  V- 26 
Figure  V- 27 
Figure  V-28 


CALI  INPUT 


REPEAT  FROM  A1  TO  B1 
FOR  ALL  TIME  STEPS 


CALL  HYORAD 


REPEAT  FROM  C?  TO  01 
FOR  ALL  CONDUITS 


• TRANSPORT  BLOCK 


CALL  OUTPUT 


COMPUTE  NOOAl 
VOLUMES 


FIGURE  V-25  Subroutine  MAIN 


FIGURE  V- 26  Subroutine  INDATA 
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FIGURE  V-27  Subroutine  INPUT 
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PROGRAM  VARIABLE  DESCRIPTIONS 


The  following  four  tables  define  variables  in  FORTRAN  storage. 

Table  V-l  defines  the  values  for  the  Display  Block,  Table  V- 2 defines 
the  values  for  the  Runoff  Block,  Table  V-3  gives  the  values  for  the  Transport 
Block,  and  Table  V-4  is  used  for  the  Transport  Quality  Block.  In  each  table 
the  COMMON  variables  are  defined  in  alphabetized  labeled  COMMON  order,  and 
then  the  remaining  variables  are  defined  for  each  alphabetized  subroutine. 
There  is  cross-referencing  in  each  table;  each  common  block  lists  the 
associated  subroutines,  and  each  subroutine  lists  the  associated  common 
blocks. 
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TABLE  V-l 

DISPLAY  BLOCK  VARIABLES 
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Variable 

Name  Description  Unit 

SUBROUTINE  CURVE 

A The  log  base  10  of  the  range  of  values  of  y 

coordinate  to  be  plotted 

FRANG  Expanded  range  (even  intervals)  of  y 

coordinates  of  curve  to  be  plotted 

K Subscript  counter 

L Subscript  counter 

M Subscript  counter 

N Subscript  counter 

NCV  Number  of  curves/plot 

NPLOT  Number  of  plots 

NPOINT  Number  of  points  on  a plot 

NPT  Number  of  points/curve  (array) 

NPTM  Numerical  value  of  NPT 

RANGE  Range  of  y values  to  be  plotted 

X X coordinate  array 

X I NT  Label  interval  for  X 

XMAX  Maximum  X value 

XMIN  Minimum  X value 

X0  Start  point  of  line  (X  coordinate) 

XSCAL  X scale  factor 

XT  End  point  of  line  (X  coordinate) 
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TABLE  v-l 

(Continued)  2 of  5 


Variable 

Name 

Description 

Unit 

V 

Y coordinates  of  curves  to  be  drawn 

YINT 

Label  interval  for  Y 

YMAX 

Maximum  Y value 

YMIN 

Minimum  Y value 

YO 

Start  point  of  line  (Y  coordinate) 

YSCAL 

Y scale  factor 

YT 

End  of  line  (Y  coordinate) 

This  subroutine  uses  the  following  common: 

COMMON/LAB/ 

SUBROUTINE  GRAPH 

IC  Calling  sequence  control  parameter 

II  Subscript  counter 

ILAB  Output  label  with  plot 

IPLOT  Array  of  nodes  to  be  plotted 

ITAB  Array  indicating  which  locations  of  the  data 

file  are  to  be  plotted 

J Subscript  counter 

JJJ  Subscript  counter 

K Subscript  counter 

L Subscript  counter 

LX  Transfer  location  from  data  file  to  plot  storage 
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TABLE  V- 1 
(Continued) 
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Variable 

Name 

M 

MC 

MM 

N 

NCURVE 

NCV 

NLOC 

NLP 

NN 

NPCV 

NPLOT 

NPT 

NQP 

NQUAL 

NR 

NSTEPS 

NTAPE 

NVAL 

TAREA 

TDELT 


Description 

Subscript  counter 

Do  loop  counter 

Subscript  counter 

Subscript  counter 

Number  of  curves  to  be  plotted 

Number  of  curves/plot 

Node  number  of  hydrograph  point 

Number  of  types  of  plot  (hydrographs  and 
pol 1 utographs) 

Subscript  counter 

Maximum  number  of  curves/plot 

Number  of  plots 

Array  containing  number  of  points  to  be 
plotted  (GRAPH) 

Number  of  quality  constituents  to  be  plotted 

Number  of  quality  constituents  on  data  file 

Subscript  counter 

Number  of  steps  in  plot 

Input  tape  number  for  plotting 

Number  of  points/data  record  on  a file 

Total  area 

Time- step  interval 
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TABLE  V-l 
(Conti nued) 


4 of  5 


Variable 

Name  Description 

XA  X increment  used  for  interpolation 

XI  Same  as  XO 

X2  Same  as  XT 

YA  Y increment  used  for  interpolation 

Y1  Same  as  Yo 

Y2  Same  as  YT 

This  subroutine  does  not  use  any  common. 

SUBROUTINE  PPLOT 

A Transfer  array  for  plotting 

I Subscript  counter 

II  Subscript  counter 

IX  Start  point  of  line 

IY  Start  point  of  line 

J Subscript  counter 

JJ  Subscript  counter 

NCT  Number  of  plots 

SYM  Plot  symbol  array 

This  subroutine  uses  the  following  common: 
COMMON/LAB/ 
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TABLE  V-I 
(Continued) 

5 of  5 

Variabl  e 
Name 

Description 

Unit 

TIMES 

Time-step  interval 

TZERO 

Zero  time 

/ . 

X coordinate  array  (GRAPH) 

Y 

Y coordinates  of  curves  to  be  drawn 

YT 

Hydrograph-pol 1 utograph  information  on 
data  file 

This  subroutine  uses  the  following  common: 
COMMON/TAPES/ 

COMMON/ LAB 

It  also  has  its  own  blank  common  which 
overlays  main  program  blank  common. 

SUBROUTINE  PINE 

AXA 

X coordinate  of  value  previously  plotted 

AXG 

X coordinate  of  value  to  be  plotted 

AYA 

Y coordinate  of  value  previously  plotted 

AYB 

Y coordinate  of  value  to  be  plotted 

IXA 

Integer  value  of  AXA 

IXB 

Integer  value  of  AXB 

IYA 

Integer  value  of  AYA 

IYB 

Integer  value  of  AYB 

N 

Subscript  counter 

NCT 

Number  of  plots 

NSYfl 

Plot  number 
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Vari able 
Name 


DEPIN 
RA ININ 


HORIZ 

IT 

TITLE 

VERT 

XLAB 

YLAB 


NAMEG 

ilGTO 


Description 


Unit 


COMMON/ INFIL/ 

Maximum  allowable  infiltration  inches 

Summation  of  infiltration  inches 

This  common  is  used  in  the  following  subroutines: 

HYDRO 

RHYDRO 

WSHED 


COM' ION/ LAB/ 


Curve  label  for  x-axis 

none 

Internal  control  variable  subroutine  HCURVE 

( PPLOT ) 

none 

Ti tie  printed  out  wi th 

graphs 

none 

Curve  label  for  y-axis 

none 

Numerical  scale  labels 

for  x-axis 

none 

Numerical  scale  labels 

for  y-axis 

none 

This  common  is  used  in 

the  following  subroutines: 

CURVE 

GRAPH 

HCURVE 

PPLOT 

C0M1  ION/ NEW/ 

Gutter  member 

none 

Gutter  member  to  which 

watershed  drains 

none 

This  common  is  used  in  the  following  subroutines: 

GUTTER 

RHYDRO 
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Variable 

Name 


BASINS 

C 

CBFACT 

CBVOL 

CDOT 

CLFREQ 

DDFACT 

DRYDAY 

GQLEN 

KLAND 

NPASS 

NQS 

PBASIN 

POFF 

PSHED 

QFACT 


Descri pti on 


Unit 


COMMON/ ABLK/ 


Equivalent  number  of  standard  catchbasins 

Removal  coefficient 

Catchbasin  factor 

Volume  of  liquid  in  catchbasin 

Internal  variable,  change  in  concentration 

Interval  between  street  cleanings 

Dust  and  dirt  factor 

Number  of  dry  days  prior  to  storm 

Number  of  feet  of  gutter  in  subcatchment 

Land  use 

Number  of  streetsweeper  passes  per  cleaning 

Number  of  quality  constituents 

Mass  in  catchbasins 

Rate  of  mass  runoff 

Mass  on  subcatchment 

Flow  factor 

This  common  is  used  by  the  following  subroutines 


none 

none 

none 

CF 

days 

lb ‘/day/ 

100  ft.  gutter 

days 

100  ft 

none 

none 

none 

grams 

gr/sec. 

grams 

mg/ gram 


BLOCK  DATA 

GQUAL 

GUTTER 

HYDRO 

QSHED1 

QSHED2 

RHYDRO 
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Variable 

Name 


INCNT 

IOUTCT 

JIN 

JOUT 

NSCRAT 


DECAY 

DELD 

DELT 

DELT2 

DFULL 


Description 


Uni  t 


COMMON/TAPES/ 

Not  used  at  present 
Not  used  at  present 

Unit  member  of  magnetic  tape  or  drum/disc  none 

used  for  non-card  input 

Unit  member  of  magnetic  tape  or  drum/disc  none 

used  for  non-card  output 

Number  of  magnetic  tape  or  drum/disc  used  none 

internally 

This  common  is  used  in  the  following  subroutines: 

GQUAL 

GRAPH 

GUTTER 

HYDRO 

QSHED1 

RHYDRO 

RECAP 

RUNOFF 


COMMON  (UNLABELED) 

Exponential 

decay  rate  for  infiltration 

1 /sec 

Instantaneous  pipe  diameter  in  radians 

radi an 

Integration 

time  ’nterval 

sec,  min 

One  half  of 

a time  step 

sec,  min 

Gutter's  maximum  depth 

inches 
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— 

Variable 

Name  Description 


Uni  t 


DUMMY* 

Dummy  Variable 

none 

FLOW 

Hydrograph  flow  value 

cfs 

FLOW 

Temporary  variable  for  printing  flow 

cfs 

GCON 

Manning's  equation  less  hydraulic  radius 

none 

3DEPTH 

Instantaneous  gutter  depth 

inches 

GFLOW 

Gutter  flow 

cfs 

GLEN 

Length  of  gutter/pipe 

feet 

Gfl 

Manning's  roughness  coefficient 

none 

GSLOPE 

Slope  of  gutter/pipe 

ft/ft 

GS1 , 

Gutter  side  slope,  left  and  right 

ft/ft 

GWIDTH 

Pipe  diameter  or  gutter  width 

ft 

HGRAPH* 

Magnitude  of  variable  to  be  printed  in  vertical 
coordinate  of  the  curve 

none 

HISTOG 

Length  of  histogram  expressed  in  time 

sec 

HTIME* 

Time  interval  to  be  printed  in  the  horizontal 
coordinate  of  the  curve 

none 

IpLOT 

lumbers  of  gutters  and  inlets  to  be  plotted 

none 

INTCNT 

Printing  counter 

none 

INTERV 

Interval  integrization  cycles  for  printed 
hydrographs 

none 

IPRNT 

Points  for  which  hydrographs  will  be  printed 

none 

I SAVE 

Points  for  which  hvdrographs  will  be  saved 

none 

NAMEW 

External  subcatchment  number 

none 

ICODE 

Plot  control  integer,  zero  means  no  plot, 
one  means  plot 

none 
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TABLE  V-2 


(Continued) 
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Variable 

Name 

Descri pti on 

Unit 

NG 

Number  of  gutters 

none 

NGTOG 

Gutter  connections 

none 

NGTOI 

Inlet  connections 

none 

N6UT 

Bookkeeping  integer 

none 

NHISTO 

Number  of  rainfall  time  intervals 

none 

NHR 

Hour 

hr 

NHYET 

Number  of  hyetograph 

none 

NIN 

Maximum  number  of  gutters  draining  to  gutter, 
and  watersheds  draining  to  gutter 

none 

NING 

Do  loop  counter 

none 

NOG 

Total  number  of  gutters/pipes 

none 

NOW 

Total  number  of  watersheds 

none 

NPG 

Control  switch  for  type  of  gutter 

1.  regular 

2.  pipe 

3.  dummy  connected  directly  to  inlet 

none 

NPLOT 

Number  of  inlets  and  gutters  to  be  plotted 

none 

NPRNT 

Number  of  time  steps  between  printing 

none 

NRAIN 

Number  of  rainfall 

none 

NRGAG 

Number  of  hyetographs 

none 

NSAVE 

Number  of  points  where  hydrographs  are  saved 

none 

NSTEP 

Number  of  time  steps 

none 

NU 

Number  of  watershed 

none 

NWTOG 

Gutter  connection 

none 
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Variable 

Name 

Description 

Unit 

NWTOI 

Inlet  connection 

none 

OUT 

Temporary  variable  for  printing  concentration 

mg/1 

OUTFLW 

Flow  out  of  the  gutter 

cfs 

PCIMP 

Percent  imperviousness  of  watershed 

Of 

h 

PCTZER 

Percent  of  impervious  area  with  zero  detention 
depth 

QUAL 

Concentration  of  quality  constituents 

mg/1 

qin 

Input  from  upstream  gutter 

cfs 

QSUR 

Surcharge 

cf 

RAIN 

Rainfall 

in/hr  and 
ft/sec 

RI 

Instantaneous  rainfall  rate 

ft/sec 

RLOSS 

Infiltation  loss,  instantaneous 

ft/sec 

SUMI 

Total  infiltration  into  ground 

cf 

SUMOFF 

Total  gutter  flow  @ inlet 

cf 

SUMQW 

Total  flow  fc.r  each  subcatchment 

cf 

SUMR 

Total  rainfall 

cf 

SUMST 

Total  surface  storage 

cf 

TAREA 

Total  area 

acres 

THISTO 

Time  of  rainfall  time  intervals 

min 

TIME 

Time 

sec 

TIME2 

Time  minus  hal f step 

sec 

TITLE** 

Description  of  problem 

none 
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TABLE  V-2 

(Continued)  7 of  17 

Variable 


Name 

Description 

Uni  t 

TMIN 

Minute 

min 

TRAIN 

Time  when  rainfall  ends 

min/sec 

TZERO 

Starting  time  of  the  hydrograph 

sec 

WAREA 

Area  of  watershed 

acres/square  foot 

WCON 

Modified  Manning's  equations,  impervious 
and  pervious  portions  of  tne  watershed 

none 

WDEPTH 

Instantaneous  depth  on  watershed 

ft 

WFLOW 

Instantaneous  flow  on  watershed 

cfs 

WLMAX 

Maximum  infiltration  rate 

in/hr 

WLMIN 

Minimum  infiltration  rates 

in/nr 

WN 

Dunriy  variable 

none 

WSLOPE 

Average  slope  of  watershed 

ft/ft 

WSTORE 

Minimum  and  maximum  storage  depth  on 
surface  of  watershed 

ft 

WUIDTH 

Average  width  of  watershed 

ft 

*USED  ONLY 

IN  HCURVE  AND  HYDRO 

**USED  ONLY  IN  QSHED1 , HYDRO,  RECAP  AND  RHYDRO 

This  common  is  used  in  tne  following  subroutines: 

GQUAL 

GUTTER 

HCURVE 

HYDRO 

QSHED1 

RHYDRO 

WSHED 
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Variable 

Name  Description 


DCORR 

DEL 

DELR 

DF 

00 

EXPON 

F 

1 

I [ID 

J 

K 

NGAG 

WAR 

UFLO 


ALFA 

BETA 

K 


SUBROUTINE  WSHED 
Time-step  water  depth 

Time-step  change  in  depth  of  watershed  flow 

Newton-Raphson  change  in  depth  for  correction 

Sum  of  volume  change  plus  flow  change  times  ti 

Instantaneous  depth 

Decay  times  time  interval 

Bookkeeping  integer 

Bookkeeping  integer 

Rain  time  period 

Bookkeeping  integer 

Bookkeeping  integer 

Hyetograph  number 

Impervious  area  of  watershed  with  immediate 
runoff 

Average  watershed  flow  during  time  interval 

This  subroutine  uses  the  following  common: 

COMMON/ INF  I L/ 

Un labeled  COMMON 

SUBROUTINE  GQUAL 

Half  of  integration  time  step 

Interval  variable 

Bookkeeping  integer 


Unit 

ft 

ft 


sq  ft 
cfs 

1 /sec 
m9 

1 -sec 


none 
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Variable 

Name 

Description 

Unit 

L 

Bookkeeping  integer 

none 

M 

Bookkeeping  integer 

none 

TEMP 

Internal  variable 

sec 
cu  ft 

V 

Storage  volume 

cu  ft 

VDOT 

Change  in  storage  rate 

This  subroutine  uses  the  following  common: 

/ ABL  K/ 

/TAPES/ 

UNLABELED  COMMON 

SUBROUTINE  GUTTER 

cfs 

AX0 

Trapezoidal  cross-sectional  area,  starting 

sq  ft 

AX1 

D 

Trapezoidal  cross-sectional  area,  final 
Computational  variable,  internal 

sq  ft 

DAX1 

DDELV 

DEL 

Change  in  trapezoidal  cross-sectional  area 
Rate  of  change  in  volume  change 
Time-step  change  in  depth  of  watershed  flow 

sq  ft 

DELV 

Average  volume  change 

cu  ft 

DF 

Sum  of  volume  change  plus  flow 
change  times  time 

cu  ft 

DFLOW1 

Change  in  flow 

cfs 

DWP1 

Change  in  wetted  perimeter 

ft 

D0 

Instantaneous  depth 

ft 

D1 

Estimated  final  depth 

ft 
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Variable 

Name  Description 


Uni  t 


F 

FLOW 

FLOW0 

FL0W1 

I 

I F LG 
INLETS 

J 

K 

MOD 

M 

NOGG 

NOUT 

NQT 

NTIMEH 

NTSCP 

riTsi 

MTS2 

NT  1 

NX 

RAD0 

RADI 


Newton- Raphson  volume  correction 

Average  flow 

Starting  flow 

Final  flow 

Bookkeeping  integer 

Surcharge  indicator 

Same  as  NSAVE,  number  of  points  where  hydrographs 
are  saved 

Bookkeeping  integer 
Bookkeeping  integer 
Library  function 
Bookkeeping  integer 

Sum  of  number  of  gutters  plus  number  of  points 
where  hydrographs  are  saved 

Output  file  variable 

Same  as  NQS, 

Hour  of  day  of  simulation  (24-hour  clock) 

Number  of  times  surcharged 
Output  file  variable 
Output  file  variable 
Output  file  variable 
Bookkeeping  integer 
Starting  hydraulic  radius 
Final  hydraulic  radius 


cfs 

cfs 

cfs 


hr 


ft 

ft 
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Variable 

Name  Description 


Unit 


BASIN 

ERROR 

I 


IE0F 

II 

IJ 

ISUB 

J 


JK 

KSPOT 

M 


N 

NRANVL 

NSPOT 

NT1 

NX 


This  subroutine  uses  the  following  common: 

COMMON/ LAB/ 

Unlabeled  Common 

SUBROUTINE  HYDRO 

Basin  number 
Name  of  error  statement 
Bookkeeping  integer 
Indicator  integer 
Bookkeeping  integer 
Bookkeeping  integer 
Bookkeeping  integer 
Bookkeeping  integer 
Bookkeeping  integer 
Bookkeeping  integer 
Bookkeeping  integer 
Bookkeeping  integer 
Rain  data  points  limiter 
Bookkeeping  integer 
Output  file  variable 
Bookkeeping  integer 

This  subroutine  uses  the  following  common: 
COMMON/ ABLK/ 
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Variable 

Name  Description 


TMX  Time  in  hours 

WP0  Wetted  perimeter,  starting 

WP1  Wetted  perimeter,  final 

This  subroutine  uses  the  following  common: 

Common / ABLK/ 

Common/NEW/ 

Common/TAPES/ 

Un labeled  Common 

SUBROUTINE  HCURVE 
I Bookkeeping  integer 

J Bookkeeping  integer 

JT  Bookkeeping  integer 

K Bookkeeping  integer 

L Bookkeeping  integer 

M Bookkeeping  integer 

N Bookkeeping  integer 

NGAGP  Number  of  graphic  points 

ORIZ  Horizontal  title  unit  for  hydrograph  in  time 

TITEL  Description  of  curve  in  horizontal  coordinates 

TITL  Description  of  curve  in  vertical  coordinates 

TMAX  Maximum  time  to  be  printed  on  curve 

VER  Vertical  title  unit  for  hydrograph 
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Unit 


hr 

ft 

ft 


hr 

hr 
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Variable 

Name  Description  Unit 

COMMON/ INF  I L/ 

COMMON/TAPES/ 

Un labeled  COMMON 

SUBROUTINE  PRESET 


This  is  a block  data  subroutine.  It  uses  the 


DD 

following  common: 

COMMON/ ABLK 

SUBROUTINE  QSHED1  & QSHED2 
Mass  of  dust  and  dirt  on  watershed 

grams 

DFACT 

Decay  term  for  mass  on  watershed 

J 

Bookkeeping  integer 

none 

K 

Bookkeeping  integer 

none 

N 

Bookkeeping  integer 

none 

NCLEAN 

Number  of  cleanings  during  dry  period 

j 

none 

R 

Flow  from  watershed 

in/sec 

RATE 

Runoff  mass  rate 

\ 

1 /sec 

REFF 

Removal  efficiency  in  decimal  percent 

none 

TGS 

Parameter  to  account  for  mass  accumulation 

none 

in  spite  of  cleaning 

This  subroutine  uses  the  following  common: 

/ABLK/ 

/TAPES/ 

Unlabeled  COMMON 
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Variable 

Name  Description  Unit 


SUBROUTINE  RECAP 


DELT 

Integration  time  interval 

sec 

J 

Bookkeeping  integer 

none 

K 

Bookkeeping  integer 

none 

M 

Bookkeeping  integer 

none 

MAX 

Bookkeeping  integer 

none 

M 

11 

Bookkeeping  integer 

none 

NLT 

Bookkeeping  integer 

none 

NPTS 

Number  of  points  to  be  plotted 

none 

nqs 

Number  of  quality  constituents 

none 

NRD 

Bookkeeping  integer 

none 

NREAD 

Same  as  NSTEP 

none 

NSTEP 

Number  of  time  steps 

none 

NTIMEH 

Hour  of  day  of  simulation  (24-hour  clock) 

hr 

NTX 

NT1 

TAREA 

Total  area 

acres 

TIMEM 

Time  of  simulation  (24-hour  clock) 

mi  n 

TZERO 

Start  time  of  hydrograph 

sec 

This  subroutine  uses  the  following  common: 

COMMON/TAPES/ 
lln labeled  COMMON 

V- 77 


TABLE  V-2 
(Continued) 


15  of  17 


Variable 

.Tame  Description 


GA 


GP 

GQ 

GR 

G V 
G1 

G2 

G3 

G4 

G5 

G6 

G7 

I 

INLETS 

IN.UM 

J 

JG 

JK 

JW 

K 

KL 


SUBROUTINE  RHYDRO 

Internal  variable,  cross-sectional  area  of 
gutter  when  full 

Internal  variable,  method  perimeter  of  gutter 
when  full 

Internal  variable,  flow  in  gutter  when  full 

Internal  variable,  hydraulic  radius  of  gutter 
when  full 

Internal  variable,  velocity  in  gutter  when  full 

Read  in  value  of  bottom  width  of  gutter  or 
pipe  diameter 

Read  in  value  of  length  of  gutter 
Read  in  value  of  invert  slope 
Read  in  value  of  left-land  side  slope 
Read  in  value  of  right-hand  side  slope 
Read  in  value  of  Manning's  coefficient 
Read  in  value  of  depth  of  gutter  when  full 
Bookkeeping  integer 

Integer  counter,  number  of  dummy  gutters 

Bookkeeping  integer 

Bookkeeping  integer 

Bookkeeping  integer 

Number  of  hyetograph 

Bookkeeping  integer 

Bookkeeping  integer 

Do  loop  counter 

V-  78 


Uni  t 


ft 

ft 

ft/ft 

ft/ft 

ft/ft 

in 
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Variable 


Name 

Description 

Uni  t 

SUBROUTINE  RECAP 

DELT 

Integration  time  interval 

sec 

J 

Bookkeeping  integer 

none 

K 

Bookkeeping  integer 

none 

M 

Bookkeeping  integer 

none 

MAX 

Bookkeeping  integer 

none 

N 

Bookkeeping  integer 

none 

NLT 

Bookkeeping  integer 

none 

NPTS 

Number  of  points  to  be  plotted 

none 

NQS 

Number  of  quality  constituents 

none 

NRD 

Bookkeeping  integer 

none 

NREAD 

Same  as  NSTEP 

none 

NSTEP 

Number  of  time  steps 

none 

NTIMEH 

Hour  of  day  of  simulation  (24-hour  clock) 

hr 

NTX 

NT1 

TAPEA 

Total  area 

acres 

TIMEM 

Time  of  simulation  (24-hour  clock) 

mi  n 

TZERO 

Start  time  of  hydrograph 

sec 

This  subroutine  uses  the  following  common: 

COMMON/TAPES/ 

(In labeled  COMMON 
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Variable 

.lame 


INLETS 


TABLE  V-2 
( Continued) 


Descri pti on 


SUBROUTINE  RHYDRO 

Internal  variable,  cross-sectional  area  of 
gutter  when  ful 1 

Internal  variable,  method  perimeter  of  gutter 
when  ful 1 

Internal  variable,  flow  in  gutter  when  full 

Internal  variable,  hydraulic  radius  of  gutter 
when  full 

Internal  variable,  velocity  in  gutter  when  full 

Read  in  value  of  bottom  width  of  gutter  or 
pipe  diameter 

Read  in  value  of  length  of  gutter 
Read  in  value  of  invert  slope 
Read  in  value  of  left-land  side  slope 
Read  in  value  of  right-hand  side  slope 
Read  in  value  of  Manning's  coefficient 
Read  in  value  of  depth  of  gutter  when  full 
Bookkeeping  integer 

Integer  counter,  number  of  dummy  gutters 

Bookkeeping  integer 

Bookkeeping  integer 

Bookkeeping  integer 

Number  of  hyetograph 

Bookkeeping  integer 

Bookkeeping  integer 


Do  loop  counter 
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Variable 

Name  Description 


Unit 


MOD 

N 

NGOTO 

NHR 

NMN 

NN 

NP 

NPTS 

NQT 

NSTP 

NTX 

W1 

W2 

W3 

W4 

W5 

W6 

W7 

W8 

W9 

WIO 


Library  function 
Bookkeeping  integer 

Gutter  number  to  which  watershed  drains 

Hour  of  start  of  storm  hr 

Minutes  of  start  of  storm  min 

Bookkeeping  integer 

Read  in  value  of  NPG 

Same  as  INLETS  or  NPRNT 

Same  as  NQS  (COMMON/ABLK/ ) 

Bookkeeping  integer 
Tape  unit  identifier 

Read  in  value  of  the  average  width  of  watershed  ft 

Read  in  value  of  the  area  of  watershed  acre 

Read  in  value  of  the  percent  of  imperviousness 


Read  in  value  of  the  slope  of  watershed  ft/ft 

Resistance  factor  for  impervious  area 

Resistance  factor  for  pervious  area 

Retention  storage  depth  for  impervious  area  in 

Retention  storage  depth  for  pervious  area  in 

Read  in  value  of  maximum  infiltration  rate  in/hr 

Read  in  value  of  minimum  infiltration  rate  in/hr 
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(Continued) 

Variable 

Name  Description 

Wll  Read  in  value  of  decay  rate  of  infiltration 

W12  Conversion  factor,  default  value  of  12 

This  subroutine  uses  the  following  common: 

COMMON/ ABLK/ 

COMMON/ INF IL/ 

COMMON/NEW/ 

COMMON/TAPES 
Un labeled  COMMON 

SUBROUTINE  RUNOFF 

J Bookkeeping  integer 

K Bookkeeping  integer 

TFIX  Tape  unit  number 

This  subroutine  uses  the  following  common: 

COMMON/TAPES/ 

Unlabeled  COMMON 
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Unit 


1/sec 


none 

none 

none 
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TABLE  V- 3 

TRANSPORT  BLOCK  VARIABLES  1 of  23 


V a r i a b 1 e 
Name 

Descr intion 

Uni  t 

COMMON /BD/ 

ANORM 

Matrix  of  normalized  wet  cross  sectional 
area  of  conduit,  based  on  shape  and  depth 

HRNORM 

Matrix  of  normalized  hydraulic  radius  of 
conduit,  based  on  shape  and  depth 

TWNORM 

Matrix  of  normalized  conduit  width  at  flow 
line,  based  on  shape  and  depth 

This  common  block  is  used  by  the  following 
subroutines 
INDATA 
DEPTH 
HYDRAD 

COMMON/ BN D/ 

JFREE 

Node  for  free  outfall 

None 

JAATE 

Node  for  non-weir  tide  gate 

None 

JTIDE 

Not  used  at  this  time 

None 

".'FREE 

Number  of  free  outfalls 

None 

NGATE 

Number  of  non-weir  tide  gates 

None 

NUDE 

Indicator  for  outfall  tide  level  control 

None 

1.  no  water  surface  at  outfall 

2.  outfall  control  water  surface  at 
constant  elevation,  A1 

3.  tide  coefficients  provided 

4.  nronran  will  compute  tide  coefficients 
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Variable 

Name 


ALPHA 
A1  ' 
A2 
A3 

M > 

A5 

\ g 

A7  „ 

DELTO 

d:lt 

DCI.T2 

ICYC 

MJSW 

NC 

NJ 

HJSW 


Description  Unit 

This  common  block  is  used  by  the  fol lowinq 
subroutines 
HOUND 
IN DATA 

C011M0N/C0NTR/ 

Title  for  print  inn  None 

Feet 

None 

None 

Coefficients  of  the  seven  term  equation  M 


for  tidal  input 

None 

None 

None 

Not  used  at  this  time 

Time  interval  of  inteoration  Seconds 

One-half  of  BELT  Seconds 

Internal  cycle  counter  None 

Number  of  tape  input  hydrographs  None 

Number  of  conduits  None 

Number  of  nodes  None 

Number  of  nodes  for  hydrogranh  input  None 

via  cards 
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Variable 

Name  Description 


Unit 


NTC 

NTCYC 

N5 

N6 

1121 


TIME 

TIME2 

TZERO 

W 


Number  of  nodal  links  includina  internal 
1 inks 

, 'lumber  of  integration  cycles 

Input  unit  number 

Output  unit  number 

Unit  number  for  input  tape  from 
RUNOFF  BLOCK 

Time  counter  for  hydrograph  input 
TIME  - DELT2 

Zero  time  for  the  analysis 

Fundamental  freguency  of  daily  tidal 
variation 


None 

None 

None 

None 

None 

Seconds 
Seconds 
Hours/Seconds 
Rad /Seconds 


This  common  block  is  used  by  the  following 
subrouti nes 
BOUND 
IN DATA 
INFLOW 
MAIN 
OUTPUT 
TIDCF 


COMMON /F.  LEV/ 


IPLT 

Plot  control  integer 

None 

ZCRN 

Plot  variable,  highest  crown  elevation 

Feet 

at  a node 

ZGRNO 

Plot  variable,  ground  elevation 

Feet 

ZINVRT 

Plot  variable,  node  invert  elevation 

Feet 
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Variable 

Name  Description 


Unit 


This  common  block  is  used  by  the  following 
subroutines 
OUTPUT 


common/hyflow/ 

ISW 

Hydrograph  input  node  number  from 

tape 

None 

JSW 

Hydrooraph  input  node  number  from 

cards 

None 

NINREC 

Counter  for  hydrograph  input  from 

tape 

None 

NSTEPS 

Number  of  input  records  on  input 
file 

hydrooraph 

None 

QCARD 

Rate  of  inflow,  from  cards 

cfs 

QTAPE 

Rate  of  inflow,  from  tape 

cfs 

TE 

Time  of  inflow  for  card  input 

Hours/Seconds 

TEO 

Previous  value  of  TE 

Hours/Seconds 

TIMEO 

TEO 

Seconds 

TP 

TZERO 

Seconds 

T2 

Time  of  inflow  for  tape  input 

T20 

Previous  value  of  T2 

Seconds 

WATSH 

Not  used  at  this  time 

None 

This  common  block  is  used  by  the 
subroutines 
INDATA 
INFLOW 

foil  owing 
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Variable 

Name  Description 


Unit 


AS 

GRELEV 

JSKIP 

JUN 
M CHAN 

din 

oiriST 

oou 

SUMAL 

Y 

YT 


COMMON/ JUNC/ 


Surface  area  of  a node 

Square 

Ground  elevation  at  a node 

Feet 

Internal  integer  control,  to  skin  node 
head  computation 

None 

Node  number 

None 

Conduits  connecting  to  a node 

Hone 

Flow  into  a node  from  an  outside  source 

cfs 

Dry  weather  flow  into  a node  from  an  outside 
source 

cfs 

Flow  from  a node 

cfs 

Sum  of  cross  sectional  area/length,  for 
all  pipes  at  node 

Feet 

Depth  of  water  at  a node  at  full 
integration  step 

Feet 

Denth  of  water  at  a node  at  half  integration 
step 

Feet 

Z Elevation  of  node  invert  Feet 

ZCROUN  Elevation  of  uppermost  conduit  crown  Feet 

at  a node 


This  common  block  is  used  by  the  followino 
subroutines 
BOUND 
HEAD 
INDATA 
INFLOW 
MAIN 
OUTPUT 
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Variable 

Name 

Description 

Unit 

COMMON/LAB/ 

HORIZ 

Horizontal  label  of  curve 

None 

TITLE 

Title  printed  out  on  curve 

None 

VERT 

Vertical  label 

None 

XLAB 

Numerical  scale  labels  for  X 

None 

YLAB 

Numerical  scale  labels  for  Y 

None 

This  common  is  used  by  the  following 
subroutines 

IN DATA 
OUTPUT 
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Variable 

Name 


AORIF 
COR  IF 
LORIF 
NOR  IF 


CPRT 

I COL 

IDUM 

INTER 

IPRT 

JPLT 

JPRT 

KPLT 

LPLT 

LTIME 

NHPRT 


Description 


Unit 


COMMON/ORF/ 

Cross  sectional  area  of  orifice 
Orifice  coefficient 
Internal  number  index 
Number  of  orifices 

This  common  is  used  by  the  following 
subrouti nes 
BOUND 
INDATA 


Square  Feet 
Square  Feet 
None 
None 


COMMON/OUT/ 

Conduit  numbers  for  detailed  Drintinq  None 

Not  used  at  this  time  None 

Dummy  array  None 

Interval  between  print  cycles 
Not  used  at  this  time 


Node  numbers  for  plotting  None 
Node  numbers  for  detailed  printing  None 
Conduit  numbers  for  plotting  None 
Number  of  conduits  for  detailed  printinq  None 
Counter  for  printed  output  None 
Number  of  nodes  for  detailed  printing  None 
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Variable 

Name  Description 


Uni  t 


NPLT 

Number  of  nodes  to  be  plotted 

None 

NPRT 

Not  used  at  this  time 

None 

NPTOT 

Total  number  of  plot  data  points 

None 

NQPRT 

Number  of  conduits  for  detailed  printing 

None 

NSTART 

First  cycle  where  saved  printing 
array  will  begin 

None 

PRGEL 

Print  matrix,  ground  elevation 

Feet 

PRTH 

Print  matrix,  water  surface  elevation 

Feet 

PRTQ 

Print  matrix,  flow 

cfs 

PRTV 

Print  matrix,  velocity 

fps 

PRTY 

Print  matrix,  water  depth  at  node 

Feet 

QPLT 

Matrix  of  flow  values 

cfs 

YPLT 

Matrix  of  water  surface  elevations 

Feet 

This  common  is  used  by  the  following 
subroutines 
INDATA 
MAIN 
OUTPUT 

COMMON/ PI  PE/ 

A 

Full  integration  step  wetted  cross  section 

Square  Feet 

AFULL 

Cross  sectional  area  of  conduit 

Square  Feet 

AT 

Half  integration  step  wetted  cross  section 

Square  Feet 

TABLE  V- 3 
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Variable 

Name 

Description 

Unit 

DEEP 

Vertical  dimension  of  conduit 

Feet 

H 

Internal  depth  variable 

Feet 

LEN 

Conduit  lenqth 

Feet 

NCOND 

Conduit  number 

None 

NKLASS 

Conduit  shape  classification 

1 . circular 

2.  rectangular 

3.  horseshoe 

4.  egnshape 

5.  baskethandle 

6.  trapezoidal 

None 

ruiric 

Modes  at  each  end  of  conduit 

None 

Q 

F 1 ov/  in  conduit  at  full  integration  steD 

cfs 

QO 

Saved  flow  values 

cf  s 

OT 

Flow  in  conduit  at  half  integration  step 

cfs 

RFULL 

Hydraulic  radius  of  conduit  when  full 

Feet 

ROUGH 

Mann  inn  coefficient 

V 

Velocity  in  conduit  at  the  full  integration 
step 

fps 

VT 

Velocity  in  conduit  at  the  half  integration 
step 

f DS 

WIDE 

Width  of  conduit 

Feet 

ZP 

Height  of  conduit  invert  above  node  invert 

Feet 

V-89 


TABLE  V- 3 
(Continued) 


10  of  23 


Variable 

Name 


JPFIIL 

LPUMP 

NPUMP 

PRATE 

VRATE 

'/WELL 


INQUAL 


Descri ption 


This  common  block  is  used  by  the  following 
subroutines 
BOUND 
DEPTH 
HEAD 
HYDRAD 
TNDATA 
MAIN 
OUTPUT 


COMMON/PUMP/ 

Internal  inteqer  switch  for  full  wet  well 

0 = full 

1 = not  full 

Internal  pump  linkage 
dumber  of  Dumps 
Pumping  rate 

Volume  for  changinq  Dump  rates 
Starting  volume  of  pump  wet  well 


This  common  is  used  by  the  followina 
subroutines 
BOUND 
IN DATA 


COMMON/QUAL/ 

Dummy  variable  for  read  purposes 


This  common 
subroutines 
IN DATA 
INFLOW 


block  is  used  bv  the  following 


V-90 


Un  i t 


None 

None 

None 

cfs 

Cubic  Feet 
Cubic  Feet 


None 
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Variable 

Name  Description 


Unit 


AA 


SXX 

SXY 

TT 

XX 

YY 


SPHI 

STHETA 


N22 


COMMON/TIDE/ 

Tidal  curve  fit  coefficients  during  least  None 

souare  process 

Matrix  used  by  best  square  process  None 

Vector  used  by  least  souare  process  None 

Clock  time  of  tidal  condition  Hours/Seconds 

Vector  used  in  least  square  tide  fit  None 

Stage  level  of  tidal  input  Feet 


This  common  is  used  by  the  following 
subroutines 
INDATA 
TIDCF 


COMMON /TRAP/ 

Side  sloDe  one  side  trapezoidal  channel  Feet/Feet 

Side  slope  other  side,  trapezoidal  channel  Feet/Feet 


This  common  is  used  by  the  following 
subroutines 
HYDRAD 
IN DATA 
DEPTH 


COMMON/TRNQAL/ 

Tape  unit  number  of  hydraulic  output  for  None 

QUALITY  TRANSPORT  BLOCK 
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Variable 

Marne  Description 


Unit 


This  common  block  is  used  by  the  following 
subroutines 
INDATA 
MAIN 


COMMON/WEIR 

COEF 

Coefficient  of  discharge  for  weir 

None 

COEFS 

Coefficient  of  discharge  for  surcharged 
condition 

None 

KWEIR 

Type  of  weir 

1 . transverse 

2.  transverse  with  tide  nate 

3.  sideflow 

4.  sideflow  with  tide  gate 

None 

LWEIR 

Internal  link  number 

None 

NWEIR 

Number  of  weirs 

None 

WLEN 

Weir  length 

Feet 

YCREST 

Height  of  weir  above  invert 

Feet 

YTOP 

Heinht  to  top  of  weir  opening,  above  weir 

Feet 

This  common  is  used  by  the  following 
subroutines 
BOUND 
INDATA 

SUBROUTINE  BOUND 

AREA 

Wetted  cross  sectional  area  of  conduit 

Sauare  Feet 

DIR 

Internal  check  on  flow  direction 

None 
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Variable 

Name 

Descri ption 

Unit 

HEADO 

Internal  variable  for  head  at  orifice 

Feet 

HEADW 

Internal  variable  for  head  at  weir 

Feet 

HLOSS 

Loss  of  head,  weir  or  tide  gate 

Feet 

HRAD 

Hydraulic  radius 

Feet 

HS 

Internal  variable  for  head  at  submerged  weir 

Feet 

HTIDE 

Head  resulting  from  tide  computations 

Feet 

HI 

Internal  variable  for  nodal  water  surface 
el evation 

Feet 

H2 

Internal  variable  for  nodal  water  surface 
elevation 

Feet 

I 

Bookkeeping  integer 

None 

J 

Bookkeeping  integer 

None 

J1 

Internal  variable  for  node  number 

None 

J2 

Internal  variable  for  node  number 

None 

K 

Bookkeeping  integer 

KW 

Internal  variable,  type  of  weir 

None 

LINK 

Internal  1 ink  number 

None 

N 

Bookkeeping  integer 

None 

POWER 

Exponent  for  we ir  equation 

None 

QINJ 

Summation  of  inflows  to  pump  wet  well 

cfs 

QORIF 

Computed  flow  rate  throuah  orifice 

cfs 

QOUT 

Computed  pumped  flow 

cfs 
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Variable 

Name 

Descri oti on 

Uni  t 

VEL 

Velocity  through  a flap  gate 

fps 

VEL1 

Velocity  through  a tide  gate 

fps 

VNEU 

Pump  wet  well  stored  volume 

Cubic  Feet 

V 2 

Reserved  for  use  as  approach  velocity 
to  weir.  No  used  at  this  time 

fps 

WIDTH 

Conduit  width  at  computed  depth 

Feet 

UK 

Weir  coefficient 

None 

Y1 

Internal  deoth  variables 

Feet 

Y2 

Internal  depth  variables 

Feet 

This  subroutine  uses  the  following  common: 
/CONTR/ 

/ J'JNC/ 

/PIPE/ 

/ORF/ 

/WEIR/ 

/PUMP/ 

/BND/ 

SUBROUTINE  DEPTH 

AREA 

Wetted  cross  sectional  area  of  conduit 

Square  Feet 

DELTA 

Ratio  of  flow  difference 

None 

HRAD 

Hydraulic  radius 

Feet 

I 

Bookkeeping  integer 

None 

OC 

Internal  variable  for  critical  flow  rate 

cfs 
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Variable 

Name 

Descri pti on 

Uni  t 

QCO 

Internal  variable  for  last  computed 
critical  flow  rate 

cfs 

ONORM 

Internal  variable  for  normal  flow  rate 

cfs 

QNORMO 

Internal  variable  for  last  computed 
normal  flow  rate 

cfs 

This  subroutine  uses  the  following  common: 
/BD  / 

/PIPE/ 

/TRAP/ 

SUBROUTINE  HEAD 

AMID 

Wetted  cross  sectional  area  of  conduit 
at  mid  conduit 

Square  Feet 

BMID 

Width  of  water  surface  in  conduit  at 
mid  conduit 

Feet 

BHH 

Width  of  water  surface  in  conduit  at 
downstream  end 

Feet 

BNL 

With  of  water  surface  in  conduit  at 
upstream  end 

Feet 

RMID 

Hydraulic  radius  at  mid  conduit 

Feet 

RNH 

Hydraulic  radius  at  downstream  end 

Feet 

YC 

Critical  depth 

Feet 

YMID 

Depth  at  mid  conduit 

Feet 

YNH 

Depth  at  downstream  end 

Feet 
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Variable 

Name  Description 


Uni  t 


YNL 

Depth  at  upstream  end 

Feet 

YNORM 

Normal  depth 

Feet 

This  subroutine  uses  the  following  common: 
/JUNC/ 

/PIPE/ 

SUBROUTINE  HYDRAD 

DELTA 

Interpolation  variable 

None 

DEPTR 

Depth  in  conduit 

Feet 

DEPTT 

Depth  in  conduit 

Feet 

FDEP 

Difference  between  water  surface  depth 
and  full  conduit  water  surface  depth 

Feet 

FDEPTH 

Ratio  of  water  surface  depth  to  vertical 
conduit  dimension 

Feet/Feet 

I 

Bookkeepinn  integer 

None 

WETPER 

Wetted  perimeter  of  conduit 

Feet 

This  subroutine  uses  the  following  common: 
/BD/ 

/PIPE/ 

/TRAP/ 

SUBROUTINE  INDATA 

D1 

Dummy  variable 

None 

D2 

Dummy  variable 

None 
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Variable 

Name  Description 


Unit 


D3 

Dummy  variable 

lone 

I 

Bookkeeping  integer 

Done 

II 

Bcokkeeninq  integer 

.lone 

1 

Bookkeeping  integer 

None 

j'i 

Internal  orint  integer 

N^ne 

JP 

Internal  integer,  node  number 

flcne 

J1 

Internal  integer,  node  number 

ilcne 

T2 

Internal  integer,  node  nunher 

None 

1/ 
l * 

Bookkeeping  integer 

None 

1/1/ 

Bookkeeping  integer 

None 

KLASS 

Conduit  shape  classification 

None 

KO 

Control  integer  on  tide  data 

ilcne 

• / 1 

•\  I 

Conduit  number 

None 

L 

Bookkeeping  integer 

None 

LDC 

Bookkeepina  integer 

ilcne 

::pt 

Bookkeeping  integer 

None 

• ! 

Bookkeeping  integer 

None 

NAME 

Tape  element  name 

None 

riCHTID 

Print  control  on  tide  data 

Hone 
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Vari able 
Name 


NH 

NHR 

NJ1 

NJ2 

NL 

nmn 

NTCYS 

N1 

TEMP 

TPT 


J 

L 


Descri ption 


Node  number,  downstream  node 
Start  time 

Control  for  internally  generated  nodes 

Control  for  internally  generated  nodes 

Node  number,  upstream  node 

Start  time 

NTCYC  plus  one 

Bookkeeping  integer 

Internal  variable,  conduit  invert  elevation 
Time  used  in  printing 


This  subroutine  uses  the  following  common: 
/BD/ 

/CONTR/ 

/TRNQAL/ 

/JUNC/ 

/PIPE/ 

/ORF/ 

/WEIR/ 

/PUMP/ 

/BND/ 

/OUT/ 

/TIDE/ 

/HYFLOW/ 

/LAB/ 

/QUAL/ 


SUBROUTINE  INFLOW 
Bookkeeping  integer 
Bookkeeping  integer 
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Unit 


None 

Hours 

None 

None 

None 

Minutes 

None 

None 

Feet 

Hours 


None 

None 
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Variable 

Name  Description 


Unit 


Q1 

Nodal  inflow  rates  from  external  sources 

cfs 

02 

Nodal  inflow  rates  from  external  sources 

cfs 

SLOPE 

Rate  of  chanqe  of  inflow  rates 

Cu. Ft. /^ec. squared 

TPT 

Time 

Hours 

This  subroutine  uses  the  followinq  common: 
/CONTR/ 

/JUNC/ 

/OUAL/ 

/HYFLOW/ 

SUBROUTINE  MAIN 

AKON 

Internal  variable 

None 

ANH 

Wetted  cross  sectional  area  of  conduit  at 
downstream  end  of  conduit 

Square  Feet 

ANL 

Wetted  cross  sectional  area  of  conduit  at 
upstream  end  of  conduit 

Square  Feet 

DELQ1 

Component  of  flow  change  during  a time  step 

cfs 

DEL02 

Component  of  flow  change  durinq  a time  step 

cfs 

DELQ3 

Component  of  flow  chanqe  durinq  a time  step 

cfs 

DELQ4 

Component  of  flow  change  durinq  a time  step 

cfs 

HRAD 

Hydraulic  radius 

Feet 

I 

Bookkeeping  inteqer 

None 

J 

Bookkeeping  integer 

None 
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Variable 

Name  Description 


Unit 


K 

Bookkeepina 

integer 

None 

L 

Bookkeeoing 

inteaer 

None 

M 

Bookkeepina 

integer 

None 

'ICY 

Bookkeepina 

integer 

None 

MCYY 

Bookkeepina 

integer 

None 

'IP 

Bookkeeping 

integer 

None 

M 

Bookkeepina 

inteqer 

None 

NDIM 

Bookkeeping 

integer 

None 

NH 

Internal  downstream  node  number 

None 

ONEW 

Last  computed  flow 

cfs 

QNORM 

Normal  flow 

in  conduit 

cfs 

SUMO 

Summation  of 

flow  at  node,  usinn 

average 

cfs 

conduit  flow 

rates 

SUHQS 

Summation  of 

flow  a*  lode,  usina 

1 ast 

cfs 

computed  flow  rates 

YMAX 

Height  at  node  from  invert  to  highest 

Feet 

conduit  crown 

This  subroutine  used  the  following  common: 

/CONTR/ 

/TRNQAL/ 

/ JUNC/ 

/PIPE/ 

/OUT/ 

/TRAP/ 
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Variable 

Name 

Description 

Unit 

I 

SUBROUTINE  OUTPUT 
Bookkeening  integer 

None 

IT 

Bookkeepina  integer 

None 

J 

Bookkeeping  integer 

None 

K 

Bookkeeping  inteoer 

None 

L 

Bookkeeping  integer 

None 

LT 

Bookkeeping  integer 

None 

LTIMEH 

Time  variable  for  printing 

Hours 

LTIME‘1 

Time  variable  for  printing 

Minutes 

MJPRT 

Node  number  for  printina 

None 

N 

Bookkeeping  variable 

None 

NKON 

Conduit  number  for  printing 

None 

TIMEO 

Start  time  of  printing 

Seconds 

This  subroutine  uses  the  following  common: 
/CONTR/ 

/JUNC/ 

/PIPE/ 

/OUT/ 

/ELEV/ 

/LAR/ 
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Variable 

Name  Description 


Unit 


SUBROUTINE  TIDCF 


DEL 

Internal  variable 

None 

DELMAX 

Maximum  difference  between  the  calculated 
and  tidal  staqe  input 

Feet 

DELTA 

Maximum  allowable  difference  between  the 
calculated  and  input  tidal  stage 

Feet 

DIFF 

Difference  between  the  calculated  and 
input  tidal  stage 

Feet 

FJ1 

Internal  variable 

None 

FJ3 

Internal  variable 

None 

I 

Bookkeeping  integer 

None 

IT 

Bookkeeping  integer 

None 

J 

Bookkeeping  integer 

None 

K 

Bookkeeping  integer 

None 

MAX  IT 

Maximum  number  of  iterations 

None 

NJ2 

Bookkeeping  integer 

None 

NTT 

Bookkeeping  integer 

None 

PERIOD 

Daily  tidal  cycle  time  in  hours 

Hours 

RES 

Accumulative  difference  between  the 
calculated  and  input  tidal  staqe 

Feet 
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Variable 

Name 


SUM 


Description 


Computed  tidal  stage 


This  subroutine  uses  the  followino  common: 
/CONTR/ 

/TIDE/ 
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Variable 

Name 

Description 

Unit 

COMMON/ BD/ 

ANORM 

Matrix  of  normalized  wet  cross  sectional 
area  of  conduit,  based  on  shape  and  depth 

None 

HRNORM 

Matrix  of  normalized  hydraulic  radius  of 
conduit,  based  on  shape  and  depth 

None 

TWNORM 

Matrix  of  normalized  conduit  width  at  flow 
line,  based  on  shape  and  depth 

None 

This  common  is  used  by  the  following  subroutines 
HYDRAD 
BLOCK  DATA 

COMMON/CNTROL/ 

ALPHA 

.Lot  used  at  present  time 

None 

BETA 

Description  of  Quality  Transport  simulation 

None 

DELT 

Time  interval  of  integration 

Seconds 

INTER 

Number  of  time  intervals  between  print 
data  points 

None 

IQCYC 

Boc..,<eeping  integer,  counter 

None 

NINREC 

Bookkeeping  integer  on  tape  reads 

None 

NPPNT 

Number  of  nodes  for  detailed  printout 

None 

*lnCYC 

Number  of  integration  cycles 

None 

10UAL 

Number  of  guality  constituents 

None 

Cycle  for  starting  printout 

None 
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Variable 

Name 

Description 

Unit 

N21 

Taoe  unit  number,  pollutograph  data 

None 

N22 

Tape  unit  number,  hydraulics  data 

None 

TE 

Current  time  from  hydraulic  input  tape 

Seconds 

TEO 

Last  current  time  from  hydraulic  input  tape 

Seconds 

TIME 

Time  from  start  of  simulation  plus  TZERO 

Seconds 

TIMEO 

Same  as  TEO 

Seconds 

TITLE 

Description  of  results  from  Runoff  Block  input 

None 

TP 

Same  as  TZERO 

Seconds 

TZERO 

Start  time  of  simulation 

None 

T2 

Current  time  from  pollutoaraoh  input  tape 

Seconds 

T20 

Last  current  time  from  pollutoaraph  input 
tape 

Seconds 

This  common  is  used  by  the  following  subroutines 
INDATA 
INPUT 
MAIN 
OUTPUT 

COMMON /COM D/ 

AFULL 

Cross  sectional  area  of  conduit 

Square  Feet 

DEEP 

Vertical  dimension  of  conduit 

Feet 

LEN 

Length  of  conduit 

Feet 

NC 

Number  of  conduits 

None 

V- 105 

i 

I 


TABLE  V- 4 

(Continued)  3 of  11 


Variable 

Name 

Descri pti on 

Uni  t 

NCCND 

Conduit  number 

None 

NJUNC 

Nodes  at  each  end  of  conduit 

None 

NKLASS 

Conduit  shape  classification 

None 

NTC 

Number  of  links  in  system  (NC  plus  internal) 

None 

Q 

Flow  in  conduit 

cfs 

QAVE 

Current  flow  in  conduit 

cfs 

RFULL 

Hydraulic  radius  of  conduit  when  full 

Feet 

V 

Velocity  in  conduit 

fps 

WIDE 

Width  of  conduit 

Feet 

This  common  is  used  by  the  following  subroutines 
HYDRAD 
IN DATA 
INPUT 
MAIN 
OUTPUT 

COMMON/JUNCT/ 

C 

Constituent  concentration 

mg/1 

JSKIP 

Control  inteaer  for  pumps 

None 

JUN 

Mode  number 

None 

MASS  IN 

Mass 

mg/1 -cu . ft. 

MCHAN 

Conduits  connecting  to  a node 

None 

NJ 

Number  of  nodes 

None 
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Variable 

Name 

Descri ntion 

Unit 

VOL 

Current  storaqe  volume  at  node 

Cu. Ft. 

V0L0 

Last  current  storaqe  volume  at  node 

Cu. Ft. 

Y 

DeDth  of  water  at  node 

Feet 

YAVE 

Current  depth  of  water  at  node 

Feet 

This  common  is  used  by  the  followina  subroutines 
INDATA 
INPUT 
MAIN 
OUTPUT 

COMMON/OUT/ 

CONST 

Constituent  name 

None 

CPRNT 

Computed  constituent  concentrations 

mg/1 

I CODE 

Plot  code  indicator  for  constituents 

ICODL 

Print  code  indicator  for  constituents 

JQPRT 

Node  numbers  for  detailed  printout 

None 

LTIME 

Counter  for  printed  output 

None 

NPLTL 

Lower  limit  of  plot  information  in 
output  storage  arrays 

NPLTM 

Upper  limit  of  plot  information  output 
storage  arrays 

QPRNT 

Nodal  flow  unbalance 

cfs 
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Variable 

Name  Description 


Unit 


This  common  is  used  by  the  following  subroutines 
IN DATA 
INPUT 
MAIN 
OUTPUT 
BLKD 

COMMON/PUMP/ 


LPUMP 

Internal  pump  linkage 

None 

NPUMP 

Number  of  Dumps 

None 

VPUMP 

Pump  wet  well  volume 

Cu.Ft. 

This  common  is  used  by  the  followinq  subroutines 
INDATA 
INPUT 
MAIN 
OUTPUT 

COMMON/RUNOFF/ 

ISW 

Node  numbers  for  pollutograph  inputs 

None 

HJSW 

Number  of  nodes  for  pollutograph  inputs 

None 

NSTEPS 

Number  of  data  points  on  pollutograph 
input  tape 

None 

WATSH 

Watershed  number 

None 

WDMDT 

Mass  input  rate 

grams/sec 

This  common  is  used  by  the  followinq  subroutines 
INDATA 
INPUT 
MAIN 
OUTPUT 
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Variable 

flame 

Descri otion 

Unit 

COfIMON/TRAP/ 

SPHI 

Side  slope  of  trapezoidal  channel,  one  side 

Feet/Feet 

STHETA 

Side  slope  of  trapezoidal  channel,  other  side 

Feet/Feet 

This  common  is  used  by  the  following 
subroutines 
HYDRAD 
INDATA 

SUBROUTINE  HYDRAD 

DELTA 

Interpolation  variable 

None 

DEPTR 

Depth  in  conduit 

Feet 

DEPTT 

Depth  in  conduit 

Feet 

FDEP 

Difference  between  water  surface  depth 
and  full  conduit  water  surface  depth 

Feet 

FDEPTH 

Ratio  of  water  surface  depth  to  vertical 
conduit  dimension 

Feet/Feet 

I 

Bookkeeping  integer 

None 

WETPER 

Wetted  perimeter  of  conduit 

Feet 

This  subroutine  uses  the  followinq  common: 
/BD/ 

/COND/ 

/TRAP/ 


V- 109 


i 


TABLE  V-4 


(Continued) 

7 of  11 

Variable 

Name 

Description 

Uni  t 

D1 ,D2  ,D3 

SUBROUTINE  INDATA 
Dummy  variable 

None 

I 

Bookkeeping  integer 

None 

J 

Bookkeepinq  integer 

None 

K 

Bookkeeping  integer 

None 

L 

Bookkeeping  integer 

None 

MIN 

Start  time 

Minutes 

N 

Bookkeeping  integer 

None 

NAMER 

Not  used  at  this  time 

None 

NAMES 

Not  used  at  this  time 

None 

NHR 

Start  time 

Hours 

NTCYC 

Total  number  of  time  steps 

None 

TD 

Dummy  variable 

None 

TMAX 

Maximum  time  simulated 

Seconds 

TMAXD 

Maximum  simulation  time  from  hydraulic  tape 

Seconds 

TMAXW 

Maximum  simulation  time  from  pollutoqraph 

Seconds 

TTMAX 

tape 

Minimum  value  between  TMAXW  and  TMAXD 

Seconds 

TTEMP 

Internal  variable 
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Variable 

Name  Description 


Unit 


This  subroutine  uses  the  following  common: 
/CNTROL/ 

/COND/ 

/juriCT/ 

/OUT/ 

/PUMP/ 

/RUNOFF/ 

/TRAP/ 


I 

J 

L 

N 

Ql 

Q2 

SLOPEQ 

SLOPEY 

TD 

W1 

W2 

Y1 

Y2 


SUBROUTINE  INPUT 
Bookkeepinq  integer 
Bookkeeping  integer 
Bookkeeping  integer 
Bookkeeping  integer 
Current  flow  rate  in  conduit 
Last  current  flow  rate  in  conduit 
Rate  of  change  of  flow 
Rate  of  change  of  head  at  node 
Dummy  variable 
Current  mass  input  rate 
Last  current  mass  input  rate 
Current  head 
Last  current  head 


None 

None 

None 

None 

cfs 

cfs 

cfs/sec 

ft/sec 

mq/l-cu.ft. 

mg/l-cu.ft. 

Feet 

Feet 
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Variable 

Name 

Descri ption 

Unit 

1 

This  subroutine  uses  the  following  common: 
/CNTROL/ 

/COND/ 

/ JUMCT / 

/OUT/ 

/PUMP/ 

/RUNOFF/ 

SUBROUTINE  MAIN 

A1 

Wetted  cross  sectional  area  of  conduit, 
upstream  end 

Square  Feet 

A2 

Wetted  cross  sectional  area  of  conduit, 
downstream  end 

Square  Feet 

GRAD 

Concentration  gradient  in  conduit 

mg/1 -ft. 

HRAD 

Hydraulic  radius 

Feet 

I 

Bookkeeping  integer 

None 

II 

Bookkeeping  integer 

None 

J 

Bookkeeping  integer 

None 

JFROM 

Upstream  node  number 

None 

JPUMP 

Pump  node  number 

None 

J1 

Upstream  node  number 

J2 

Downstream  node  number 

K 

Bookkeeping  integer 

None 

L 

Bookkeeping  integer 

None 

MCY 

Bookkeeping  integer 

None 

k. 
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Variable 

Name  Description 


Uni  t 


f) 

Bookkeeping  integer 

None 

Ml 

Bookkeeping  integer 

None 

SU’IQ 

Sum  of  flows  at  a node 

cfs 

VPUMPO 

Pump  wet  well  volume 

Cubic  Feet 

This  subroutine  uses  the  following  coimion: 
/CNTPnL/ 

/COND/ 

/JtlNCT/ 

/OUT/ 

/PUMP/ 

/RUNOFF/ 

SUBROUTINE  OUTPUT 

I 

Bookkeeping  integer 

None 

IT 

Bookkeeoinq  integer 

None 

0 

Subscript  integer 

None 

KC 

Bookkeeping  integer 

None 

L 

Bookkeeping  integer 

None 

LL 

Bookkeening  integer 

None 

LT 

Bookkeeping  integer 

None 

LTIMEH 

Time  for  print  purposes 

Hours 

LTIMENI 

Time  for  print  purposes 

Minutes 
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Variable 

flame 


Description 


This  subroutine  uses  the  followina  common: 
/CNTROL/ 

/COND/ 

/JUNCT/ 

/OUT/ 

/PUMP/ 

/RUNOFF/ 
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CHAPTER  VI 

CRITERIA  APPLICABLE  TO  THE  RIBCO  STUDY  AREA 


INTRODUCTION 

The  urban  stormwater  models  will  be  used  by  counties,  munici- 
palities, local  agencies,  and  consultants  as  tools  for  evaluating  existing 
and  future  drainage  systems,  and  for  comparing  alternative  courses  of 
action.  Presented  in  this  section  are  specific  rainfall,  infiltration, 
and  receiving  water  criteria  that  are  representative  of  the  Green  and 
Cedar  River  Basins  and  should  be  used  as  input  to  the  models. 

In  addition  to  the  above  criteria,  recommendations  are  presented 
for  the  planning  of  storm  drainage  facilities.  These  recommendations 
are  based  upon  the  design  standards  presently  being  followed  by  most  of 
the  governmental  agencies  within  the  RIBCO  study  area. 

Both  the  criteria  and  the  recommendations  were  reviewed  and 
approved  by  the  Urban  Runoff  and  Basin  Drainage  Technical  Review  Sub- 
committee and  by  RIBCO. 


RAINFALL 


The  time  distribution  of  rainfall  in  a storm  is  of  considerable 
importance  in  calculating  the  resulting  runoff.  This  is  especially  true 
when  utilizing  the  urban  stormwater  models  to  simulate  runoff. 

Rainfall  records  were  collected  and  analyzed  for  the  stations 
located  in  and  adjacent  to  the  Cedar  and  Green  River  Basins.  The  results 
indicate  that  a generalized  rainfall  rate  versus  time  relationship  is 
applicable  to  all  lands  within  the  Basins  at  elevations  below  500  feet 
that  are  located  west  of  a north-south  line  through  Issaquah  and  Enumclaw. 
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The  criteria  presented  herein  are  not  applicable  to  lands 
above  elevation  500  feet,  as  a large  portion  of  the  precipitation  in  the 
upper  areas  of  the  Cedar  River  and  Green  River  Basins  is  in  the  form  of 
snow,  and  frozen  ground  condi tons  are  prevalent. 


The  planning  of  storm  drainage  facilities,  as  discussed  in 
subsequent  pages,  is  based  upon  the  runoff  resulting  from  the  rainfall 
having  a 10-year  recurrence  interval.  Therefore,  rainfall  for  the  10- 
year  recurrence  interval  is  as  follows: 


Total  Rainfall 

(inches) 

Within  the  City 

Remainder  of 

Duration 

of  Seattle 

Study  Area 

10  minutes 

0.27 

0.33 

20  minutes 

0.36 

0.44 

30  minutes 

0.41 

0.52 

1 hour 

0.55 

0.67 

2 hours 

0.74 

0.88 

4 hours 

1.00 

1.16 

6 hours 

1.20 

1.40 

12  hours 

1.78 

1.98 

24  hours 

2.40 

2.80 

For  all  sub-basins  within  the  study  area  a four-hour  duration 
storm  will  generate  the  highest,  or  almost  the  highest,  peak  rate  of 
runoff.  It  is  recommended,  therefore,  that  a four-hour  duration  storm 
be  used  in  the  computer  model.  If  desired,  larger  or  shorter  duration 
storms  can  be  used  and  the  results  compared  against  the  results  of  a 
four-hour  storm. 

U.S.  Weather  Bureau  technical  papers  can  be  used  for  obtaining 
rainfalls  with  other  than  a 10-year  recurrence  interval.  The  rainfall 
distribution  versus  time  relationship  for  storms  with  greater  recurrence 
intervals  will  be  approximately  the  same  distribution  as  that  listed 
above  for  the  10-year  recurrence  interval. 
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INFILTRATION 


The  surface  infiltration  rates  of  the  soils  within  the  Green 
and  Cedar  River  Basins,  except  forested  lands  at  high  elevations,  have 
been  measured  by  the  USDA  Soil  Conservation  Service  and  are  appropriate 
for  use  in  the  stormwater  models. 

The  Soil  Conservation  Service  infiltration  rates,  high,  medium 
or  low,  for  each  soil  type,  have  been  compiled  on  maps  at  a scale  of  1"  = 
2000  ft.  by  the  Puget  Sound  Governmental  Conference.  Both  initial  (maxi- 
mum) infiltration  rate  and  base  (minimum)  infiltration  rate  are  used  in 
the  computer  model.  The  rate  of  decay  is  0.00115  for  all  soil  types. 
Infiltration  of  rainfall  on  impervious  areas  is  to  be  considered  negligi- 
ble. 


Infiltration  rates  by  land  use  classification  is  as  follows: 


Intake  Rate 

(Inches 

per  Hour) 

Land  Use 

from  Maps 

Maximum 

Mi nimum 

Parks/Dedicated  Open  Space 

High 

3.0 

0.7 

Agriculture 

Medium 

2.5 

0.4 

Unused  Land 

Low 

2.0 

0.2 

Single-Family  Residential 

High 

2.5 

0.4 

Med i urn 

2.0 

0.2 

Low 

2.0 

0.2 

Multi-Family  Residential 

High 

2.0 

0.2 

Commercial /Services 

Medium 

2.0 

0.2 

Government  and  Education 

Low 

2.0 

0.2 

Industrial 

As  the  land  use  becomes  more  intense  intake  rates  are  reduced  from  the 
natural  rates  (i.e.  those  indicated  for  Parks/Dedicated  Open  Space, 
Agriculture,  Unused  Land).  This  reduction  in  intake  rates  is  due  to  a 
greater  compaction  of  the  surface  soils  as  a result  of  clearing  and  grad- 
ing operations  associated  with  the  construction  of  homes,  stores,  indus- 
trai 1 taci lities,  etc. 
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ENVIRONMENTAL  PLANNING  FOR  THE  METROPOLITAN  AREA  CEDAR-GREEN  RI--ETC (U) 
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In  addition  to  the  intake  rate,  certain  soils  within  the  study 
area  are  limited  as  to  the  total  volume  of  water  that  can  be  infiltrated. 


The  following  soils,  as  identified  by  the  Soil  Conservation  Service, 


have  a maximum  to 

two  inches  of 

infiltration 

capaci ty : 

Symbol 

Soil  Type 

Symbol 

Soil  Type 

Br 

Briscot 

Tu 

Tukwi la 

Nk 

Nooksack 

Bh 

Bell  ingham 

No 

Norma 

Bu 

Buckley 

Or 

Orcas 

0s 

Oridia 

Sh 

Sammamish 

0 

Oval  1 

Sk 

Seattle 

Pu 

Puget 

Sm 

Shalcar 

Su 

Sultan 

So 

Snohomish 

Wo 

Woodinvilli 

Sr 

Snohomish 

(thick) 

RECEIVING  WATERS 

All  storm  drainage  facilities  in  the  Green  and  Cedar  River 
Basins  empty  into  one  of  the  following  six  receiving  water  bodies.  There- 
fore, all  facilities  must  be  able  to  discharge  storm  runoff  against  the 
natural,  or  man-controlled,  water  surface  elevations  indicated  herein. 


Lake  Washington 


Water  Surface  Elevations  (USGS  datum) 
15.0  ft. 


Lake  Sammamish 


32.5  ft. 


Puget  Sound  8.2  ft. 

Green  Ri ver/Duwamish  River  The  water  surface  profile  labeled 

"12,000  cfs  at  Auburn"  on  28  x 40"  plan 
and  profile  sheets  2 of  12  through  12 
of  12.  This  represents  a recurrence 
interval  of  approximately  300  years  and 
includes  SCS  pumping  system.  These 
drawings  were  prepared  by  the  Corps  and 
are  entitled  Green  & Duwamish  Rivers, 
Wash.,  Survey  Report  Profiles,  May  1968. 
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Water  Surface  Elevations  (USGS  datum) 


Cedar  River  The  water  surface  profile  labeled  "100 

year  flood"  on  28"  x 40"  plan  and  pro- 
file sheets  1 of  7 through  7 of  7. 

These  drawings  were  prepared  by  the 
Corps  and  are  entitled  "Water  Surface 
Profiles  for  Flows  Limited  to  the  Flood- 
way," dated  1 August  1970. 

Sammamish  River  The  water  surface  profile  labeled  "100 

year"  on  11"  x 22"  profile  sheets,  two 
total.  These  drawings  were  prepared  by 
the  Corps  and  are  entitled  "Sammamish 
River  Profiles  Mile  0.00  to  Mile  14.95," 
dated  May  1974. 

NOTE:  Corps  of  Engineers’  datum  is  6.33  ft.  below  the  USGS  datum. 

(e.g.  12.00  ft.  Corps  = 5.67  ft.  USGS) 


STORM  DRAINAGE  FACILITIES 


A.  General 


All  storm  drainage  facilities  should  be  planned  to  conform 
to  the  recommendations  set  forth  herein. 


All  storm  drainage  facilities  shall  be  planned  to  accommo- 
date the  runoff  resulting  from  the  rainfall  having  a 10-year  recurrence 
interval . 

Watercourses  and  storm  sewers  having  a drainage  area  of 
over  four  square  miles  are  classed  as  major  facilities. 

Watercourses  and  storm  sewers  having  a drainage  area  of 
less  than  four  square  miles  are  classed  as  minor  facilities. 

A given  watercourse  or  storm  sewer,  therefore,  will  be 
classed  as  minor  in  its  upper  reaches,  then  changed  to  the  major  classi- 
fication at  a point  where  the  drainage  area  exceeds  four  square  miles. 
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Depth  of  flow  in  gutters  shall  not  exceed  0.4  feet.  Road- 
side ditches*  are  allowed,  except  that  they  shall  not  be  used  where  the 
flow  is  greater  than  that  which  could  be  carried  in  a standard  gutter 
flowing  0.4  feet  deep  on  the  same  slope  as  the  road  profile  slope.  Where 
the  discharge  exceeds  gutter  capacity,  a storm  sewer  shall  be  provided. 
Roadside  ditches  shall  be  planned  so  that  the  water  surface  will  be  at 
or  below  the  outside  edge  of  the  road  shoulder  and  below  adjacent  ground 
level . 


Insofar  as  is  practicable,  catch  basins,  manholes,  inlet 
structures,  etc.,  shall  conform  to  standard  plans  of  the  public  agency 
with  jurisdiction  where  the  improvement  is  located. 

B.  Manning's  "n"  Values 

Manning's  "n"  values  (coeffi  of  roughness)  for  plan- 

ning shall  be  as  follows: 


1. 

Concrete,  steel  troweled  or  smooth- 
form  finish 

n = 

.013 

2. 

Concrete  pipe,  precast  of  cast-in-place 

n = 

.014 

3. 

Concrete,  wood  flat,  or  broomed  finish; 
including  pneumatically  applied  mortar 

n = 

.015 

4. 

Asphaltic  concrete 

n = 

.017 

5. 

Corrugated  metal  pipe 

n = 

.024 

6. 

Sack  concrete  riprap 

n = 

.030 

7. 

Grouted  rock  riprap 

n = 

.040 

8. 

Loose  rock  riprap 

n = 

.050 

For  grassed  channels,  natural  channels,  or  materials  not 
covered  above,  appropriate  "n"  values  shall  be  determined. 


♦Stream  channels  parallel  to  roads  are  not  classed  as  roadside  ditches. 
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C.  Channels* 


Channels,  as  referred  to  herein,  exclude  the  main  stems  of 
the  Cedar  River,  Green  River/Duwamish  River,  and  the  Sammamish  River. 

Major  channels  shall  be  planned  with  a minimum  freeboard 
between  the  water  surface  and  the  top  of  bank  of  1.0  foot.  Minor  channels 
may  be  planned  without  any  freeboard.  The  use  of  levees  for  control  of 
channel  flow  will  be  discouraged. 

For  natural  watercourses  flow  may  be  allowed  in  the  floodway 
area  above  the  defined  banks.  However,  that  flow  will  be  contained  with- 
in a defined  overflow  area,  and  freeboard  as  specified  above  shall  be 
provided  between  the  water  surface  and  adjacent  ground  elevation  or  fin- 
ished grade  elevation. 

Culverts  along  major  channels  shall  be  planned  to  flow  full, 
but  not  under  pressure,  and  shall  cause  no  encroachment  on  the  specified 
minimum  freeboard  in  the  upstream  channel  or  waterway.  Culverts  along 
minor  channels  may  be  planned  for  full  capacity  and  pressure  flow,  but 
shall  cause  no  encroachment  on  the  specified  minimum  freeboard  in  the 
upstream  channel  or  waterway. 

Minimum  bottom  width  of  artificial  channels  shall  be  2.0 
feet.  Channels  may  be  lined  with  concrete,  grouted  rock  riprap,  sack 
concrete  riprap,  or  air-blown  mortar.  Grassed  channels  or  loose-rock 
riprapped  channels  shall  have  side  slopes  not  steeper  than  2 to  1 . 

D.  Storm  Sewers 


Major  facilities  placed  within  a storm  sewer  system  shall 
be  planned  to  flow  full,  but  not  under  pressure.  Minor  facilities  placed 


*In  conjunction  with  the  criteria  set  forth  herein,  it  is  recommended  that 
all  local  agencies  require  building  sites  to  be  one  to  two  feet  above 
planned  water  surfaces  elevations. 


in  storm  sewers  may  be  planned  for  full  conduit  capacity  and  pressure 
flow.  At  inlets  and  non-pressure  type  manholes  within  a storm  sewer  system, 
the  water  surface  shall  be  1.0  foot  or  more  below  the  gutter  or  inlet 
surface  elevation. 

E.  Storage  Facilities 


Storage  facilities  can  be  natural  lakes  or  ponds,  they  can 
be  all  or  partially  man-made,  and  they  can  be  located  either  offstream 
or  onstream. 


Embankments  and  holding  ponds  located  near  or  on  channels 
must  be  specifically  sited  and  designed.  However,  small-sized  (100  acre- 
feet  or  smaller)  holding  ponds,  other  than  natural  areas,  will  be  consi- 
dered as  follows: 

Depth  of  water  above  adjacent  ground:  8 ft.  maximum 

Perimeter  embankment:  12  ft.  top  width 

3 horz.  to  1 vert, 
side  slopes  con- 
structed of  earth 
material 

The  perimeter  embankments  shall  be  grass  lined  on  both  the 
pond  side  and  the  outside. 

All  holding  ponds  shall  have  an  uncontrolled  overpour  spill- 
way to  keep  the  pond  from  overtopping  the  embankment.  The  top  of  the  em- 
bankment shall  be  1.5  feet  above  the  maximum  water  surface  elevation 
within  the  holding  pond. 

Outlets  shall  be  provided  in  all  holding  ponds  to  release 
incoming  flows  to  the  downstream  channel  at  retarded  rates,  but  not  greater 
than  the  capacity  of  the  downstream  facilities.  Outlets  also  can  be  used 
to  release  temporarily  stored  waters,  and  to  allow  the  holding  pond  to  be 
emptied. 
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F.  Pumping  Facilities 


Pumping  facilities  may  be  used  with  detention  and  retention 
storage  ponds.  In  planning  retention  storage  facilities,  the  use  of  gravity 
flow  to  release  as  much  water  as  possible  after  the  storm  passes  is  most 
desirable.  However,  pumping  facilities  are  usually  located  in  low  portions 
of  levee-protected  areas  where  pumping  is  required  to  lift  the  water  to 
higher  elevations  such  as  discharging  to  a river. 

After  a retention  storage  pond  volume  is  selected,  pumps  shall 
be  sized  to  pump  all  subsequent  inflow  water  from  the  pond  to  the  receiving 
area. 

Pumps  shall  be  capable  of  pumping  against  the  most  severe 
static  head  (usually  low  pond  water  level  to  peak  river  stage). 
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CHAPTER  VII 

COST  DETERMINATION  PROGRAM 


INTRODUCTION 

This  chapter  describes  a computer  program  which  is  designed 
to  facilitate  the  calculation  of  various  costs  associated  with  urban 
drainage  problems.  The  primary  function  of  the  cost  analysis  program  is 
to  calculate  preliminary  capital  costs  of  construction  of  storm  drainage 
facilities  and  the  generalized  cost  of  purchasing  land  required  for  the 
facil ities. 

The  program  is  formulated  to  include  four  broad  categories  of 
facilities  construction  (see  Figure  VII-1).  These  are: 

1 . Open  channel s 

2.  Pipe  inlets  and  outlets 

3.  Holding  ponds 

4.  Pipes 

This  program  is  intended  to  be  used  in  conjunction  with  the  stormwater 
simulation  models,  and  can  be  employed  to  estimate  preliminary  costs  of 
various  drainage  alternatives. 

The  program  was  developed  for  use  in  preliminary  cost  estimating. 
Several  mathematical  computations  have  been  shortened  to  approximate 
versions  for  convenience.  The  total  estimated  capital  cost  is  heavily 
dependent  upon  the  unit  costs  and  these  formula  modifications  have  little 
effect  upon  the  total  costs. 


UNIT  COSTS  FOR  THE  SEATTLE  AREA 

The  model  is  programmed  to  calculate  construction  quantities 
and  then  to  apply  unit  costs  to  determine  totals.  The  following  unit 
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costs  have  been  compiled  from  actual  construction  costs  in  the  Seattle 
area.  The  costs  are  based  upon  the  Engineering  News  Record  (ENR) 
construction  cost  index  of  1450,  because  the  basic  portions  of  the  pro- 
gram related  to  pipes  and  unit  costs  were  originally  developed  for  use 
on  a METRO  project  in  1971. 


1.  Open  Channel  Excavation  $3.30  per  cubic  yard 

(off-site  haul  of  less  than  one  mile) 
includes:  clearing  and  grubbing 

coffer  dam  and/or  diversion  pipe 

fish  protection 

erosion  control 

final  trimming 

clean-up 


2.  ,ning 

ies  thick,  reinforced 

3.  Loose  Rock  Rip-Rap 

8 to  15  inch  diameter  rock,  placed  an 
average  of  12  inches  thick 
2. 

4.  Pipe  Costs 

A.  Class  III  Concrete  Pipe 
Class  IV  Concrete  Pipe 
Class  V Concrete  Pipe 

B.  Excavation 

C.  Gravel  backfill 


$100.00  per  cubic  yard 


1. 

$12.50  per  square  yard 


Cost  in  dollars 
per  foot  of  length 

0.0796  D ]-50 

0.0736  D ^ 

0.0728  D 1-58 

$2.50  per  cubic  yard 

Cost  in  dollars 
per  foot  of  length 


4.0 

"27 


D 


Combining  A,  B and  C together  with  cost  of  installation,  earth  backfill, 
catch  basins  or  manholes  at  500  feet  spacing,  site  restoration  and  clean-up 
the  following  total  installed  cost  formulas  were  derived: 


1.  Based  upon  small  quantities  placed  primarily  at  isolated  locations 
along  natural  streams. 

2.  D = Diameters  are  in  inches 

B = Backfill  depths  are  in  feet 
Z = Trench  depths  are  in  feet 
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Class  III  Concrete  Pipe 
Class  IV  Concrete  Pipe 
Class  V Concrete  Pipe 


Cost  in  dollars  per  foot  of  length 
for  trench  depths  12  feet  of  less 


Class  III  Concrete  Pipe 
Class  IV  Concrete  Pipe 
Class  V Concrete  Pipe 

D.  Imported  Earth  Backfill 

for  Z ^12  ft. 

for  l >12  ft. 

E.  Imported  Gravel  Backfill 

F.  Sheet  Piling 

for  Z £10  ft. 

for  Z >10  ft. 

G.  Dewatering 

H.  Pavement  and  Utility  Restoration 
27  inch  diameter  pipe  or  smaller 
greater  than  27  inch  diameter 

5.  Inlets  and  Outlets 

includes:  headwall 

two  wingwalls 
loose  rock  rip-rap 

6.  Holding  Ponds 

A.  Embankment  construction 
(local  material) 

includes:  clearing  and  grubbing 

foundation  excavation 
compaction 
final  trim 
grass  seeding 


0.431  D 


0.492 

0.586 


1.15  z0.118 

1.15  20 . 11 8 

1.15  -,0.118 


Cost  in  dollars  per  foot  of  length 
for  trench  depths  greater  than  12  feet 


1 15 
0.181  D1,  3 

1 15 
0.207  D1-  10 

0.246  D1,15 


z0.45 

z0.45 

z0.45 


Cost  in  dollars  per  foot  of  length 
(T§  + 3)  ^ (1-75) 

0.03  D0,78  B 
0.08  D0,78  B 


20 

20  + (.70  Z-43)(Z  - 10) 
0.73 


1.11  D 
0.32  D 


0.23 

0.62 


$40  per  inch  of  pipe  dia. 


$2.90  per  cubic  yard 
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8.  Spillway  (uncontrolled 
Capaci ty 
100  cfs  or  less 
110  to  400  cfs 
450  to  1 ,000  cfs 
over  1 ,000  cfs 


overpour) 

$4,000  each 
$8,000  each 
$12,000  each 
$16,000  each 


C.  Outlet 

includes:  24  inch  pipe  $4,000  each 

inlet 

control  gate 
outlet 


The  program  applies  the  foregoing  unit  costs  to  the  specified 
quantities  and  calculates  a subtotal.  To  this  subtotal  the  program  adds 
50  percent  for  contractor's  profit,  engineering,  legal,  and  contingencies 
to  arrive  at  a total  capital  cost  of  construction. 


All  unit  costs,  except  land  costs,  are  based  upon  the  ENR 
construction  cost  index  of  1450.  Adjustments  to  past  or  future  years 
can  be  made  using  past  or  projected  cost  indices  as  summarized  below. 

Date  EUR  Construction  Cost  Index 


1913 

100 

Rational 

Dec.  1950 

560 

Seattle 

Dec.  1955 

690 

Seattle 

Dec.  1960 

850 

Seattle 

Dec.  1965 

1030 

Seattle 

Dec.  1970 

1410 

Seattle 

Mar.  1971 

1450 

Seattle 

June  1973 

1760 

Seattle 

Dec.  1973 

1840 

Seattle 

Land  costs  represent  an  average  of  costs  throughout  the  study 
area  to  purchase  the  land,  plus  50  percent  for  severance  and  acquisition. 
These  are  based  upon  actual  fall  1973  costs^'  within  the  Green  and  Cedar 
River  Basins.  The  unit  land  costs  included  in  the  program  are  the  following 

1.  Data  from  METRO'S  right-of-way  records. 
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Agricul tural , vacant 
Single  family  residential 
Multiple  family  residential 
Commercial  and  industrial 


SO. 12  per  square  foot 
1 .05  per  square  foot 
1 .65  per  square  foot 
2.25  per  square  foot 


It  should  be  noted  that  these  costs  are  not  based  upon  the  ENR 
construction  cost  index,  and  will  not  increase  or  decrease  in  the  same 
manner  as  that  index.  If  new  land  costs  are  deemed  necessary  the  computer 
program  will  need  appropriate  modification. 


DESCRIPTION  OF  THE  COMPUTER  PROGRAM 

This  section  describes  the  logic  and  input  preparation  fo>"  the 
COST  program.  A functional  description  of  the  various  computer  routines 
plus  a detailed  explanation  of  the  data  input  is  provided.  This  explanation 
is  intended  to  provide  the  user  with  a basic  understanding  of  the  program, 
but  is  not  intended  to  provide  the  detail  necessary  to  make  program 
modifications.  If  modifications  are  desired  the  user  should  contact  the 
authors  of  this  report.  A program  listing  is  supplied  at  the  back  of  this 
document. 

PROGRAM  LANGUAGE 

This  model  was  originally  programmed  for  the  UNIVAC  1108  in 
FORTRAN  V.  This  version  of  the  FORTRAN  compiler  is  essentially  compatible 
with  IBM  FORTRAN  level  G and  versions  of  the  model  have  been  executed  on 
IBM  equipment.  This  FORTRAN  is  also  compatible  with  the  Extended  FORTRAN 
compiler  used  on  CDC  6000  series  equipment. 

OPERATING  REQUIREMENTS 


This  program  is  executed  as  an  in-core  unit,  and  has  the  following 
requi rements: 


High  speed  core 


* 14,200g  words 

- 6 , 000 1 o words 

* 24,000]q  bytes 

1 card  reader  (logical  unit  5) 

1 line  printer  (logical  unit  6) 


PROGRAM  LOGIC 
Program  Cost 

The  flow  of  logic  for  the  cost  program  is  indicated  in 

Figure  VII-2.  Program  execution  is  initiated  in  the  executive  routine, 

COST,  which  reads  and  prints  an  overall  run  comment  card,  the  specified 

ENR  index,  and  then  information  which  directs  control  to  each  of  the 

four  component  cost  subroutines.  A grand  total  for  all  facilities  will 
be  accumulated  for  the  cost  run  until  a card  is  encountered  which  directs 
the  program  to  print  the  total  and  reset  the  grand  total  to  zero.  This 
feature  enables  the  user  to  make  several  grand  total  calculations  for 
different  problems  in  the  same  computer  run. 

As  far  as  the  COST  program  is  concerned  the  order  of  calculation 
is  unimportant.  The  program  simply  accepts  the  input  data  in  the  order 
presented  and  then  processes  it  through  the  appropriate  subroutine  in  a 
serial  fashion.  Any  subroutine  can  be  called  any  number  of  times. 

Subroutine  CHANNEL 

Subroutine  CHANNEL  computes  the  costs  associated  with  open 
channel  construction.  The  routine  first  computes  the  total  excavation 
and  land  quantities  and  then  applies  the  unit  cost  factors.  Concrete  or 
rip-rap  lining  can  be  included  in  the  cost  calculation  on  an  optional 
basis.  A summary  of  costs  is  printed  each  time  CHANNEL  is  called.  There 
is  no  limit  on  how  many  channels  may  be  included  in  one  computer  run. 
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PIPES 


Subroutine  PIPES  computes  the  installation  cost  of  pipe  excava- 
tion and  placement  from  the  given  unit  costs.  The  total  capital  costs  is 
obtained  from  the  sum  of  the  costs  of  the  pipe,  excavation,  installation, 
gravel  backfill  to  a depth  of  six  inches  above  the  top  of  the  pipe,  catch 
basins  or  manholes,  and  cleanup.  Options  are  imported  earth  backfill, 
imported  gravel  backfill,  sheet  piling,  dewatering,  and  pavement  and  utility 
restoration.  No  land  costs  are  included  as  pipes  are  assumed  to  be  located 
in  public  ri ghts-of-way . A maximum  of  100  pipes  may  be  costed  in  any  com- 
puter run. 


. ■:  routine  OUTLET 

This  subroutine  computes  the  total  capital  cost  of  a headwall, 
two  wingwalls,  and  rock  rip-rap  for  an  inlet  or  an  outlet.  No  land  costs 
are  included  as  inlets  and  outlets  to  pipes  are  assumed  to  be  located  in 
public  rights-of-way.  There  is  no  limit  to  the  number  of  inlets  or  outlets 
which  may  be  included  in  any  computer  run. 

Sui  r .it  ine  PONDS 


The  capital  cost  of  holding  ponds  is  calculated  in  this  subroutine. 
The  cost  calculation  includes  the  cost  of  the  embankment,  an  uncontrolled 
spillway,  an  outlet,  and  land  cost.  Any  number  of  ponds  may  be  input  to  a 
single  computer  run. 

CARD  l HP UT  SPECIFICATIONS 


This  section  outlines  the  exact  specifications  for  defining 
card  input  to  the  COST  program.  In  the  job  stream  cards  (1)  and  (2)  below, 
must  appear  in  each  run.  Subsequent  to  that  the  specifications  for  ponds, 
channels,  inlets/outlets  or  pipes  may  be  in  any  order  or  not  included  at 
all.  Whenever  the  work  PRINT,  beginning  in  card  column  6,  is  encountered. 
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the  current  grand  total  for  the  run  is  printed,  set  back  to  zero,  and 
control  sent  to  the  top  of  the  program  to  read  new  entries  on  cards  (1) 
and  (2)  and  to  loop  through  the  entire  sequence  again. 

Input  for  all  Buns 


Card 

Number 

Card 
Col umn 

FORTRAN 

Format 

FORTRAN 

Name 

1 

1-  8 

20A4 

TITLE 

2 

1-10 

F10.0 

ENR 

Input  for  Channels 

Card 

Card 

FORTRAN 

FORTRAN 

Number 

Col umn 

Format 

Name 

1 

1-  5 

15 

N 

7-14 

A4 

IDX 

N 

1-  8 

18 

J 

more 

cards 

9-16 

F8.0 

D1 

17-24 

F8.0 

B1 

25-32 

F8.0 

SI 

33-40 

F8.0 

D2 

41-48 

F8.0 

B2 

49-56 

F8.0 

S2 

57-64 

F8.0 

XL 

Description 
Any  comment 

The  Engineering  News  Record 
Cost  Index  for  this  run 


Description 

The  number  of  channels  to 
be  costed 

The  word  CHANNELS 

The  channel  ID  number 

Depth  of  existing  channel  (ft) 

Bottom  width  of  existing  channel  (ft) 

Side  slope  of  existing  channel  (H/V) 

Depth  of  new  channel  (ft) 

Bottom  width  of  proposed  channel  (ft) 

Side  slope  of  new  channel  (H/V) 

Channel  length  (ft) 
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Card 

Number 

Card 
Col umn 

FORTRAN 

Format 

FORTRAN 

Name 

Description 

1 

1-  5 

15 

N 

The  number  of  holding  ponds 
to  be  costed 

N 

groups  of 
cards  A-B 

7-14 

A4 

IDX 

The  word  PONDS 

A 

1-  5 

15 

J 

Pond  ID  number 

6-15 

FI  0.0 

SVOL 

Pond  volume  (acre-feet) 

16-25 

F10.0 

AREA 

Pond  area  (acres) 

26-35 

F10.0 

SPILL 

Spillway  capacity  (cfs) 

36-45 

no 

ICLS 

Land  class: 

1 = Ag.  or  vacant 

2 = Single  family  residential 

3 = Multi-family  residential 

4 - Commercial  and  Industrial 

B 

1-  5 

15 

repeat  IR 
group 

Embankment  section  ID  number 

6-10 

F5.0 

1 to  14  H 
times 

Embankment  section  height  (ft) 

• ■ f ■< 

11-20 

Inlet  3/0: 

F10.0 

it  lets 

on  each  XL 
card 

Embankment  section  length  (ft) 

Card 

Card 

FORTRAN 

FORTRAN 

Number 

Col umn 

Format 

Name 

Description 

1 

1-  5 

15 

N 

The  number  of  INLETS/OUTLETS 
to  be  costed 

7-14 

A4 

IDX 

The  word  OUTLET 
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Card 

Number 

Card 
Col  umn 

FORTRAN 

Format 

FORTRAN 

Name 

Description 

65-72 

18 

ICLS 

Land  category: 

1 = Ag.  or  vacant 

2 = Single  family  residenti 

3 = Multi-family  residentia 

4 = Commercial  and  Industri 

73-80 

Input  for  Pipes 

18 

ILYN 

Channel  lining: 

0 = no  lining 

1 = concrete  lining 

2 = rip-rap  lining 

Card 

Card 

FORTRAN 

FORTRAN 

Number 

Col  umn 

Format 

Name 

Description 

1 

1-  5 

15 

N 

The  number  of  pipes  to  be 
costed  (<  100) 

7-14 

A4 

1DX 

The  word  PIPES 

N 

more 

1-10 

110 

NO 

The  pipe  ID  number 

cards 

11-20 

FI  0.0 

DIA 

Pipe  diameter  (in) 

21-30 

FI  0.  0 

LENGTH 

Pipe  length  (feet) 

31-40 

F10.0 

DEPTH 

Average  depth  to  invert  (ft) 

41-45 

15 

INDEX 

Enter  1 for  imported  backfill; 
blank  otherwise 

46-50 

15 

INDEX 

Enter  1 for  gravel  backfill; 
blank  otherwise 

51-55 

15 

INDEX 

Enter  1 for  sheet  piling, 
blank  otherwise 

56-60 

15 

INDEX 

INDEX 

Enter  1 for  dewatering, 
blank  otherwise 

Enter  1 for  restoration, 
blank  otherwise 
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EXAMPLE  CALCULATIONS 


An  example  calculation  for  each  type  of  cost  calculation,  pipes, 
ponds,  channels  and  inlets/outlets  is  included  in  this  section.  The  card 
data  used  in  the  examples  is  listed  and  shown  at  Table  V 1 1 - 1 . Note  that 
all  the  data  is  entered  in  one  computer  run  and  the  PRINT  card  separates 
the  calculation  of  the  separate  items.  If  the  PRINT  card  was  omitted 
the  comment  and  ENR  index  cards  would  be  deleted  for  all  cases  except 
the  first,  and  the  run  results  would  be  essentially  the  same.  The  PRINT 
card  was  included  here  to  allow  a different  comment  to  be  printed  for 
each  example  output. 

Table  VII-2  shows  the  costs  associated  with  the  enlargement  of 
two  open  channels.  No  channel  lining  costs  have  been  included  in  this 
example.  Table  V 1 1 -3  shows  the  results  for  the  construction  of  an  inlet 
on  a 24-inch  pipe  and  an  outlet  on  a 48-inch  pipe.  There  are  no  options 
in  the  INLET/OUTLET  cost  calculation.  The  costs  associated  with  two 
holding  ponds  are  displayed  in  Table  VII-4.  As  can  be  seen  the  first 
pond  has  two  different  embankment  sections  while  the  second  is  defined 
by  three.  The  final  example.  Table  VII-5,  shows  the  cost  of  placing  two 
pipes,  one  42-inch  and  one  30-inch.  This  example  includes  gravel  backfill 
and  the  option  for  including  pavement  and  utility  restoration. 

These  four  examples.  Tables  VII-2  to  VII-5,  cover  the  complete 
range  of  the  COST  program's  capability  and  should  be  thoroughly  understood 
by  any  potential  user  of  the  program. 
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EXAMPLE  COST  CALCULATION  FOR  PIPES 
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CHAPTER  VIII 
PROGRAM  LISTINGS 


RUNOFF  BLOCK 


The  runoff  block  program  listing  is  presented 
on  the  following  sheets  identified  as  Pages  3 
through  37. 
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The  cost  determination  program  listing  is 
presented  on  the  following  sheets  identified 
as  Pages  1 through  12. 
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